VotumE 17 JANUARY, 1950 


CONTENTS 


Predictions of Supersonic Airplane Performance HAROLD LuSKIN 


The Calculation of Supersonic Downwash Using Line Vortex Theory 
HAROLD MIRELS AND RUDOLPH C. HAEFELI 


The Pulse Method for the Determination of Aircraft Dynamic Performance 
R. C. SEAMANS, JR., B. P. BLASINGAME, AND G. C. CLEMENTSON 


The Theoretical Lateral-Stability Derivatives for Wings at Supersonic Speeds 
ARTHUR L. JONES 


Influence of Acceleration on Aerodynamic Characteristics of Thin Airfoils in Super- 
sonic and Transonic Flight C. S. GARDNER AND H. F. LupLorr 


Readers’ Forum: ‘Supersonic Airfoils Simplified,’’ by James E. Broadwell; ““Un- 
steady Flow Theory in Dynamic Stability,” by John W. Miles; “Stability of 
Sandwich Plates in Combined Shear and Compression,” by P. P. Bijlaard; “On 
the Reduction of Unsteady Supersonic Flow Problems to Steady Flow Problems,” 
by John W. Miles 


Published by the 
INSTITUTE OF THE AERONAUTICAL SCIENCES, INC. 


— 
BRIAT 

ING Y 

4 

. 
; 

13 

22 
47 


JOURNAL OF THE 
AERONAUTICAL SCIENCES 


UsRARY 
/ 

SS b 


Editor: Hucu L. Drypen, Director, N.A.C.A. 
Associate Editor: Ropert R. DEXTER Associate Editor: Berneice H. Jarcx 


The Institute of the Aeronautical Sciences publishes monthly the JourNAL oF THE 
AERONAUTICAL SCIENCES AND THE AERONAUTICAL ENGINEERING REVIEW. 


The JourNAL OF THE AERONAUTICAL SCIENCES prints papers of scientific and technical 
interest to Institute members and subscribers. In general, it provides an opportunity 
for specialists to publish articles disclosing new knowledge and new applications in the 
field of the aeronautical sciences. The Journal restricts itself to the aeronautical as- 
pects of advances in science and engineering, and papers printed have either been sub- 
mitted directly to the Journat or have been previously presented at a meeting of the 
Institute. 

The AERONAUTICAL ENGINEERING REviEw keeps members and subscribers informed 
regarding the contents of new books and current periodicals in the field of aeronautics. 
It reviews new developments in the aircraft industry and reports on progress from house 
organs and catalogs published by aircraft and engine manufacturers and other com- 
panies engaged in aeronautical activities. It publishes papers on applied engineering 
and aircraft production. News of Institute activities is published in the Review. 


Statements and opinions expressed in the JourNAtL are to be understood as individual 
expressions of the authors and not those of the Institute. 


SUBSCRIPTION RATES 


Journal of the Aeronautical Sciences 


Published Monthly 
Foreign Countries Including Canada 


United States and Possessions (American Currency Rate) 


Aeronautical Engineering Review 
Published Monthly 


Foreign Countries Including Canada 


United States and Possessions (American Currency Rates) 


1 

Single 0.50 Single 0.50 
Notices of change of address should be sent to the Institute at least 30 days prior to 
actual change of address. 
Manuscripts for publication, proofs, and all correspondence should be addressed to the 
Editorial Office of the JouRNAL. 
Correspondence regarding subscriptions may be addressed to Publication Office, 20th 
and Northampton Sts., Easton, Pa., or to the Editorial Office of the JouRNAL OF THE 
AERONAUTICAL SciENcES, 2 East 64th Street, New York 21, New York. 

Copyright, 1950, by the Institute of the Aeronautical Sciences, Inc. 


Entered as Second Class Matter at the Post Office, Easton, Pa., May 1, 1937. Acceptance for mailing at a special rate of 
Postage provided for in the Act of August 24, 1912. Authorized April 29, 1937. 


. 
| 
hee) 
a 
; 
4 
: 
ag 
4 
yy 
J 
4 
| 
Were 
ae 


President 
w. A. M. BURDEN 
Vice-Presidents 
J. L. ATwooD 


J. McCaRTHY 
G. S, SCHAIRER 
Eimer A. SPERRY, JR. 
Director 
§ PAUL JOHNSTON 
Assistant Director 
James L. STRAIGHT 


Pp. R. BASSETT 

Rex B. BEISEL 

WiLLIAM BOLLAY 

VictoR E. CARBONARA 
CuarLes H. CoLvin 
SHERMAN M. FAIRCHILD 


WALTER BEECH 
FraNK W. CALDWELL 
S. DAMON 


R. H. FLEET 
R. PAauL HARRINGTON 


Carvin M. Bo_steR 
Prank O. CARROLL 

B. U. ConpoNn 

Harry F. GUGGENHEIM 
Louis A. JouNSON 


H. H. ARNOLD 
Prank W. CALDWELL 
Six Grorrrey HAVILLAND 
ames H. DooLitTTLe 

WALD W. DouGLas 
L. DrypEN 


. ALDRIN 
Hart O. Berc (Deceased) 
W. A. M. BurDEN 


Donatp W. Douctas 
Autan D. Emit 
Suerman M. Farrcuiip 


Ina H. Aspott 
Joun D. AKERMAN 

B. E. ALpRIN 

H. Juttan ALLEN 
Harry G. ARMSTRONG 
Kart ARNSTEIN 

L. Atwoop 

au. S. BAKER 

S. BARNABY 
P.R. Bassett 
Wetiwoop E. Beat. 
Rex B. 

G. M. BeLLanca 

Don R. BERLIN 

Henry A. BaRLINER 
Maurice A. Brot 
Wittram 
Grorce W. Brapy 
Lyman J. 

W. G. BromBacHerR 
G. Brown 
P. BURGESS 
Kewnetn 
Cuartes H. CHATFIELD 
ALLAN CHILTON 

Francis CLAuseR 
FRANKLIN R. 
Cartes H. Corvin 

RL T. Compton 
Emerson W. Conon 
Onn W. CrRowLey, JR. 
PH S. DAMON 

Luis pz FLorez 
Suita J. pz FRANCE 
J.P. Den Hartoo 
Water S. 

S. DRAPER 
Ivan H. Driccs 
Wrram K. Exe. 


INSTITUTE 


KENNETH S. M. DAVIDSON 


of the AERONAUTICAL 


Past-Presidents 


J. C. HUNSAKER 


Cuarves L. LAWRENCE 
DonaLp W. DouGLas 


GLENN L. MarRTIN 


R. 


T. P Wricut 


Georce W. Lewis (Deceased) 
James H. Doo.ittLe 
FRANK W. CALDWELL 


L. 
L. Drypen 
H. Freer 


Cuarves H. Corvin 
ARTHUR E. RAYMOND 


P. R. Bassett 
Jonn K. Norturope 


COUNCIL 


LuIS DE FLOREZ 

R. G. Fo_som 

C. C. FurNas 
Leroy R. GRUMMAN 
M. HAwKINS 


R. M. Hazen 
HALL L. HIBBARD 

J. C. HUNSAKER 

TH. VON KARMAN 
Tuomas A. KNOWLES 
WiLuraM LITTLEWoop 


ADVISORY COMMITTEE 


J. A. HERLIHY 
WILLSON H. HUNTER 
Joun C. LESLIE 

J. S. McDONNELL 
Joun K. NorRTHROP 


GROVER LOENING 
GLENN L. MARTIN 
CLARK B. MILLKAN 
Tuomas A. MorGAN 
ARTHUR E. RAYMOND 
F. W. REICHELDERFER 


HONORARY MEMBERS 


Marion M. Katrez 

. A. KimBati 
THEODORE LONNQUEST 
Joseru J. O'CONNELL, JR. 


W. F. Duranp 
Str RIckarD Fairey 
Sir Roy Feppen 


Sir FrRepgeRIcK HaNDLey Pace 


J. C Hunsaker 


RicHarp FarrRey 
R. H. Freet 
Lester D. GARDNER 


Florence GuGGENHEIM (Deceased) 


C. D. Hanscom 
Joun P. V. HEINMULLER 


Ovtver P. 

C. L. Ectvept 
FREDERIC FLADER 
CHARLES FROESCH 

Jack Frye 

GARLAND FULTON 

C. C. Furnas 

Lester D. GARDNER 
Ricnarp C. GazL_ey 
E. GLunarerr 
Me N. Goucu 
Leroy R. GruMMAN 
Harorp R. Harris 

R. M. Hazen 

Epwarp H. HeIngMANN 
S. D. Heron 

Hatt L. 

L. S. Hopss 

Paut E. Hovcaarp 
HovGaarp 
RoBerRT INSLEY 
EasTMan N. Jacosps 

C. L. JoHnson 

S. JOHNSTON 
KaPLan 
ALEXANDER KARTVBLI 
Paut H. KeMMER 
CARLETON KEMPER 

J. H. KinDELBERGER 
Otro E. KircHNER 
ARTHUR L. 
ALEXANDER KLEMIN 

W. B. KLEMPERER 
Paut KOLLSMAN 

Orto C. Kopren 

B. V. Korven-Krovukovsky 
S. M. Kraus 

ARNOLD M. Kurtue 

I. M. Lappon 


HONORARY FELLOWS 


BENEFACTORS 


AMERICAN FELLOWS 


Howarp T ORVILLE 
Georce B. Pecram 
ALFRED M PRIDE 
Donatp L Putt 


Sir Jones 
Gienn L. Martin 
CrarRK MILLIKAN 

J. Laurence PritcHaRD 
D. R. Pye 


ALEXANDER KLEMIN 
Paut Ko_usMan 
C. LANDAUER 
LANDAUER 
GROVER LOENING 


Epmunp C. Lyncu (Deceased) 


E. S. Lanp 

CuHarves L. Lawrence 
Jerome LEDERER 
Joun G. Lee 

Joun C. 
CuHarves A. LINDBERGH 
WituraM LitTLEWoop 
GROvER LOENING 
ALBERT E. LomBARD 
F. A. Loupen 

W. Lovecaces, II 
EuGenkt E. Lunpguist 
R. D. MacCart 

Joun R. MarKHam 
ERLE MARTIN 

Wm. H. McAvoy 

C. J. McCartuy 

J. S. McDonngeLL 
Ropert J. 
Joseru S. 
S. NILes 
Joun K. Nortrurop 
ARTHUR Nutr 
SHATSWELL OBBR 
RoBErRT R. OSBORN 
W. Bartey OswaLp 
Grorce A. Pace, Jr. 
F. W. PAWLOwsSKI 

F. W. Pennoyer, JR. 
A. A. PRIESTER 
ARTHUR E. RayMonpD 
F. W. REICHELDERFER 
G. Rep 

H., J. E. Rei 

H. C. RicHaRDsON 

L. B. RICHARDSON 
Russet G. ROBINSON 
LeonaRD E. Root 


SCIENCES 


Secretary 
Rosert R. DEXTER 
Treasurer 
D. OsBoRN 
Controller 
JoserH J. MAITAN 
Counsel 


ALLAN D. EmIL 


L. B. RICHARDSON 
RussE__ G. ROBINSON 
D. R. SHOULTS 
ERNEstT G. STouT 
GEorGE H. TWENEY 
H. Russe_t YouNG 


W. C. ROCKEFELLER 
T. CLAUDE RYAN 
EDWARD P. WARNER 
BURDETTE S. WRIGHT 
T. P. Wricut 


F. W. ReicHeLperrer 
W. Rentzer 

W. Stuart SyMINGTON 
Hoyt S VaNDENBERO 


Icor I. Srkorsky 

Sir Georrrey TAyLor 
Tu. von KARMAN 
EDWARD P. WARNER 
T. P. Wricut 


F. W. Macin 
GLENN L. MARTIN 


Grorce J. Meap (Deceased) 


S. A. Reep (Deceased) 


LauRANce ROCKEFELLER 


Epwarp Warner 


C. E. RosENDAHL 
Appison M. RorHrock 
W. R. Sears 

Gacen B. SCHUBAUER 
WiiuraM T. SCHWENDLER 
Francis R. SHANLEY 
James M. SHOEMAKER 
Icor ALExIs SrKORSKY 
ABE SILVERSTEIN 

C. RIcHARD SODERBERG 
SELDEN B. SPANGLER 
Ev_Mer A. Sperry, JR. 
Joun Stack 

RoBERT M. STANLEY 
B. Stout 
Harry A. Sutton 

C. Fayette TayLor 
Epwarp S. Taytor 
B. Tayior 

TH. THEODORSEN 

T. E. Trcurncuast 

S. TrmosHENKO 

F. L. THoMPsoN 

L. B. TUCKERMAN 
Raven H. Upson 

J. ParKerR VAN ZANDT 
J. G. Vincent 
RICHARD VON MISES 
Frank L. WATTENDORF 
P. V. H. Weems 

Frep E. Weick 

R. D. WEYERBACHER 
Joun B. WHEATLEY 
DonaLp H. Woop 
Ropert J. Woops 
Joun E, YounGER 
ALBERT F. ZAHM 
J. ZaAND 


LAURANCE S. ROCKEFELLER 


‘ 
= 
it 
| 
- 
= 
AS ~ 


Corporate Members 
of the 
Instetute of the Aeronautical Sciences 


Academy of Aeronautics, Inc. 
Casey Jones School of Aeronautics 


Alvin P. Adams and Associates 
Adams Rite Manufacturing Company 
Aerojet Engineering Corporation 
Aerolab Development Company 
Aeroproducts Division, General Motors Corpora- 
tion 
Aeroquip Corporation 
Agawam Aircraft Products, Inc. 
Airborne Accessories Corporation 
Airborne Instruments Laboratory, Inc. 
Aircraft Radio Corporation 


AiResearch Manufacturing Company, Division of 

The Garrett Corporation 
Allis-Chalmers Manufacturing Company 
Allison Division, General Motors Corporation 
Aluminum Company of America 
American Airlines System 
American Bosch Corporation 
American Overseas Airlines, Inc. 

Iceland Airport Corporation (Subsidiary) 
American Phenolic Corporation 
Associated Aviation Underwriters 
Atlas Supply Company 
Baker Steel & Tube Company 
Beech Aircraft Corporation 
Bell Aircraft Corporation 
Bendix Aviation Corporation 

Bendix International Division 

Bendix Products Division 

Bendix Radio Corporation 

Eclipse-Pioneer Division 

Friez Instrument Division 

Pacific Division 

Scintilla Magneto Division 

The BG Corporation 
Boeing Airplane Company 
Seattle Division 
Wichita Division 
Breeze Corporations, Inc. 
Aircraft Standard Parts Company, Inc. 
Aldrich Company 
Anderson Stove Company, Inc. 
Federal Laboratories, Inc. 
Foundry Service, Inc. 


Brooks and Perkins, Inc. 


Cal-Aero Technical Institute, A Division of Gran; 
Central Airport Company 


Canadair, Ltd. 

Cessna Aircraft Company 

Chandler-Evans Division, Niles-Bement-Pon 
Company 

Chase Aircraft Company, Inc. 

Chase National Bank of the City of New York 

The Cleveland Pneumatic Tool Company 
Automatic-Aircraft Division 

Clifford Manufacturing Company, Division of 
Standard-Thomson Corporation 


Consolidated Vultee Aircraft Corporation 
Fort Worth Division 
Stinson Aircraft Division 

Continental Motors Corporation 


Curtiss-Wright Corporation 
Airplane Division 
Propeller Division 
Wright Aeronautical Corporation Division 
Dayton Aircraft Products, Inc. 
Doak Aircraft Company, Inc. 
Doelcam Corporation 


Douglas Aircraft Company, Inc. 
El Segundo Plant 
Long Beach Plant 
Santa Monica Plant 

The Dow Chemical Company 

Dzus Fastener Company, Inc. 

Eastern Air Lines, Inc. 

Eaton Manufacturing Company 
Dynamatic Corporation (Subsidiary) 
Eaton-Wilcox-Rich Ltd. (Subsidiary) 

Edo Corporation 

Elastic Stop Nut Corporation of America 

Electrol Incorporated 

Esso Standard Oil Company 

Ethyl Corporation 


Fairchild Camera & Instrument Corporation 
Fairchild Aerial Surveys, Inc. 


Fairchild Engine and Airplane Corporation 
Al-Fin Division 
Fairchild Aircraft Division 
Fairchild Personal Planes Division 
Fairchild Pilotless Plane Division 
NEPA Division 
Ranger Aircraft Engines Division 
Stratos Corporation (Subsidiary) 


federal 
Feder 
Feder 
(Pz 
The Fit 
forni: 
Fletche 
Genera! 
Airct 
Avia’ 
Globe 
The B. 


| The Ge 


Goo 
W. &I 
Harvey 
Industr 
Insurat 

The 


Melet: 
Minne 

Bro 
Natio’ 
North 

Not 
North 
North 
Pan A 
The P 
Pesco 

tior 


. 
i 
™ 
Insut 
Phil: 
The In 
Irving 
Jack & 
Farle } 
Walter 
Kollsn 
Lavell 
Lear Ii 
‘ 
‘ Link f 
Lockh 
Ai 
irq 
Loc! 
Loc! 
Longii 
Marm 
ott 
roo: 
Rot 


of Grand 


1ent-Pond 


v York 


‘ision of 


sion 


ration 


federal Telephone and Radio Corporation 
Federal Electric Manufacturing Company, Ltd. 
Federal Telecommunication Laboratories, Inc. 

(Partial Subsidiary) 

The Firestone Tire & Rubber Company of Cali- 
fornia 

Fletcher Aviation Corporation 

General Electric Company 
Aircraft Gas Turbine Division 
Aviation Divisions 

Globe Corporation, Aircraft Division 

The B. F. Goodrich Company 

The Goodyear Tire & Rubber Company 
Goodyear Aircraft Corporation 

Grumman Aircraft Engineering Corporation 

W. & L. E. Gurley 

Harvey Machine Company, Inc. 

Industrial Sound Control 

Insurance Company of North America Companies 
The Alliance Insurance Company of Philadel- 


hia 
eke Insurance Company of North Amer- 
ica 
Insurance Company of North America 
Philadelphia Fire and Marine Insurance Com- 
pany 
The International Nickel Company, Inc. 
Irving Air Chute Company, Inc. 
Jack & Heintz Precision Industries, Inc. 
Johns-Manville Sales Corporation 
Earle M. Jorgensen Company 
Walter Kidde & Company, Inc. 
Kollsman Instrument Division, Square D Com- 
an 
Lavelle Aircraft Corporation 
Lear Incorporated 
Link Aviation, Inc. 
The Liquidometer Corporation 
Lockheed Aircraft Corporation 
Airquipment Company (Subsidiary) 
Lockheed Air Terminal (Subsidiary) 
Lockheed Aircraft Overseas Corporation 
(Subsidiary) 
Lockheed Aircraft Service, Inc. (Subsidiary) 
Longines- Wittnauer Watch Company, Inc. 
Marman Products Company, Inc. 
The Glenn L. Martin Company 
GLM Chemical Division 
Rotawings Division 
McDonnell Aircraft Corporation 
Meletron Corporation 
Minneapolis-Honeywell Regulator Company 
Brown Instrument Company (Subsidiary) 
National Credit Office, Inc. 
North American Aviation, Inc. 
North American Aviation, Inc., of Texas 
Northrop Aircraft, Inc. 
Northwest Airlines, Inc. 
Pan American World Airways System 
The Parker Appliance Company 
Pesco Products Division, Borg-Warner Corpora- 
tion 


Phillips Petroleum Company 
Pure Oil Company 
Radioplane Company 
Republic Aviation Corporation 
Rheem Manufacturing Company 
E. V. Roberts and Associates 
A. V. Roe Canada Limited 
John A. Roebling’s Sons Company 
Rohr Aircraft Corporation 
Ryan Aeronautical Company 
Saval, Inc. 
Schrillo Aero Tool Engineering Company 
Servomechanisms, Inc. 
Shell Oil Company, Inc. 
Simmonds Aerocessories, Inc. 
A. O. Smith Electrical Manufacturing Company 
Socony-Vacuum Oil Company, Inc. 
Solar Aircraft Company 
Sperry Gyroscope Company Division of The 
Sperry Corporation 
Square D Company 
Standard Oil Company of California 
Standard Oil Company (Indiana) 
The Steel Products Engineering Company 
Stewart-Warner Corporation, South Wind Divi- 
sion 
Sturgess, Inc. 
Summers Gyroscope Company 
Telecomputing Corporation 
H. I. Thompson Company 
Thompson Products, Inc. 
Tinnerman Products, Inc. 
Toolko Engineering Company 
Union Carbide and Carbon Corporation 
Bakelite Corporation 
Electro Metallurgical Company 
Haynes Stellite Company 
The Linde Air Products Company 
National Carbon Company, Inc. 
United Air Lines, Inc. 
United Aircraft Corporation 
Chance Vought Aircraft Division 
Hamilton Standard Propellers Division 
Pratt & Whitney Aircraft Division 
Sikorsky Aircraft Division 
United Aircraft Export Corporation 
United Aircraft Service Corporation 
United States Aviation Underwriters, Inc. 
Vard Inc. 
Vickers, Inc. 
Western Gear Works 
Westinghouse Electric Corporation 
Aviation Gas Turbine Division 
Industrial Electronics Division 
Marine & Aviation Sales Department 
Small Motor Division, Aviation Sales 
Weston Electrical Instrument Corporation 
Wyman-Gordon Company 
Young Radiator Company 


Aerodynamics 
Francis Clauser 
Howard Emmons 
R. Paul Harrington 
R. T. Jones 

Carl Kaplan 

W. B. Klemperer 
Otto C. Koppen 
John G. Lee 

H. W. Liepmann 
Cc. C. Lin 

H. F. Ludloff 
John R. Markham 
Clark B. Millikan 
Shatswell Ober 
W. Bailey Oswald 
Allen E. Puckett 
Elliott G. Reid 

H. J. F. Reid 
Russell Robinson 
G. B. Schubauer 
W. R. Sears 

John Stack . 

Th. von K4rm4n 
Richard von Mises 


Air Transport 
Charles Froesch 
J. A. Herlihy 
Otto E. Kirchner 
John C. Leslie 


Aircraft Design 
W. E. Beall 

Rex B. Beisel 
W. K. Ebel 

Hall L. Hibbard 
C. L. Johnson 

I. M. Laddon 


Aircraft Industries 
Marcus 


EDITORJAL COMMITTEE 


HuGu L. Drypen, Chairman 


C. J. McCarthy 
Fred E. Weick 
Robert J. Woods 


Physiologic Problems 


H. G. Armstrong 

Louis H. Bauer 

Eugene Du Bois 

W. Randolph Lovelace, II 
Ross A. MacFarland 

A. D. Tuttle 


Aircraft Propulsion 


George W. Brady 
Frank W. Caldwell 
L. S. Hobbs 

E. E. Johnson 
Joseph Keenan 
R. P. Kroon 
Arthur Nutt 
Hans Reissner 
Abe Silverstein 

C. Fayette Taylor 
E. S. Taylor 

P. Taylor 


Structures and Materials 


W. B. Bleakney 
E. W. Conlon 
J. P. Den Hartog 


E. undquist 


C. F. Nagel, Jr. 
J. S. Newell 
Alfred S. Niles 
W. Ramberg 
E. Reissner 

F. R. Shanley 
H. W. Sibert 
C. R. Strang 
S. Timoshenko 
H. S. Tsien 
Chi-Teh Wang 


Meteorology 


C. F. Brooks 
D. M. Little 
E. J. Minser 
F. W. Reichelderfer 
H. C. Willett 


Fuels and Oils 


A. L. Beall 

J. H. Doolittle 
Graham Edgar 
R. T. Goodwin 
S. D. Heron 

A. M. Rothrock 


Instruments 


Allan G. Binnie 
W. G. Brombacher 
Charles H. Colvin 
B. M. Craig 

B. Del Mar 

C. S. Draper 

O. E. Esval 

W. L. Howland 


Association of America, Inc.—CHARLES 


American Association for the Advancement of Science—J. C. 


HUNSAKER 


American Society of Mechanical H. CoLvin 


American Standards Association—P. R. BASSETT 


Engineering Society of Detroit—Donatp C. Hunt, GEORGE 


TWENEY 


Greater New York Safety LitTLEwoop 


Daniel Guggenheim Medal Board of Award —PRESTON R. BASSETT, 


J. E. Lindberg 
David W. Moore 
George H. Purcell 
W. A. Reichel 

C. F. Savage 

O. Hugo Schuck 
C. L. Seward, Jr. 
R. Sylvander 

D. K. Warner 

P. F. Weber 


Rotating Wing Aircraft 
Robert Coleman 

K. H. Hohenemser 
Alexander Klemin 

Ralph McClarren 

R. H. Miller 

R. H. Prewitt 

Igor I. Sikorsky 

Robert A. Wolf 


Vibration and Flutter 


Lee Arnold 

M. A. Biot 
William Bollay 
Chieh-Chien Chang 
Martin Goland 

J. P. Den Hartog 
Manfred Rauscher 


Flight Testing 


W. K. Ebel 

R. R. Gilruth 
Melvin N. Gough 
Robert Stanley 
E. G. Stout 


1.A.S. REPRESENTATIVES TO OTHER ORGANIZATIONS 


W. A. M. Burpen, S. PAUL JOHNSTON 
Institute of Radio Engineers—L. M. Hutu 
Manufacturers Aircraft Association—ARTHUR E. RAYMOND 


National Aeronautic Association—WILL1IAM R. ENYART 


National Fire Protection Association—E. E. ALDRIN 


National Research 
DRYDEN 


National Safety Council—JEROME LEDERER 
Society of Automotive Engineers—STEPHEN J. ZAND 


Council Engineering Division—Hucu L. 


VOL 


| 


Prev! 
airplan 
as bein; 
type of 
tors in 
for sup 
It is po 
safety 
creasin 


, 


1, sho" 
the U1 
is seer 
centas 
and, fi 
charac 
dustry 
declin 
stand 
break 
natur 
figure 
which 
move 

Pro 
The a 
sever: 
growt 
tion t 
level. 
brilliz 
point 
chang 


Pres 
Meeti 


. 
> 
ye at 
L. H. Donnell 
E. C. Hartmann 
N. J. Hoff 
a 
# 
2 
¢ 
: 
§ 


raft 


SETT, 


JOURNAL THE 
AERONAUTICAL SCIENCES 


VOLUME. 17 


JANUARY, 1950 


NUMBER 1 


ABSTRACT 


Previous writers on the history and prediction of the trend in 
airplane maximum speeds have sometimes represented the trend 
as being linear with time. This paper discusses the more general 
type of growth curve in relation to airplane speeds. Limiting fac- 
tors in airplane performance are studied, and a prediction is made 
for supersonic developments expected to be possible by 1960. 
It is pointed out that increases in speed and improvements in the 
safety record have been obtained simultaneously, in spite of in- 
creasing wing loadings. This trend is expected to continue. 


THE GrRowTH CURVE 


Soe NATURE OF THE GROWTH of things is in many 
cases! given by a curve of a well-known type. Fig. 
1, showing the history of steam railroad road mileage in 
the United States,*» * is an example of such a curve. It 
is seen that progress begins at roughly a constant per- 
centage rate; then tapers off; reaches an upper limit; 
and, finally, begins a decline. In some cases additional 
characteristics are found. Fig. 2, from a competing in- 
dustry, shows a curve of growth and greater eventual 
decline,! the high L/D of horse-drawn vehicles notwith- 
standing! Fig. 3 shows a growth in two segments, the 
break corresponding to a revolutionary change in the 
nature of the growing thing.' The last two of these 
figures are plotted on semilogarithmic coordinates, on 
which the initial curvature to be noted on Fig. 1 is re- 
moved. 

Progress in aeronautics also has its growth curve. 
The aeronautical engineer, however, is likely to spend 
several years working on a design at one fixed level of the 
growth curve. This has the result of dulling the percep- 
tion to the time derivative of the curve existing at that 
level. The failure of past predictions of even our most 
brilliant engineers will later be discussed in proof of this 
point. It is believed to be important that the sense of 
change, the awareness of growth, be developed.’ 


Presented at the Aerodynamics Session, Annual Summer 
Meeting, I.A.S., Los Angeles, July 21-22, 1949. 
* Aerodynamics Research Engineer, Santa Monica Plant. 


Predictions of Supersonic Airplane 
Pertormance 


HAROLD LUSKIN* 
Douglas Aircraft Company, Ine. 


To establish an appreciation for the long-term growth 
trend, a few examples thought to be of special interest 
will be drawn from fields of transportation other than 
the aeronautical. 

The real beginning of our railroad system took place 
in 1830. To provide transportation for the rapidly 
growing nation, there had been, in the decades just 
previous to this, considerable development of travel on 
our rivers and in man-made barge canals. Water travel 
was not, however, a perfect method of transportation, 
and in its heyday “500 people annually lost their lives in 
river boat accidents, and Charles Dickens was able to 
say that the high-pressure river boats ‘conveyed that 
kind of feeling to me which I should be likely to experi- 
ence, I think, if I had lodgings on the first floor of a 
powder mill.’ ’’* Considering the population size, and 
the low speed of the river boats, the fatalities per pas- 
senger-mile for river boats must have been enormous 
compared to those for modern air travel. The gradual 
ascendency of the railroads over river boats and barges 
was not accomplished without difficulty, however. 
An amusing editorial view of the competition offered 
by the upstart rail industry can be gained from the fol- 
lowing item from the Ploughman’s Gazette Courier, Can- 
ton, Ohio, August 15, 1830:T 

“Returning Travellers report that the people of the 
East are still suffering from the Railroad mania. Dis- 
regarding the warnings of Providence, they have ex- 
tended the line of the Railroad thirteen miles to Elliott’s 
Mills west of Baltimore. 

“One traveller recently told of having seen with his 
own eyes a single horse drawing two waggons containing 
41 persons at a speed of 11 miles per hour on the portion 
of track which is completed. He heard, however, of 
several serious accidents caused by the horse stumbling 
and being run over by the heavy waggons, which were 
derailed and upset. 


t Furnished by courtesy of the Association of American Rail- 


roads. 
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“Little confidency is now felt in Mr. Peter Cooper’s 
plan for propelling the carriages by means of a steam en- 
gine. 

“Except as a novelty, however, Railroads will never 
take the place of canals. Railroads are untried for long 
distances in any country, and for short distances they 
are still in the experimental stage. The longest in exist- 
ence is the London and Manchester which is but 40 
miles in length. 

“Any farmer can build a canal boat of 25 tons capac- 
ity from materials growing on his own land. To carry 
this amount on a Railroad would require eight waggons 
and a locomotive costing $4,000 at the least, and no 
farmer could afford this equipment. 

“His canal boat would carry livestock, hay, firewood, 
large trees for ship buildings, boards, planks and grain. 
Railroads can not do this. Imagine a load of hay com- 
ing along a Railroad. The sparks would set it afire be- 
fore it started. 

“The farmer who is watchful of his interest will not 
be misled by talk about Railroads, which would be of no 
service to him and destroy his crops and his stock —but 
will endeavor in all possible ways to promote the build- 
ing of canals.” 

In the first example, we saw that an economically 
valuable means of transportation successfully grew 
through a period of bad safety. In the second, we saw 
lack of vision. 


THE GROWTH OF AIRPLANE SPEEDS 


The 45 years that have elapsed since the first powered 
airplane flight have provided a wealth of statistical 
aviation information for study purposes. Certain of 
these data were used in a recent paper by Clauser,° in 
which an analysis of the history of aircraft power 
plants was made. The history of aircraft maximum 
speeds, measured in formal fashion over a course at sea 
level, has also been frequently discussed, and the data 
are again reproduced in Fig. 4. A typical growth curve 
has been faired through the data on propeller-recipro- 
cating power-plant powered airplanes. Various predic- 
tions®* by skillful engineers are shown, the rule 
being that the predictions are conservative. There is 
no intent here to belittle these predictions. The point 
can be made, on the other hand, that these men were so 
familiar with the difficulties facing designers of the day 
that further rapid advances must have seemed hardly 
likely. The prediction of Heinemann made in 1947 has 
yet to be thoroughly put to test. Developments in the 
next few years’ flights will begin to provide the desired 
comparison. A change in the method of determining 
the speed record is also needed, particularly in allow- 
ing higher altitudes to be used. The maximum speed 
of transport airplanes has followed a regular trend and 
is shown on Fig. 5. 

No attempt will be made to predict the future growth 
curve of airplane maximum speeds in the present paper. 


It is felt that the increases will be more rapid, however, 
than the average rate for the past 45 years and that a 
new segment of the growth curve—that for jet airplanes 
—is beginning. This would produce a speed-growth 
curve similar in shape to that of Fig. 3. The first por- 
tion of such a curve is faired through the data of Fig. 4, 
using the standard form for the growth curve. » There 
are insufficient data available at this time on jet speed 
records to attempt to fit the new portion of the growth 
curve. 

A more fundamental reason for a new trend in the 
growth curve than the empirical data of Fig. 4 is the 
new type of thrust and resistance variations to be ex- 
pected from modern jet units and supersonic air frames. 
The point was made in connection with Fig. 3 that a 
great change in the nature of the growing thing being 
studied could lead to a new growth period. That sucha 
change is now taking place in airplane design is certain. 
For the space of time given by the data of Fig. 4 to 
which the growth curve is fitted, the speed of each and 
every airplane was determined by the subsonic inter- 
section of drag and thrust curves of a similar nature. 
Gradual advances changed the constants determining 
the high speed but did not have the effect of changing 
the nature of the laws governing the performance. De- 
sign improvements were the cause of these gains. They 
include the transition from biplanes to externally braced 
monoplanes to cantilever monoplanes (a reduction in 
structural thickness to span ratio considerably more 
drastic than any proposed today for supersonic flight) ; 
elimination of open cockpits; retractable landing gear; 
cowling and cowl flaps for engines; variable-pitch pro- 
pellers; and, probably more important, a steady in- 
crease in specific ratings and maximum ratings of the 
reciprocating engine. Fig. 6 has been prepared to show 
several aspects of this process of improvement. The 
curves for the subsonic aircraft have similar shapes but 
cannot be clearly compared on the same scales because 
of the minute thrust and resistance of the Wright broth- 
ers’ first airplane. Also shown in Fig. 6 are the types of 
thrust and drag variations to be expected from future 
supersonic airplanes. The great values of thrust and the 
increase in thrust due to ram for jet engines compensate 
for the increase in drag with speed. The subsonic inter- 
section of the thrust and drag curves is eliminated, thus 
opening up an extensive new speed region, heretofore 
beyond reach. This completely new speed régime will be 
the cause of the new rapid development of the speed 
growth curve. It is interesting to note that the first pe- 
riod of growth was preceded by a period in which the 
thrust required for flight had to be obtained from means 
other than the installed power plant. A component of 
the weight obtained by gliding downward or energy 
taken from the atmosphere was used to provide power. 
The supersonic era is also being ushered in by “simple” 
means of obtaining power, including rocket motors and 
dives. Supersonic aircraft of economic and military use- 
fulness will result from the development of suitable air- 
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breathing engines, just as in the case of subsonic air- 
planes. 

An unusual speed prediction was recently made using 
data collected by Prof. Hornell Hart, of Duke Univer- 
sity.'" His data, shown in Fig. 7, were faired by Time 
Artist R. M. Chapin, Jr., to show a speed of 1,700 m.p.h. 
by 1960. The growing things studied by Hart were not 
just airplane speeds but were velocities of several man- 
carrying vehicles. Such an approach is novel in the lit- 
erature studied during the preparation of this paper. 
Hart states that the airplane speed curve appears to 
have several subcurves. The jet airplanes, as pointed 
out above, may well make an additional subcurve. An 
appreciation of the growth of airplane speeds is ob- 
tained by noting that the airplane speed record curve 
does not cross the horse curve until 1909 and the rail- 
road curve until 1915. 


SPEED, CEILING, WING LOADING, AND SAFETY 


It has just been pointed out that a new speed growth 
curve for jet-propelled aircraft is starting. This new 
period of growth will carry speeds well into the super- 
sonic-speed régime. The airplanes of this future era will 
certainly fall into design and performance patterns, 
much as has been the case with subsonic aircraft. 
Rough calculations will be made in an attempt to pre- 
dict the simpler of these patterns. The performance 


characteristics, speed and altitude, will be studied first. 

Before examining the supersonic case, however, the 
reasoning will be applied to subsonic airplanes. These 
will be assumed to have conventional characteristics. 
There is, therefore, a limit to maximum speed set by the 
drag rise occurring in the transonic region, to overcome 
There is 


which sufficient power cannot be obtained. 
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also a lower bound to the minimum speed arising from 
the limiting lifting area which can be fabricated from, 
given weight of structural material. This is a conse. 
quence of the strength of materials. Stated differently 
the lower speed limit is a result of the fact that the lem 
est possible wing loading is of necessity greater than the 
minimum unit wing weight. 

It will be of some interest to calculate, for the syb. 
sonic region, the wing loading for maximum lift to drag 
ratio. Assuming a parabolic drag polar, the maximum 
lift to drag ratio corresponds to a condition where the 
induced and parasite drags are equal.'® 


Cpp = D ‘wA.R., (1) 
in which 
Cpp = parasite drag coefficient 
Crr/p = lift coefficient for maximum lift to drag ra- 
tio 
A.R., = effective aspect ratio 


To obtain the desired rough approximations for the per- 
formance limits, we shall assume that the values of 
Cpp and A.R., are constant for all subsonic airplanes. 
Then, the variables remaining are speed, altitude, and 
wing loading. For horizontal unaccelerated flight the 
lift will equal the weight, and Eq. (1) is equivalent toa 
relation between the speed, atmospheric density, and 
wing loading for maximum L/D. 


= (pV2/2) A.R.g (2) 


While it is not necessary that an airplane operate ac- 
cording to Eq. (2), the penalties for great deviations 
may be large, since the parabolic drag polar has the 
property 


(D/L) | 4 (3) 
(D, 2 Ch Cr D 
Since 
(W/S)/(pV2/2) 


(4 
[(W/S)/(pV?/2) 
it is possible to easily compute the penalty in lift to drag 
ratio which results from deviating from Eq. (2). Ata 
given speed and altitude the variation of minimum drag 
to lift ratio with wing loading is given in Table 1. 


TABLE 


 (D/L)/(D/L) min 
1.0; 1.000 1.000 
0.8; 1.250 1.025 
0.6; 1.667 1.133 
0.4; 2.500 1.480 
0.2; 5.000 2.600 
0.1; 10.000 5.050 
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Fig. 8 has been prepared to show the performance 
limits under discussion. Optimum wing loading lines 
according to Eq. (2) are shown, that for | lb. per sq.ft. 
being used to indicate the approximate minimum speed 
limit. For these curves, A.R.,. = 10 and Cp, = 0.02 
have been used. 

The ceiling is determined by the fact that the maxi- 
mum power-plant output falls off with increasing height 
at a different and greater rate than the minimum power 
required. It is well known that at full throttle the alti- 
tude variation of power output for a reciprocating type 
aircraft engine approximately follows the empirical re- 
lation 


P/P' sr = (o = 0.1) 0.9 (5) 


in which P’s; is the full throttle output that would be 
obtained at sea level. If the engine is restricted to part 
throttle at sea level, with an output Ps,, then 


P 


= o — 0.1 
Psi 


(6) 
will hold above the critical altitude. o is the at- 
mospheric density ratio. Ps, is taken here to be the 
constant power rating that holds below the critical alti- 
tude, at which ¢ = ocr. The weight of the complete 
power-plant installation will be assumed to be 


= (7) 


which will further be assumed to be some fraction \ of 
the total airplane weight WV. Then, 


For a given weight of airplane and for some values of 
the constants K and \, Eq. (8) gives the altitude varia- 
tion of maximum brake horsepower above the critical 
altitude. 

Since the drag-weight ratio for minimum power re- 
quired is 


D/Wmp = (9) 

the power to weight ratio for minimum power is 

P/W = (V/375y) (10) 
and at the ceiling, 

Vavg] 16Cpp _»® — 0-1 

3759035 AR,  K oce — 0.1 an) 
This ceiling limitation has been plotted on Fig. 8, using 
n = 0.85, Cyp = 0.02, A.R.. = 10, A/K = 0.13, and 
ocr = 0.4. The wing loadings for the points on this 
curve are 1/3 times greater than those shown for maxi- 
mum L/D. 

For a jet-propelled subsonic airplane, the ceiling will 
approximately correspond to the condition that the min- 
imum drag to lift ratio will equal the ratio of maximum 
thrust to gross weight. For that condition 


(D/W) min. = (D/L) min. = 2“ Cpp/#A.R.. (12) 


The thrust output of a gas turbine, above about 35,000 
ft., follows the approximate relation 


If c is the gas-turbine installed weight per unit of sea- 
level thrust and ) is the fraction of the gross weight al- 
lotted to the power plant, 


T/W = 1.2 (d/c)o (14) 


Equating to the minimum D/L, 
1.2 (A/c)o = (2/n)W Cop/wA.R.g (15) 


Here » represents the installation efficiency factor and 
takes into account duct losses and, in addition, the vari- 
ation of thrust with speed. The fact that Eq. (15) does 
not depend on speed, as did Eq. (11), is a result of the 
assumption of constant thrust for the gas turbine. 

Eq. (15), defining the ceiling of subsonic jet air- 
planes, is plotted on Fig. 8, using A/c = 0.5, » = 0.90, 
Cpp = 0.02, and A.R., = 10. Both of the relations, 
Eqs. (11) and (15), are unusual because they define ceil- 
ings in terms of the main aerodynamic, power plant, 
and weight parameters in general fashion. 

Simple considerations have now shown that subsonic 
airplanes operate in a limited region of the speed-alti- 
tude diagram. The minimum speed is determined by 
structural weights; the maximum speed, by the drag 
associated with compressibility; and the ceiling, by the 
drop-off in power-plant output with altitude. It must 
be kept in mind that the limits shown are not those for 
any single airplane but rather are those of various air- 
planes whose wing loadings vary with speed and alti- 
tude. Since the regions bounded by these limits are 
reasonably close to the performances actually obtained 
by subsonic airplanes, the same approach will now be 
used in the supersonic régime. Some characteristics 
of the supersonic drag polar will be discussed in addi- 
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tion, since no textbooks have yet appeared giving the 
information needed here. 
Subsonic airplane performance calculations are based 
on the type of drag variation of the elliptic plan-form 
wing. This gives a constant drag coefficient plus one 
varying with the square of the life coefficient. Drags 
other than that of the wing are included in the two wing 
terms. For supersonic performance calculations, the 
two-dimensional wing will serve as the basis. Its drag 
coefficient!’ is 

t/c)? 
Cp = Kpr 1 + WANE —1+ (16) 
The first term is the profile drag associated with the dis- 
tribution of normal pressure at zero lift. Its magnitude 
depends on the airfoil profile shape. The second term 
represents the drag associated with normal pressures at 
an angle of attack, independent of profile shape. The 
third term is the skin friction. The variation of profile 
drag with the plan-form variables can be introduced 
into the first term, since such drags are also dependent 
on the factor (t¢/c) 2// Mt? — 1. By combining this term 
with the skin friction term, a profile drag coefficient is 
obtained similar to the subsonic one. The second term 
can also be corrected for plan form and other effects to 
give a coefficient depending on C,” just as in subsonic 
performance calculations. The induced drag correc- 
tions can be grouped to form an induced drag efficiency 
factor, e. 
Because of the rearward shift of the wing center of 
pressure at transonic speeds, the effect of tail loads on in- 
duced drag can have an important effect on e. Simple 
calculations show this factor to be 


1— 2d/l) + d/l)? 
1 — (de/da) (d/h) + (Sw/ Si) (d/l)? 


d 
(17) 


W/qS 


Here, e has the definition 


Co, = (Cr2/4e VM? — 1; Cy = (18) 


and é,is the portion of e resulting from the trimming tail 
load. d/l, is the distance from the center of gravity to 
the wing-fuselage center of pressure, measured in tail 
lengths with both dimensions positive when the center 
of pressure and the tail are behind the c.g. de/da is the 
wing downwash derivative at the tail. S,/.S, is the ra- 
tio of horizontal tail area to wing area. The importance 
of this effect can be seen by noting that values of d/l, = 
0.1 are not only possible but are probable. The trim- 
ming tail load thus can add 25 per cent to the induced 
drag. 

The effect of the wing plan form on e is obtained for 
sharp leading-edge wings from the slope of the lift 
curve, since then 


(19) 


Co; = C, 
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The plan-form contribution to e is 
= dimensional (20) 


The final value of e is the product of ép, @, and som 
other factors that take into account the variation of 
fuselage drag with attitude, etc., just as in the sub. 
sonic case. 

It is suggested that, as a consequence of the natural 
division of the supersonic drag into terms independent 
of lift and dependent on lift squared, the terminology of 
subsonic performance analysis be retained. The firsj 
term may thus be called the parasite drag; the second 
can be called induced drag. This is believed to be more 
practical than adaptation of new terms having to do 
with wave drag. It is obvious that, with the parabolic 
drag polar, not only the subsonic terminology can be 
retained but also all the formulas having to do with 
maximum lift-drag ratio, minimum power, etc., as dis- 
cussed in reference 16. 

Using suitable values in the relations just developed, 
the ceiling characteristics of supersonic aircraft have 
been calculated and are shown in Fig. 8. It is seen that 
ceilings of over 50,000 ft. can be obtained and the opti- 
mum wing loadings tend to be high. There is also 
shown in this figure the variation of speed and altitude 
for which, on an Air Force hot day, a stagnation tem- 
perature of 300°F. will be reached. There can be no 
doubt of the serious nature of this temperature effect. 
It will force structural weight increases and make nec- 
essary the provision of huge refrigerators for cooling the 
crew, communicating equipment, instruments, engine 
accessories, batteries, tires, and all sorts of other equip- 
ment. It is estimated that a penalty of 10 per cent of 
weight-empty will be required to design for continuous 
operation with a 300°F. stagnation temperature. 
Further increases in speed will rapidly increase this pen- 
alty because of the fact that increase in temperature de- 
pends on the square of the velocity. 

Developments beyond those just discussed are cer- 
tainly going to take place. A crude estimate for these 
has been obtained in the following way. First, the sub- 
sonic airplanes were assumed to have afterburners 
installed, and, second, the power-plant thrust-to-weight 
ratio has been doubled. The doubled thrust was shown 
by Clauser® to be possible over a 7-year period. Fig. 9 
shows the results of these calculations, together with 
peak performances of 1903 and 1949. Rocket-powered 
aircraft, of course, should be able to achieve much 
higher speeds and higher ceilings than those shown on 
Fig. 9. The calculations made here are confined solely 
to airplanes using air-breathing engines, since only 
these vehicles will have usable endurance, at least with 
presently foreseeable rocket fuels. 

The calculations also show some of the other charac- 
teristics of the supersonic airplanes. On Fig. 10, for 
example, are plotted the maximum lift-drag ratios and 
lift coefficients for maximum lift-drag ratios which will 
be properties of these airplanes. While the magnitude 
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of the lift-drag ratio is small compared with those for 
conventional aircraft, it should be remembered that 
range depends on the product of lift-drag ratio and veloc- 
ity, so that the great speed of a supersonic airplane in 
some measure compensates for the reduced aerodynamic 
eficiency. Increases in size will also lead to better aero- 
dynamic efficiency. In addition, the normal rate of en- 
gineering improvement will provide better performance 
than that predicted here. 

The wing loading for best lift-drag ratio for these 
small supersonic airplanes is shown in Fig. 11. It can 
be seen that enormous wing loadings are required to 
achieve maximum lift-drag ratios at sea level. Re- 
quirements for maneuverability, take-off, and landing 
will prohibit the use of such high wing loading. Penal- 
ties as shown in Table 1 will then reduce the per- 
formance at sea level, so that actual lift-drag ratios 
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will no doubt be in the order of magnitude of one-half. 
The possibility of obtaining great endurance at super- 
sonic speeds at low altitudes thus appears to be remote. 

The general trends of the wing loadings shown are in 
the direction of higher loadings than those now in use. 
While it would be desirable to have efficient fast aircraft 
with low wing loadings, it is nevertheless true that the 
trends in aircraft design have indicated that higher per- 
formance, improved safety, and higher wing loading can 
go hand in hand. On Fig. 12, for example, is shown the 
wing loading and speed of fighter airplanes. To aid in 
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the extrapolation, a value for the X-1 airplane has also 
been shown, based on data given by Aviation Week.'* 
A steady increase in wing loading can be seen to have 
taken place and certainly will continue to take place, but 
it will be accompanied by higher speeds. 

Much the same sort of wing-loading trend can be 
shown for transport aircraft. Tremendous improve- 
ments in their safety record have nevertheless been 
achieved. On Fig. 13 are plotted the data that show the 
safety record of the United States domestic scheduled 
air-line operators. I predict that, in spite of the fact 
that airlines are now operating aircraft with unprece- 
dentedly high wing loadings, the safety record will con- 
tinue to show great improvement. There may be, in the 
next decade, several years in which there will be no pas- 
senger fatalities on the domestic scheduled air lines. 
The curve on Fig. 13 has been extended into the future, 
using a constant percentage rate of improvement. A 
tendency to level off may begin to appear but cannot be 
predicted by any means known to the author. 
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CONCLUSIONS 


Progress in the aeronautical sciences has been conti, 
uous and rapid since the flights of the Wright Brother 
in 1903. The engineering technologies, including je 
propulsion and high-speed aerodynamics, are sufficiently 
advanced to allow soon the exploitation of the super. 
sonic region, first with military aircraft and later with 
passenger-carrying vehicles. Top speeds, cruising 
speeds, and ceilings will be increased greatly, and ther 
will also be a continuation in the improvement of flying 
safety. These goals are a challenge to those working in 
the aeronautical sciences, a challenge that will surely be 
met. 
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The Calculation of Supersonic Downwash 
Using Line Vortex Theory 


HAROLD MIRELS* AND RUDOLPH C. HAEFELI* 
Lewis Flight Propulsion Laboratory, N.A.C.A. 


ABSTRACT 


A straightforward procedure, based on lifting-line theory, is 
presented for obtaining the downwash behind lifting surfaces 
at supersonic speeds. The customary assumptions of an inviscid 
fluid and small perturbation velocities are made. 

In particular, the vortex field associated with a supersonic 
lifting surface is examined; the vertical perturbation velocity field 
due to a line vortex segment of constant slope is presented in ex- 
plicit form and is utilized to discuss some of the properties of line 
vortices in a supersonic free stream; and a bent lifting-line 
method and an unnent lifting-line method (horseshoe vortex sys- 
tem) are proposed for determining the downwash field behind 
swept and unswept wings, respectively. The lifting-line methods 
are applied to triangular and rectangular wings and give results 
that are in close agreement with the exact linearized solutions. 

This paper is an abbreviated version of reference 1. 


INTRODUCTION 


es METHODS, based on linearized theory, have 
been presented for obtaining the downwash_be- 
hind supersonic wings. These employ conical super- 
position,” doublet distributions,** and vortex distribu- 
tions. Each method provides integral expressions for 
downwash which are generally tedious to evaluate. 
The complexity of these expressions indicates the need 
for a straightforward procedure for obtaining reason- 
ably accurate, if not exact, downwash solutions. 

A logical approach to the development of a simplified 
supersonic downwash theory is to derive the supersonic 
analogues of the line vortex procedures that have proved 
valuable in incompressible flow theory. Certain dif- 
ferences exist, however, between the properties of vor- 
tices in an incompressible flow field and vortices in a 
supersonic flow field. These differences must be in- 
vestigated before an extension of subsonic techniques 
is possible. The present paper has two main objec- 
tives: (1) The downwash field induced by a line vortex 
segment of constant slope, obtained from an integral 
expression provided by Robinson,°® is presented in ex- 
plicit form and is utilized to discuss some of the prop- 
erties of line vortices in a supersonic free stream; (2) 
lifting-line methods for computing supersonic downwash 


Presented at the Aerodynamics Session, Annual Summer 
Meeting, I.A.S., Los Angeles, July 21-22, 1949. 

* Wind Tunnel and Flight Division. 

t An alternate approach is to derive the velocity field due to the 
line vortex by using an equivalent surface distribution of dou- 
blets. Schlichting? obtained the flow field about a horseshoe vor- 
tex in a supersonic stream in this manner. 


are presented, and several calculations, based on these 
methods, are compared with the exact linearized solu- 
tions. 

Zero thickness wings (lifting surfaces) are considered 
throughout. The wings are assumed to induce small 
perturbation velocities in an inviscid fluid so that the 
linearized equations of motion are applicable. 


SYMBOLS 


Only important symbols are defined here. 


HHS fT = Cartesian coordinate system (Fig. 1) 
Yor 20) 
u,v, W = perturbation velocities 
¢ = components of vortex vector (three-dimensional 
field) 
9’, = components of vortex vector (vortex sheet) 
b = wing span 
Cp = pressure coefficient 
— VY; — X,) 
— m Xi)? + (1 — 
M = free-stream Mach Number 
m = slope of line vortex or lifting line (dy/dx) 
ri = VX? — BY? — 
Xi = x— X; 
Y; 
= 3; 
a = angle of attack 
3 = V/M?-1 
€ = integration interval 
K = circulation 
¢ = perturbation velocity potential 
Ag = or — 
p = free-stream density 
r = bound circulation 
Dus = bound circulation at midspan 
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Subscripts 

a,b = points of intersection of forward Mach cone 

with lifting line (or edge of vortex sheet) 

B = bottom surface of s = 0 plane 

l = 1,2,3... points on line vortex 

l = plan-form leading edge 

T = top surface of zs = 0 plane 

t = plan-form trailing edge 

0 = variable of integration 
Superscripts 


slope at bend in line vortex when approaching 
from negative y-direction 

slope at bend in line vortex when approaching 
from positive y-direction 


m = 


VORTEX FIELD RELATIONS 


The definitions of vorticity and circulation apply 
equally to subsonic and supersonic flow fields. Several 
defining relations for a volume distribution of vorticity, 
a surface distribution of vorticity, and a line vortex are 
presented in the following section. 

The velocity field is assumed to consist of a supersonic 
free-stream velocity U plus small perturbation velocities 
u,v,andw. Invector form 


g=(U+ u)ti+ vj + wk (1) 


The coordinate system is illustrated in Fig. 1. The 
vortex vector field @ is defined as the curl of the velocity 


vector. Thus, 
®=VXg=t +ajt+ (2) 
where 
ow dv Ou Ow Ou 
a * ~ ox oy 


The magnitude of @ is referred to as vorticity. Vortex 
lines are lines that are tangent at all points to the local 
vortex vector and are determined from 


dx/§ = dy/n = dz/f (3) 


Vortex lines cannot terminate within a flow field but 
must either form a closed curve or extend to infinity 
or a boundary of the field. Circulation is the line in- 
tegral about a closed curve 
K= udx+ody+wds (4) 
If vorticity exists as a surface of velocity discon- 
tinuity in the z = 0 plane, then outside this plane &, », 
and ¢ are all zero. In the zs = 0 plane, ¢ = 0, whereas 
£ and 7 are infinite. However, the limits 


= lim édz, n’ = lim ndz 
dz—>0 dz—>0 
are finite and equal? 
t’ = vp — 07, n = Up — Up (5) 


The subscripts 7 and B designate velocities on the top 
and bottom surface of the z = 0 plane, respectively. 
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The vortex lines through all points on an infinitely 
small closed curve bound a vortex “‘tube.’”’ The cir. 
culation about any path encircling the tube is constant 
and equals K = wo, where w is the vorticity and ois the 
normal cross-sectional area at any point along the 
tube. A line vortex is generated by allowing the cross 
section of the vortex tube to approach zero while maip. 
taining wo constant. Thus, the distinction between 
a line vortex and a vortex line is that the line vortex 
has an assigned constant circulation A, whereas the 
vortex line indicates the direction of the local vortex 
vector at each point along its length. 


LINEARIZED SUPERSONIC WING THEORY RELATIONS 


The vortex field generated by a lifting surface in a 
supersonic stream will now be considered. 


Velocity Potential 


If the boundary conditions for a lifting surface are 
specified in the z = 0 plane, the u and v velocities are 
antisymmetric and the w velocities are symmetric 
about this plane.® Thus, the velocities on the top and 
bottom surface of the z = 0 plane (considering each 
surface separately) are related by = V7 = 
andw, = wg. The discontinuities in u and v constitute 
a vortex sheet. If it is assumed that the field, excluding 
the vortex sheet, is everywhere irrotational, a perturba- 
tion-velocity potential ¢ can be so defined that 


lp = —dx + —dy+ —d 


(6) 
=udx+vdy+wdsz 


When the undisturbed flow field upstream of the wing 
is considered to be of zero potential and the boundary 
condition requiring uy = ug = O off the wing (uz = 
—uz ~ 0 implies lift) is applied, the potential in the 
z = O plane may be obtained by integrating, along 
lines of constant y, 


x x 
¢r = f ur dx, ¢s = fi Uy dx (7) 
xl xl 


where x, is the equation of the leading edge as a function 
of y. From Eq. (7) it may be concluded that: 

(1) or = ¢p = O everywhere in the z = 0 plane ex- 
cept behind the wing leading edge. 


LINES OF CONSTANT POTENTIAL FOR RECTANGULAR WING 
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(2) Ata given span station, g7 and ¢, remain con- 
stant for all values of x downstream of the trailing edge. 
Lines of constant potential for a rectangular wing are 


shown in Fig. 2. 


Vortex Lines 


The equation for the vortex lines, from Eq. (3), is 


n' dx — t’dy = 0 


Substituting, 
—(vr Vp) _9 Wr On 
oy oy 
9 = up = = 2 
Ox Ox 


dx + Ovr “n= vn dx + 

Ox oy Ox oy 
Comparison with Eq. (6) shows that Eq. (8) represents 
lines of constant potential. Thus, vortex lines coincide 
with the lines of constant potential in the z = 0 plane 
(Fig. 2). 


Circulation 


The circulation included between two points x),; and 
%,J2 On a Wing can be obtained by the integration indi- 
cated in Eq. (4). The path of integration is arbitrary, 
except that the path should cross the z = 0 plane only 
at the two specified points. Taking the integral along 
the top and bottom surface of the z = 0 plane, 


x 


X2,y2 
k= f (up dx + vp dy) + (up dx + v, dy) 


1M X2,V2 
(¢7,2 ¢r,1) + (¢p,1 
= Ag — Ay, (9) 


where Ay equals gy — ¢, and represents the jump in 
potential at the point. The quantity Ag is, in fact, 
the doublet strength® so that the net circulation between 
two points in the s = O plane equals the difference be- 
tween the doublet strengths at these points [Eq. (9)]. 
The equivalence of a doublet distribution and a vortex 
distribution indicates that the flow about a lifting sur- 
face can be calculated on either basis. Schlichting’ 
computed the velocity field due to a supersonic horse- 
shoe vortex by using an equivalent surface distribution 
of doublets. 


Circulation and Lift 


The wing lift per unit span is given by the chordwise 
integration 


L(y) = : pur (Cy, — Cp, r)dx (10) 


LINE VORTEX SEGMENT EXTENDING FROM (x,,y,) TO (xo ,y¥o) 


SLOPE 


Ko - x \ 
Yo 
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(x,y,z) 


FIG. 3 


Substituting the linearized expressions for pressure co- 
efficient 


U dx U ox 
and integrating, Eq. (10) becomes 
L(y) = pU ant 
= pUT (11) 


where Ag, is the doublet strength at the trailing edge 
and I is the bound circulation. Also, 


dV /dy = d(Ag,)/dy = —(&’), (12) 


Eq. (12) relates the shed vorticity to the bound circula- 
tion. Eqs. (11) and (12) are familiar incompressible 
flow relations. 


UpwasH DUE TO SUPERSONIC LINE VORTEX 


Line Vortex Segment 


Complicated velocity fields can be generated by the 
linear superposition of relatively simple fields. It will 
therefore prove useful to examine the upwash field due 
to a line vortex segment of constant slope m which ex- 
tends from x,y; to X2,y2 in the z = 0 plane. 

The upwash at x,y,z, a point whose forward Mach 
cone contains the entire line vortex segment (Fig. 3), is 


w= (K/2r)[G2(m) — Gi(m)]} (13) 


where 


G,(m) = - mX) (BMV Xd) 
— mX,)? + (1 — B'm?)Z;"] 


The limits are chosen so that ye > y,. The strength K 
is positive when the vortex vector (applying right-hand 
rule) has a component in the positive y-direction. Fora 
line vortex parallel to the x-axis, the strength is posi- 
tive when the vortex vector is in the positive x-direction. 


a 
7 
Xo 
(x2, 
j 
j 
(7) 
ne ex- 
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INFINITE LINE VORTEX 


(x, yi) 


(x, y,2) 


(b) SUBSONICALLY 
INCLINED |Bm|<! 
FIG. 4 


(a) SUPERSONICALLY 
INCLINED 


Eq. (13) applies for any line vortex segment of constant 
slope. If the line vortex segment terminates on, or is out- 
side of, the forward Mach cone from x,y,z, the corres pond- 
ing limit [G,(m) term] in Eq. (13) is neglected.' 


Infinite Line Vortex 


Supersonically Inclined (|8m| > 1).*—An infinite line 
vortex inclined supersonically to the free stream is indi- 
cated in Fig.4a. Both limits in Eq. (13) are outside the 
forward Mach cone and therefore neglected. Thus, the 
upwash and, in fact, all the perturbation velocities in- 
duced by this line vortex are zero. This result agrees 
with oblique airfoil theory, since the perturbation ve- 
locities are zero downstream of the two Mach waves 
from the trailing edge of a two-dimensional airfoil in- 
clined supersonically to the free stream. (The line 
vortex may be considered as the limit of a two-dimen- 
sional airfoil whose chord goes to zero at constant lift.) 

It may appear that the circulation integral [Eq. (4) ] 
taken about the supersonically inclined line vortex will 
yield zero and thus violate the statement that the 
circulation is equal to the strength of the line vortex. 
However, the proper value of circulation will be obtained 
if the necessary precautions are taken when the path of 
integration crosses the Mach waves from the line vortex 
(that is, take the line integral about a lifting surface 
whose chord goes to zero at constant lift). A value of 
K/2 is obtained when the path crosses a Mach wave. 

Subsonically Inclined (| < 1).—An infinite line 
vortex inclined subsonically has one limit that extends 
to infinity but remains within the forward Mach cone 


(Fig.4b). From Eq. (13), the upwash for 0 < Bm < 1is 
w= —(K/2r) G,(m) 


_K V1 = B'm(¥, — mX,)_ 


- — (14) 
(Y; — mX,)? + (1 — 


where x;,y; is a point on the line vortex. The upwash 
is infinite along the line vortex, and the flow circulates 
about the vortex in a manner similar to that for a vor- 


* A line vortex is supersonically or subsonically inclined as the 
component of the free-stream velocity normal to the line is super- 
sonic or subsonic, respectively. 
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tex in an incompressible field. For m = 0, Eq, (14) 
becomes 


w = + 2)] 


which is identical with the expression for upwash due ty 
an infinite line vortex parallel to the stream (along y = 
y;) in an incompressible field. 

These results indicate that the behavior of the jp. 
finite line vortex is completely different from that of the 
incompressible flow vortex when Bm| > 1. However, 
when | Bm! < 1, both vortices have similar flow fields 
and are, in fact, identical for m = C. 


Bent-Line Vortex 


A bent-line vortex can be considered as the super. 
position of two line vortex segments (Fig. 5). The up- 
wash, utilizing Eq. (13), is 


w = (K/2r)[Gi(m-) — Gi(mt)] (15) 


where m,~ and m,* designate the slopes of the line vor- 
tex before and after the bend at x,y,. The term 7 ap. 
pears as a factor in the denominator of Eq. (15), so 
that the upwash exists only in the downstream Mach 
cone from x,y; and is infinite on the cone surface (ex- 
cept in the z = 0 plane). 

If m~ = o and m+ = 0, a horseshoe vortex results 
(Fig. 6a). The upwash is zero, except within the Mach 
cone from the bend. The stream-line pattern in a 
plane normal to the free stream (section AA) is shown 
in Fig. 6b.7. At the cone surface, the stream lines are 
radial to the center. The perturbation velocities are 
infinite at the cone surface except in the z = 0 plane. 
The spanwise distribution of upwash in this plane 
is shown in Fig. 6c for s = O and z = 0.3R (R is radius 
of Mach cone at section AA). 


LiFTING-LINE THEORY 


The vertical perturbation velocity field due to a bent- 
line vortex is given by Eq. (15). The vortex distribu- 
tion associated with a lifting surface has been discussed. 
These relations will be used to develop lifting-line 
methods for computing supersonic downwash assuming 
the load distribution for the lifting surface is known. 


BENT LINE VORTEX 


FIG.5 
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PERTURBATION VELOCITIES DUE TO HORSESHOE 
VORTEX 
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(a) HORSESHOE VORTEX 


(b) STREAMLINE PATTERN (SECTION AA) ” 


w 
| 


| z=.3R 


(c) SPANWISE DISTRIBUTION OF UPWASH (SECTION vA) 
FIG.6 


The load distribution for a large variety of plan forms 
can be found by the methods of Evvard’ and Heaslet 
and Lomax.® 

A lifting line concentrates the chordwise loading into 
aline. Thus, the bound circulation is represented by a 
line vortex of variable strength ([ = Ag,), whereas 
the trailing vortex sheet maintains the same vorticity 
but now originates at the line rather than the wing trail- 
ing edge. 


Bent Lifting Line 


A bent lifting line approximating the section centers 
of pressure (Fig. 7a) seems to be a reasonable representa- 


LINE VORTEX REPRESENTATION OF BENT LIFTING LINE 


SECTION CENTERS 
OF PRESSURE 


(b) LINE VORTEX REPRESENTATION 
FIG. 7 
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tion of a swept wing. The vortex distribution asso- 
ciated with the bent lifting line can be represented by 
the superposition of bent-line vortices as indicated in 
Fig. 7b.* The vortex sheet is considered to consist of 
elemental bent-line vortices each of strength dK = 
+(dI'/dy,)dy,. The contribution, to the downwash 
field, of each elemental bent line vortex is [from Eq. 


(15) | 


dw = [G,(m) G,(0) | dy, 


2 dy, 


where the subscript o designates a variable of integra- 
tion and refers to the coordinates at the bend. The 
vertical perturbation velocity field due to the bent 
lifting line of Fig. 7 is the sum of the contributions of the 
elemental vortices and the bent-line vortex of strength 
I, (bound circulation at midspan) and equals 


| 
w= = [G,(m~) G,(m*) ~ 
“7 


: [G.(m) — dy, (16) 


IMPROPER INTERVALS OF DOWNWASH INTEGRAL 


(xo,¥a) 


(a) INTERSECTION OF FORWARD 
MACH CONE WITH LIFTING (b) INTERSECTION OF FORWARD 
LINE MACH CONE WITH EDGE OF 


Yq VORTEX SHEET 


(x,y, 0) 


(c) POINT ON VORTEX SHEET 
FIG. 8 


The integration is conducted along the lifting line in- 
cluded in the forward Mach cone. The nonintegral 
term of Eq. (16) is readily evaluated. The integral 
term may be evaluated by analytical or mechanical 
means. When mechanical methods are used, the singu- 
larities in the integrand must be isolated. Suitable 
procedures for isolating the commonly encountered 
singularities are as follows: 

(a) Singularity due to intersection of forward 
Mach cone with lifting line. The integrand of Eq. 
(16) is infinite at the intersection of the forward Mach 
cone and the lifting line (for zs 0). The contribution 
éw of the integral, for the interval y, — & < Yo SM 
(Fig. 8a), to the vertical perturbation field may be 


written 


* Suggested by H. S. Ribner. 
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Mb 
— [G.(m) — G,(0) /dy,)dy, (17a) 
24 J dy, dy, Ve 
dv \ dy, 
—(G, — G,(0) |r 17 


where ¢, is a convenient length. The first integral in Eq. (17b) is proper and may be mechanically evaluated (the 
integrand is zero at y, = y,). The second integral is 


dy, 2m (1 — B’m”*) 
, O< < 1) 
m +/26,/(mx y), (Bm = 1) 


dye sin~ Us (bm > 1) 
To V ‘NG 2[m(x — B’my) — (1 — B?m?)y| 
where m is the slope of the lifting line at y,._ A similar procedure applies at y,. 
(b) Singularity in dI'/dy, at wing tip. The vertical perturbation velocity at a point whose forward Mach cone 
intersects the edges of the trailing vortex sheet (Fig. 8b) is obtained by integrating between the limits y, = +b/2. 
If dI'/dy, is singular at these limits (a half order singularity is common), the singularity may be isolated by a pro- 


cedure similar to that used in Eq. (17b). Thus, the integral in Eq. (16) may be written, for the interval 
[(b/2) <y¥, < 


Il 


b/2 av 
iw = — — ({G,(m) — G,(0)] — — G,(0)]}y,=6/2) —- dy. + 
T Sb dy, 


/2—e 


1 
{[G.(m) — (T) (18) 


The integrand in Eq. (18) equals zero at y, = 6/2. 

(c) Singularity at points on vortex sneet. The Cauchy principal value of the integral [Eq. (16) ] is required 
for determining w at a point on the trailing vortex sheet. Consider the interval y — « < y, < y + «€ (Fig. 8c). 
The integral, in expanded form, is 


y+ 2dr 
= = VX, — = dy, (19) 
Qe — mX,) Wo 
If dI'/dy, can be approximated by the first two terms of the Taylor’s expansion 


a Vo 
and eis sufficiently small that 
meV X,? — — moX,) 
Eq. (19) may be approximated by 


Jy—e dy, yo=y dy,” yo=y Y, dy,” 


The value ¢ = 0.05 is satisfactory for most computations. 


Unbent Lifting Line (Horseshoe Vortex System) 


The unbent lifting line (Fig. 9) appears to be a reasonable representation for an unswept wing. The vertical 
perturbation velocities are obtained from Eq. (16) by setting m = © and equal, in the expanded form, 


Yb 
ye ro(x? — + 2°) dy, 


The singularities are isolated as follows: 
(a) Singularity due to intersection of forward Mach cone with lifting line. 
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w= -— + — — + — — sin =. 
L(x? — (Y,2 + + 2°) dy, Br \xY, dyo/y, =y, 2(% — 


(b) Singularity in dI'/dy, at the tip. 


xY,(r.% — [ — 62") | 
Qn (r,(x? — B?2?) (Y,? + 27) — B2”)(Y,? + 27) Jy, 20/2) 


r(x* 6°?) (Y, 24 2°) = yo = b/2 — 


2r 


(c) Singularity at points on vortex sheet. 
dw = = 


The chordwise location of the unbent lifting line yielding the best average agreement with the exact linearized 
solution is still to be determined. According to the techniques used in incompressible flow theory, an unbent lift- 
ing line at the wing center of pressure should be a good approximation. It may be possible to determine the best 
location for each general class of plan forms by comparing the lifting-line solution with the exact linearized solu- 


tion for several representative plan forms. 


Downwash at Infinity 
The velocity field an infinite distance downstream of a lifting line is obtained by allowing x to become infinite 


in Eq. (16). The result is 
Ww Y,? + dy, Yo ( 


Eq. (22) is identical to that for the velocities an infinite distance behind a wing in an incompressible flow field. This 
result has been pointed out by others.*:* The evaluation of Eq. (22) is relatively simple and may be used to ap- 
proximate the downwash several chords behind a supersonic wing. 


ized solution for each wing. In order to simplify the 
—_ expressions, M = /2 (that is, 8 = 1) will be assumed. 
A bent lifting line and an unbent lifting line will be Triangular Wings.—The spanwise distribution of 


used to obtain solutions for the downwash in the y = 0 circulation for a triangular wing of span b and root 
plane behind triangular wings having subsonic leading chord c,, having subsonic [(b/2c,) < 1] leading edges, 
edges. The chordwise distribution of wing loading sug- js [from Eq. (24), Heaslet and Lomax*] 
gests the use of the bent lifting line. The unbent lifting = : 
line will be used for purposes of comparison. An un- P= Ay, = (aUb/E)V1 — (4/b*)y? 
bent lifting line will be used to compute the downwash where E is the complete elliptic integral of the second 
bent lifting line will be placed at the chordwise location 
giving the best average agreement with the exact linear- 

lifting line approximation of Fig. 10a, is 


4aU 


[G.(m* @ =0 + x 


” [G,(m+) — G,(0 
f IG 2 Vo dy, (23) 
0 


LINE VORTEX REPRESENTATION OF UNBENT LIFTING LINE 


9 9 

V1 — (4/b*)y,? 
x) | | | a where (m+) = b/c,. Eq. (23) is obtained from Eq. 
( SUAAAL (16) utilizing the properties of an even function. The 
vertical perturbation velocity field behind an unbent 


lifting line having the same loading (Fig. 10b), is [from 


A | Eq. (21)] 
4aU 
w= —— X 


+ — g2 2 d 
(b) LINE VORTEX REPRESENTATION f % — 3°) Yor (24) 
0 r(x? — + 27) v1 - (4/ 
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(a) BENT LIFTING LINE 
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LIFTING LINE REPRESENTATIONS OF TRIANGULAR WING 
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(b) UNBENT LIFTING LINE 


FIG. 10 
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Rectangular Wings.—The spanwise distribution of 
circulation for a rectangular wing of aspect ratio 
b/c, > 2 is [from Eq. (20), Evvard®], for [(b/2) - 


Cr] Mo b, 2, 
+ 


tal’ b\/b 
A¢, (>. + Cr — 


| ‘ 
c, tan \ (6/2) — Yo / 
Vo +c, — (b/2) 


and, for 0 Z y, Z [(b/2) — ¢,], 
r = 


The vertical perturbation velocities in the y = 0 plane 
behind an unbent lifting line having this loading is given 
by 


° 
= 
(b/2) — cr — 37)(y,? + 27) 


ly, + (b 


(25) 

Com putations.—Eq. (23) was evaluated to determine 
the downwash along the line y = 0, z = 0 and along the 
line y = 0, s = 0.10 for triangular wings of aspect ratios 
2b/c, = 1.6 and 3.2. Eq. (24) was evaluated along the 
line y = 0, s = O for the same wings. Eg. (25) was 
evaluated at points along the line y = 0, z = 0 for ree- 
tangular wings of aspect ratios b/c, = 2 and 4. The 
integrations were performed mechanically for Eqs. 
(23) and (25) and analytically for Eq. (24). The re- 
sults are compared with the exact linearized solutions’~* 
in Figs. 11 to 13. 


DISCUSSION 


The bent lifting-line solutions for the triangular wings 
are shown in Figs. lla and 11b. The agreement with 
the exact solution is very good for all points except 
those within a one-half chord of the trailing edge. The 
agreement indicates that the bent lifting line is a good 
average representation of the triangular wing. The 
method should give even better agreement when used to 
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represent a swept wing, such as the one in Fig. 7a, be- 
cause the bent lifting line would then more closely ap- 
proximate the actual chordwise distribution of 
vorticity. 

The discontinuity in the curves for z = 0.1 b desig- 
nates large negative (upwash) values that become in- 
finite on the Mach cone from the tip of the lifting line. 
Lifting-line solutions always yield this infinity when the 
loading at the tip has a one-half order singularity. 
The infinity may or may not exist in the exact linear- 
ized solution for the flow field about the original wing, 
depending on the plan-form configuration at the wing 
tip. Infinite upwash will exist on the downstream 
Mach cone from a wing tip formed by the intersection 
of a subsonic leading edge and a supersonic trailing 
edge, provided the subsonic edge is not in the free- 
stream direction at the intersection.! 

The unbent lifting-line solutions for the triangular 
wings are presented in Fig. 12. The lifting line was 
placed at the three-quarter chord point to give the best 
average agreement with the exact linearized solution. 
The center of pressure for these wings is at the two- 
thirds chord point. The agreement is not so good as 
that obtained with the bent lifting line and indicates 
that the bent lifting line is more suitable for computing 
downwash behind triangular wings. The smaller as- 
pect ratio wing (2b/c, = 1.6) is in better agreement with 
the exact linearized solution because of the rapidity with 
which the downwash approaches the asymptotic value 
at infinity. 

The unbent lifting-line solutions for the rectangular 
wings are presented in Fig. 13. The downwash ob- 
tained with the lifting line at the one-half chord point 
isin excellent agreement with the exact linearized solu- 
tions. The location of the center of pressure is at the 
4/9 and the 10/21 chord points for aspect ratios 2 and 
4, respectively. These calculations suggest that the 
best chordwise position for the unbent lifting line is at, 
or slightly downstream of, the wing center of pressure. 

A more accurate estimate of the downwash in the 
z = () plane, close to the trailing edge, may be obtained 
by judiciously fairing the curve obtained by the lifting- 
line method to the known value of downwash at the 
trailing edge. At a subsonic trailing edge satisfying 
the Kutta condition, —w/aU is unity, whereas at a 
supersonic trailing edge, —w/al’ can be computed 
using the method of Lagerstrom and Graham.” 


DOWNWASH 


CONCLUDING REMARKS 


The downwash computations based on lifting-line 
theory are in close agreement with those based on exact- 
linearized theory. Thus, lifting-line methods can be 
used to obtain the linearized solution for the downwash 
behind supersonic wings. It should be noted, however, 
that linearized theory assumes the trailing vortex sheet 
to be in the zs = 0 plane and thereby neglects the dis- 
placement and distortion of the vortex sheet. In addi- 
tion, the load distributions indicated by linearized 
theory may be modified by viscous effects and by larger 
magnitudes of perturbation velocities than are per- 
mitted by linear theory. An experimental program, 
such as that conducted by Silverstein, Katzoff, and 
Bullivant'® for wings in subsonic flight, is ultimately 
required to determine the necessary modifications to 
theory which will result in good agreement with prac- 


tice. 
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The Pulse Method for the Determination 


of 


Aircraft Dynamic Performance 


R. C. SEAMANS, JR.,* B. P. BLASINGAME,?7 AND G. C. CLEMENTSON?t 
Massachusetts Institute of Technology and U.S.A.F. Institute of Technology 


SUMMARY 


In this paper a method is described for evaluating an aircraft 
dynamic test involving the response of the aircraft to an applied 
displacement of a control surface. 
can be expressed mathematically by means of a concept called a 
performance dperator' and that with selected types of control sur- 
face motion of common occurrence the response can be related to 
the applied displacement more simply by meaiis of a performance 
function. These two concepts are used extensively by Dr. C. S. 
Draper and his associates in the Instrumentation Laboratory of 
the Massachusetts Institute of Technology. 

A general description is given of the performance operator and 
the performance function, following which a technique is de- 
scribed which permits the performance function of an aircraft to 
be determined from knowledge of its response to a pulse of less 
than 1-sec. duration. The performance function of a U.S. Air 
Force B-25 medium bomber, obtained by this technique, is cor- 
related with the performance function computed from stability 
derivatives. 

It is seen that the conventional analysis gives satisfactory re- 
sults at low Mach Numbers and in the absence of aeroelastic ef- 
fects. However, it is anticipated that, at high Mach Numbers and 
when aeroelastic effects cannot be ignored, experimental tech- 
niques such as outlined in this paper will be required to obtain the 
aircraft performance functions. 


It is shown that the response 


INTRODUCTION 


_ INCREASED AERODYNAMIC DRAG resulting from 
transonic and supersonic speeds can severely limit 
the range capabilities of an aircraft. For this reason 
the stabilizing surfaces of transonic and supersonic 
aircraft are reduced to the minimum size that consider- 
ations of stability wil! permit. Therefore, it is almost 
a matter of course that an appraisal of the dynamic 
characteristics of a high-speed airplane is needed when 
the airplane is still in the design stage. Later, it must 
be expected that the prototype will be given not only 
the conventional flight tests but also dynamic tests 
that will enable artificial stability derivatives to be 
harmonized with the air frame and from which auto- 
matic features for landing, navigation, and_ special 
military purposes can be successfully incorporated. 
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The performance of an aircraft is difficult to analyze 
by strictly mathematical procedures because its re. 
sponse depends upon compressible fluid mechanics and 
involves, even for a rigid air frame, several degrees oj 
freedom. For these reasons a method has been sought 
which will enable the complete performance character- 
istics of an aircraft to be determined from relatively 
simple practical tests. With this in view, the Instru- 
mentation Laboratory, working under a contract with 
the U.S. Air Force, made extensive longitudinal and 
lateral tests on a B-26 (designated A-26 prior to Octo- 
ber, 1948) medium bomber. Assistance in the tests was 
given by the Cornell Aeronautical Laboratory, who 
installed the necessary instrumentation and control 
equipment, conducted the required flights, and reduced 
the resulting flight data. 

In the tests, the response of the airplane was deter- 
mined as the result of separately applied displacements 
of the elevator, rudder, and ailerons. In each case, 
both sinusoidal and step displacements were used in 
order to investigate both steady-state and transient 
effects. The response to a sinusoidal input, as deter- 
mined for several air speeds and center-of-gravity loca- 
tions, agrees reasonably well with the response pre- 
dicted by theoretical methods in both the longitudinal 
mode and the lateral mode at the low Mach Numbers 
used. 

The dynamic performance resulting from step dis- 
placements agreed well with theory in the longitudinal 
tests; the agreement was reasonably good in the lateral 
tests for which the rudder motion provided the excita- 
tion. However, considerable difficulty was experienced 
with the aileron step displacements because, with only 
one degree of aileron, the bank angle reaches 30° before 
the airplane reaches a constant rate of roll. It was sug- 
gested by Major Clementson® that the difficulty could 
be overcome by using a pulse displacement rather than 
a step displacement. With the pulse method, the 
control surface is quickly returned to its equilibrium 
position after initiation of the disturbance. The 
method is simple to use, adequate data for transient 

evaluation can be obtained from only a few flights, and 
the further advantage is offered that the airplane is 
operated at nearly equilibrium conditions. 

In June, 1948, a B-25 medium bomber was made 
available to the Instrumentation Laboratory by the 
U.S. Air Force to make possible an investigation of the 
pulse method. For the purpose, the Laboratory has 
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installed a stop mechanism in the airplane so that the 
pilot can apply elevator pulses. There are two stops 
in the mechanism: One of them is located so that when 
the control column is held against it the pilot can trim 
the elevator for the desired air speed and altitude con- 
figuration; the second stop is adjusted before flight to 
give the desired elevator travel. The pilot introduces 
the pulse by moving the control column from the stop 
giving trimmed flight to the other stop and back again 
to the initial stop in a selected interval of time. The 
stop mechanism is provided with a quick disengage as a 
safety measure and to permit recovery from long-period 
phugoid oscillations that may be incurred at low air 
speeds. 

Instruments have been installed for recording elevator 
position, rate of pitch, and normal acceleration. Dur- 
ing a flight of 2-hour duration, there is time to apply 
several pulses at each of three air speeds and to make 
static calibrations of the recording equipment both 
before and after obtaining the pulse data. The accu- 
racy of the indications of the elevator pulse displacement 
and response is greater than with indications involving 
a step input, because all of the required quantities can 
be considerably greater in magnitude and at the same 
time the total change in the pitch angle can be consid- 
erably reduced. 

The procedure for reducing the pulse data into a form 
similar to that obtained directly from sinusoidal tests 
involves a Fourier Transformation. Computing ma- 
chines, such as the Corodi Analyzer, are designed for 
this specific mathematical application, and other com- 
puting machines may be adapted to this purpose. An al- 
ternate, approximate method is outlined in this paper. 
Using this approximate method, a pulse response may 
easily be transformed by graphical vector summation. 
A computing machine based on this approximate 
Fourier Transform has been constructed. The ma- 
chine permits the real and imaginary parts of this vector 
summation to be recorded continuously as functions of 
frequency, the new independent variable. 


(I) DEFINITION OF THE PERFORMANCE OPERATOR 
AND THE PERFORMANCE FUNCTION 


The airplane in flight can be considered a single 
operating component, which receives certain inputs and 
which produces outputs that are functions of the 
inputs. For present purposes, the principal inputs are 
the control surface displacements and the power-plant 
thrust; the principal outputs are the angle of attack, 
aerodynamic yaw, air speed, angular velocity, and 
linear acceleration. In general, the output quantities 
differ from the input quantities in dimension, magni- 
tude, and integral or differential order. The relation- 
ships of the outputs to the inputs can be determined 
from the kinematic and dynamic equations of motion 
of the airplane. All of these equations can be sum- 
marized in the general form: 


AIRCRAFT DYNAMIC PERFORMANCE 


23 


Qiout) = GinQout) (1) 


The term $P} in this equation designates the perform- 
ance operation required to convert the input to the 
output quantity. The subscript, A, indicates the 
component or system under investigation, in this case 
the airplane, A. The subscripts gj, and Qoue within 
the parentheses show by their designations and order 
the desired input and output quantities. 

As an example, consider the performance of an air- 
plane required when investigating stabilization about 
its roll axis. The input to the airplane in this case is 
the aileron displacement 6,, and the resulting output is 
the angle of roll, g. The airplane performance operator 
that describes the relationship is written [5,9]. 

It is possible to describe the lateral response to 
aileron displacement in terms of three linear differential 
equations with constant coefficients, provided the 
airplane is restricted in motion to small deviations from 
equilibrium. But this restriction does not hold when 
the airplane is in a diving turn. In this case at least 
eight nonlinear equations are required to specify the 
motion of the airplane. 

With the performance operator it becomes possible 
at least to express the angle of roll in terms of the aileron 
motion. In the form of Eq. (1), we can write: 


¢ = $PA(say) 5a (2) 


It must be emphasized that the performance operator 
notation primarily summarizes the operations required 
to establish the output quantity for a given input 
variation. It yields a solution of the performance 
equations only when the input quantity is specified, and 
then only by using conventional mathematical pro- 
cedure. 

The performance of a large number of useful com- 
ponents can be expressed by a single linear differential 
equation having constant coefficients. This equation 
may be written: 


dqiou ) 
dqiin) 


This is one of the possible kinds of relationships between 
the output and the input implied in the performance 
operator definition of Eq. (1)—1.e., 


by d ad? 


ve 21 (out) a2 (out) (4) 
ay dt a J 
The dynamic characteristics are described mathemati- 
cally by Eq. (4) and specify a particular component 
when the coefficients (a’s and 6’s) are known. The 
output response depends upon these dynamic char- 
acteristics and upon the time variation of the input. 
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Several types of inputs are considered in this paper. 
For the special class of sinusoidal input functions, cer- 
tain simplifications in the relation between output and 
input quantities result. It is convenient to use the 
complex exponential notation, e” for these types of 
functions. This term has a unit magnitude and im- 
plies rotation in the complex plane at an angular 
frequency, w, expressed in radians per unit time. The 
angular position of the exponential term, measured 
with respect to the real axis, is equal to the product of 
the angular frequency and the elapsed time, /. It can 
be seen that the complex exponential is periodic, having 
a repetition rate of w/27 and a repetition time, called 
the period, of 27/w. 

An input quantity that has a complex-plane rotation, 
similar to e* but with an initial phase angle (PA) (jn) 
and an amplitude gijn),, can be expressed by 


= eJlwt + (PA) (in)] (5) 


The complex form of Eq. (5) is particularly convenient 
when differential equations are to be solved, because 
the operation of either differentiation or integration 
involves only the multiplication or division by jw. 
In addition, any periodic quantity can be expressed as 
an infinite series of complex exponentials rotating 
both positively and negatively at frequencies that are 
integral multiples of a fundamental frequency, «). 
This series, called a Fourier Series, pairs the exponen- 
tials in such a way that their summation is a real quan- 
tity. Sine and cosine functions can be expressed using 
only two terms of the Fourier Series. More compli- 
cated functions, such as the square wave and the 
sawtoothed wave, require an infinite number of terms 
for mathematical exactness, but usually a limited num- 
ber is satisfactory for engineering purposes. 

It is known from the elementary theory of linear 
differential equations that, when the input is described 
by a complex exponential, the steady-state output has the 
same angular frequency but has a different amplitude 
Qout)a and phase angle, (PA) (out)—that is, 


Q(out) = Qout)a eJlwt + (PA) (out) (6) 


It can be seen from a comparison of Eqs. (5) and (6) 
that the ratio of the steady-state output to the input 
is independent of time because the term e” is factor- 
able. This ratio is termed the performance function,* 
(PF)( )(4intout), SO that by definition 


(PF) \(4inZout) = (out) — (PA ‘(in)] 


(AR) )(ZinYout ) (4) 


In Eq. (7), the ratio of the output amplitude to the 
input amplitude is called the amplitude ratio, AR, and 
the difference in phase between the output and the in- 


* The performance function is equivalent to the frequency 
response or transfer fu-iction used by other authors. 
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put is called the phase angle, PA. The scheme oj 
subscripts established with the performance operator 
is used for designating the performance function for a 
particular application. The term vector performance 
operator, which was used in the previous article! has 
confused the distinction between an operator and 
function. The distinction can be explained by noting 
from Eqs. (5), (6), and (7) that, with a complex expo. 
nential input, 


= (PF) + (PA) inl 
and from Eqs. (1) and (5), 
Yiout) = \(4intout) + (9) 
so that 
e'"(PF)( )(aintout) = ¢P$( )(aintout) —(10) 


Eq. (10) shows that the product of the complex expo- 
nentia' e? and the performance function resultsfrom the 
operation, ¢P}, carried out on the complex exponential. 
The performance function affects the amplitude and 
phase of the complex exponential representing the 
output quantity but does not operate on the exponential 
in the truly mathematical sense. 

The performance function depends upon the angular 
frequency. The relationship for a linear differential 
equation with constant coefficients, determined by 
combining Eggs. (3), (5), and (8), is 


(PF)( )(¢invout) = 


be + ba + bs (jue)? + + 
do + ay (jw)! + az (jw)? + 


Ay ( jw)” 


The nth integer appearing in the denominator of Eq. 
(11) is usually greater than the mth integer shown in the 
numerator, in which case the performance function 
approaches zero as the angular frequency approaches 
an infinite value. 


(II) RELATION OF THE PERFORMANCE FUNCTION TO A 
GENERALIZED RESPONSE 


The dynamic characteristics of an operating com- 
ponent can be completely specified either in terms of 
its performance function or in terms of its response to 
a unit impulse (U7R). When the response to a unit 
impulse is known either by direct measurement or by 
analysis, the response to any input may be evaluated 
by use of the real convolution integral, a procedure 
commonly used to determine airplane structure loads 
due to prescribed control deflections. This integral,’ 
sometimes identified with Duhamel, states the relation 
between the input and output of a component or system 
as 


t 
Yiout)(t) = [(UTR)(t — )(@indout) dt (12) 
0 


or, substituting +r’ = ¢ — 7 into Eq. (12) and then drop- 
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ping the primes, 


t 

= any — )(@intout) (13) 

The response of a component can be obtained by con- 
volving the input quantity and the unit impulse re- 
sponse according to either Eq. (12) or (13). It should 
be noted that the output quantity can be determined at 
any time, ¢, and, hence, these equations can be used to 
study either transient or steady-state effects. Steady- 
state conditions are approached as the time becomes in- 
finitely large. 

Since this mathematical operation combines the in- 
put quantity with component characteristics to give 
the output quantity, it can be seen that 


t 
0 


Computing machines for operating on the input in the 
manner of Eq. (14) are in existence as, for example, the 
Cinema Integraph. An approximate arithmetic proce- 
dure proposed by Tustin‘ requires a summation based on 
the response to triangles of finite duration. 

Generally, in the study of closed-loop dynamic sys- 
tems, benefits result from a transformation that permits 
the use of a complex variable, rather than time, as the 
independent variable. There are a number of trans- 
formations, such as the Fourier Transform and the La- 
Place Transform, which can be utilized. The LaPlace 
Transform has the most general application, and a 
great many engineering treatments of its development 
are available. The LaPlace Transform is used in this 
section to show the relation between the performance 
function and the unit impulse response of a component 
or system, although it is a more powerful method than 
is needed for much of the analysis discussed in this pa- 
per. The special properties of the pulse allow the Four- 
ier Transform to be used in place of the LaPlace Trans- 
form in a large number of cases. 

The complex variable, s, of the LaPlace Transform 
has a real part, a, and an imaginary part w—that is, 

= a + jw—and the transform is defined as 


ILT][q(t)] = g(t)e~ “dt (15) 


For the Fourier Transform, the lower limit of the inte- 
gral in Eq. (15) is changed to — , and the complex 
variable is taken simply as jw. 

A theorem® that has important application states that 
convolution can be reduced to multiplication after trans- 
forming the dependent variables in the manner of Eq. 
(15). Application of this theorem to either Eq. (12) or 
(13) gives 

[LT] [Qcout) (¢)] 16 

[LT] [UTR] ) (4inJout [LT] (16) 

Eq. (16) can be used whenever the input quantity and 
the impulse response are transformable—that is, when- 


ever their insertion in Eq. (15) leads to a convergent 
integral. Since the complex variable, s, includes a fi- 
nite, real part that may be made large, a divergent input 
quantity is transformable provided it is of exponential 
order. Similar restrictions are placed on the transfor- 
mation of the unit impulse response. As a consequence, 
the response of any stable or unstable component hav- 
ing linear performance equations with constant coef- 
ficients can be determined from Eq. (16). 

The relationship of the performance function to the 
unit impulse can be established by the convolution in- 
tegral of Eq. (13). Because this integral is developed 
for general input functions, it is valid for the complex 
exponential input of Eq. (5). The transient response to 
a periodic input of this type is of little importance and is 
not taken into account in the original definition of the 
performance function. The steady-state response to the 
complex exponential input is: 


[UTR] (7) )(4in@out)@T 


= iin) + (PA x 
0 


[(UTR)(t)]¢ ) dt (17) 
A comparison of Eqs. (8), (15), and (16) shows that 


(PF); )(Zindout) = [LT] [UTR] )(Gin@out) | s = jw 
[LT] [Q(out) (¢)] 
[LT] (¢)] S=jw (18) 


When a mathematical expression for the unit impulse 
response exists, the performance function can be deter- 
mined from its transform by substituting jw for the 
complex variable, s, as indicated in Eq. (18). 

However, it is usually possible to short-cut the La- 
Place Transform when the performance function is 
evaluated experimentally. For the case of physically 
realizable components, the output quantity is zero prior 
to ¢ = 0, when the pulse occurs at some positive time. 
Hence, the lower limit of integration implied by Eq. 
(18) may be replaced by — © instead of zero without 
loss of generality. Then, since the condition s = jw is 
required, 


(19) 


This equation shows that the performance function 
equals the ratio of the Fourier Transform of the output 
quantity to the Fourier Transform of the input quan- 
tity. Eq. (19) may be used to obtain values for the per- 
formance function at any angular frequency, , pro- 
vided the integrals contained in the equation are con- 
vergent. 

The performance function carries the dimensions of 
the input and the output quantities. It is convenient 
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to define a parameter, independent of frequency, which 
carries these dimensions. This parameter, named the 
sensitivity S()(gingour), equals the value of the per- 
formance function at zero frequency when such a value 
exists, so that from Eq. (19) 


S( )(4in%out) = (out) (t) an/ f (t) dt (20) 


The sensitivity can be determined from an experimental 
response record by dividing the area under the output 
curve by the area under the input curve. The sensitiv- 
ity equals the ratio of bo to ao, the coefficients contained 
in Eq. (11). 

The sensitivity of certain components may have a 
value that is inconvenient to use. The sensitivity of 
the longitudinal response of the airplane is not usually 
considered for autopilot design because it includes the 
effect of both the long-period and short-period modes. 
In such cases, a reference sensitivity S()(¢ingout)(ref.) is 
defined which carries the dimensions of the input and 
the output quantities. The value of the reference sensi- 
tivity is selected for maximum convenience. 

If the input and output are known as analytical math- 
ematical functions, the performance function may be 
evaluated by direct application of Eq. (19). However, 
frequently the input and output are available only as 
recorded plots on oscillograph paper. The problem, 
then, is one of approximating these curves so that an 
evaluation of the performance function can be obtained. 


(III) DETERMINATION OF THE PERFORMANCE 
FUNCTION FROM PULSE RESPONSE DATA 


In Section (II) the performance function is shown to 
be the ratio of the Fourier Transforms of the output and 
input quantities. Special computing machines for 
taking the exact Fourier Transformer may not be avail- 
able, in which case methods presented in this and the 
following part of the paper can be employed. It is felt 
that, even when special computing machines are acces- 
sible, the methods outlined herein may prove more 
economical and nearly as accurate. 

A convenient method of approximating any curve® is 
by the use of a series of pulses occurring at definite time 
intervals. A useful pulse for this purpose is the unit 
triangular pulse [(U7)(¢)], defined in Fig. 1, which 
leads directly to approximating any curve by a series of 
chords of the curve. As shown, the pulse occurs during 
the interval from — Ar to Ar and is zero at all other 
times. The pulse attains its maximum value at zero 
time and, consequently, is symmetrical about the ori- 
gin. Its height is unity, and therefore the value of the 
pulse at any time is nondimensional. 

Shown in Fig. 1 with the unit triangular pulse is a 
similar but more general type called simply a triangular 
pulse [(7)(¢)]. Like the former, it has a time duration 
of 2 Ar and is symmetrical about the origin. Unlike the 
unit triangular pulse, its height, 7,, is not fixed but can 


THE UNIT TRIANGULAR 
~ PULSE HAS UNIT 
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BY DEFINITION, THE UNIT AND GENERAL PULSES ARE 
SYMMETRICAL AND HAVE A COMMON TIME DURATION 
OF 2 AT 

CONSEQUENTLY 


T(t-T) =T, [(UT)(t-T)| 


Fic. 1. Definition of unit triangular pulse and general triangular 
pulse. 
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Fic. 2. The approximation of a time function by a series of 
triangular pulses. 


be selected to give the best possible approximation of a 
given quantity. 

When a triangular pulse occurs not at time zero but 
at any time, 7, the unit triangular pulse is designated 
as [(UT)(t — 1)] and the general triangular pulse is 
[((T)(t — 7)]. Since the two pulses differ only in 
height, they are related by the expression 


((T)(¢ — 7)] = T.[(UT)(t — 7)] (21) 
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Fig. 2 shows how the triangular pulse is used to ap- 
proximate a general quantity g(t) which varies with 
time. The initial time 7» is fixed at a time when the 
quantity q(t) is zero, and the later instants of time 7, 
7», 73 are arrived at by successive increments of the con- 
stant time interval Ar. The corresponding values of 
g(t) are q(71), q(72), g(73). In terms of the correspond- 
ing unit triangular pulses, the function q(¢) is given at 
time ¢ by: 

q(t) = ((UT) (t — 71)] + @(72) (t — 72)] + 

(22 

The sum of two linear functions of some independent 
variable is itself a linear function of the variable. Hence, 
for example, the series of triangles in the region between 
7, and 7; is the straight line joining the ordinates g(72) 
and g(r3), since one of the linear functions occurring in 
this interval is zero at 72 and the other is zero at 73. By 
applying a similar reasoning to the remaining time in- 
tervals, it can be seen that Eq. (22) expresses q(t) as the 
sum of the successive chords. It is clear that the accu- 
racy of the approximation can be made arbitrarily good 
by reducing the time duration of the pulse Ar. 

Now the approximate Fourier Transform of the func- 
tion of time can be written as the summation of the 
Fourier Transform of each of the triangular pulses in 
the time series. In order to simplify the resulting sum- 
mation, reference is made to a Fourier Transform theo- 
rem. This theorem can be stated as follows: “If the 
Fourier Transform F(w) of a function of time, f(/), exists, 
then the Fourier Transform of f(t — 7) is e~ 7*’F(w).” 
A delay, 7, in the time domain is equivalent to a 
phase-angle lag in the frequency domain equal to the 
product of the time delay and the angular frequency. 
In terms of a unit triangle, 


[FT][(UT)(t — = (23) 
The series of triangular pulses used to approximate 
the quantity, g(t), are evenly spaced along the time 
scale, as shown in Fig. 2. Since the first matching tri- 
angle is delayed by 7, + Ar and the time interval be- 
tween successive triangles is A7, the mth triangle is 
delayed 7) + nAr with respect to zero time. Conse- 
quently, transformation of Eq. (22), followed by sub- 
stitution of Eq. (23), yields 
[FT] [g(t)] = + + + 
+. je (24) 


The exponential terms appearing in the series of Eq. 
(24) have magnitudes that depend upon the value of the 
quantity q(t) at 71, 72, 73, etc. The vectors are rotated 
negatively from the real axis of the complex plane by 
an angle proportional to the integer, m, the time incre- 
ment, Av, and the angular frequency, w. Although it 
might appear that any angular frequency can be in- 
serted in Eq. (24), the approximation is not satisfactory 
when w Ar is in excess of 7/2 rads. If this restriction 
does not allow a high enough angular frequency, the 
interval Av must be made smaller. 
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The Fourier Transform of the unit triangle must also 
be included in Eq. (24) when the complete transform of 
the quantity is desired. This transformation is dis- 
cussed in Appendix A, and is plotted as a function of 
wAr in Fig. 16. Since the unit triangular pulse, 
(UT)(t), is symmetrical about the origin, it is an even 
function, and its transform has zero phase angle. 

When Fourier Transforms are used to determine the 
performance function, it is only necessary to know the 
transform of the output relative to the input. Thus, 
when the same 7») and Ar are chosen for each quantity, 
the combination of Eqs. (19) and (24) gives 

The summations of Eq. (25) can be effected either (1) by 
direct graphical addition of the complex quantities or 
(2) by resolving the complex quantities into their real 
and imaginary parts before carrying out the addition.* 
A Fourier synthesizer based on method (2) has been 
built for this purpose. The synthesizer consists of elec- 
tromagnetic resolvers, the rotors of which are geared in 
the ratio 1:2:3: .. Potentiometers are used to in- 
troduce the correct amplitudes, and the outputs are 
summed electrically and recorded. The details of the 
synthesizer are given in Appendix B. 

An example is included to illustrate the direct graphi- 
cal addition of the complex quantities. A critically 
damped, second-order system is selected for this pur- 
pose. A critically damped system has a damping ratio, 
DR (ratio of actual damping to that required for critical 
damping), of 1. One other parameter is needed to de- 
fine the dynamic characteristics of the second-order sys- 
tem. The undamped natural frequency, w,, is a conve- 
nient parameter to use. The undamped natural fre- 
quency is the frequency at which the system would os- 
cillate in radians per unit time if, after removing its 
damping, the system were excited by an impulse. For 
this example, the undamped natural frequency is 7/10 
rad. per sec. 

The response of this second-order system to a unit 
triangle is shown in Fig. 3. The time duration of the 
pulse is 2 sec., and the resulting response returns to 
equilibrium in about 25 sec. The maximum excursion 
of the response is only 11.4 per cent of the input, but, 
because the static sensitivity of this system is unity, 
the time integrals of the input and the output are equal. 

The determination of the performance function by 
graphical summation is shown in Fig. 4. Values for the 
addition are obtained from Fig. 3 using a time interval, 
Ar, of 1 sec. Consequently, the input and output ordi- 
nates used in Eq. (25), are the values of these quantities 
at 1, 2,3,...sec. However, since the input quantity 
is zero at 2 sec. and all subsequent times, Eq. (25) re- 
duces in this example to: 


* This latter method is well adapted to computing machine 
operations since the functions cos(mwAr) and sin(mwAr) may 


be tabulated on card files. 
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a triangular pulse. 


—jnwdr 


(PF), )(4intout) = 
a= 
The upper limit of the summation is the 25th ordinate in 

as much as the pulse response returns to zero in 25 sec. 
The evaluation of the performance function is shown 
for a frequency of seven-tenths of the undamped na- 
tural frequency. Consequently, each consecutive term 
in the summation is rotated 0.077 rad. or 12.6° with re- 
spect to its preceding term. The angular orientation of 
the terms up to and including the 25th are laid off on the 
rose of Fig. 4. The values of the output quantity at the 
times 7, are read directly from the curve of Fig. 3. It 
can be seen that the summation has a magnitude of 0.67 

and a phase lag of 82.6°, so that 


(PF)( )(dintout) = [0.67 e~ = 0.67 e 


Hence, at a frequency of seven-tenths of the undamped 
natural frequency, the amplitude ratio of a critically 
damped second order system is 0.67 and the phase angle 
is —70° (the minus sign indicating a lag). 

Nondimensional plots of the amplitude ratio and the 
phase angle of a second-order system are presented in 
Fig. 14. The frequency is expressed nondimensionally 
as the ratio of the forced frequency, w, and the natural 
frequency of the system, w,, making 


(26) 
It can be seen that the result of the graphical summation 


of Fig. 4 agrees with the generalized performance func- 
tion data of Fig. 14. 


FR = frequency ratio = w/w, 


(IV) EXTENSION OF THE PULSE METHOD TO LIGHTLY 


DAMPED SYSTEM 


The approximate method for determining the per- 
formance function outlined in Section (III) is suited to 
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Fic. 4. 


components having a well damped or slightly oscillatory 
response. Direct application of this method to a com- 
ponent having a markedly oscillatory response becomes 
prohibitively tedious, inasmuch as many terms are re- 


quired before the series of Eq. (25) converges. An ex- 
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ample of this type is the lateral response of an airplane 
that has an oscillatory mode often referred to as ‘‘Dutch 
Roll.” When a highly oscillatory mode occurs, its pe- 
riod of oscillation and rate of decay can usually be 
measured directly. Hence, an equation for the oscilla- 
tory mode Qiout)(ose) (¢) can be written from which the 
Fourier Transform can be obtained directly. The re- 
maining portion of the response Qout)(rem.) (/) can be 
transformed by the approximate method previously 
demonstrated. The validity of this procedure is impli- 
cit in the assumption that the performance of the com- 
ponent can be represented by a linear differential equa- 
tion, so that 


(out) (t) = Q(out)(osc.) (t) + (out) (rem.) (t) (27) 
and, from Eq. (19), 


[FT] (ose) 
[FT] (¢)] 
[FT] [Qcouty(rem) 
[FT] 


During the time interval in which the input pulse is 
applied, the system has a forced response and a tran- 
sient response. However, the instant the inputreturns to 
zero, the output is a transient that is completely depend- 
ent upon the characteristics of the system. The tran- 
sient is the sum of exponentials of the form e”, where 7 
is a root of the system characteristic equation. When 
any of the roots have positive real parts, the system is 
divergent and unstable; when all roots have negative 
real parts, the transient is stable; when roots with imag- 
inary parts exist, they appear in conjugate pairs. Such 
a pair of exponentials gives an oscillatory mode, which 
can be written for > 0 


(PF): 


—(DR)wpt 
Q(out) (ose.) (t) { ja @ ‘ } x 


sin — (DR)% + (PA) (29) 


It can be seen that the real part of the roots leading to 
the oscillatory mode of Eq. (29) is equal to —(DR)w, 
and the imaginary parts are +0,V1 — (DR)*. The 
damping ratio DR and the angular undamped natural 
frequency w, are convenient parameters even though 
the system is third-order or higher. By using these 
parameters, nondimensional plots originally intended 
for second-order systems can be utilized directly. 

One such plot permits the damping ratio to be found 
for a particular transient. The heights of successive 
peaks are measured, and the transient peak ratios can 
be determined. When the mth peak is compared with a 
reference peak, m,, the transient peak ratio, [(7P)R], 
according to Eq. (29), is 


((TP)R] = — (DR)*) (30) 
Consequently, the ratio is independent of the natural 


frequency, depending only upon the damping ratio 
and the difference in the order of the peaks, m — m, 


(m — m, = 1 for adjacent peaks). The curves of Fig. 5, 
which are based on Eq. (30), are useful when evaluating 
the damping ratio. 

The discussion of the previous paragraph presupposes 
a datum from which the peak heights can be ascertained. 
With a second-order system, the datum is the final 
steady-state response, but with a higher-order system, it 
must be estimated. For example, the datum shown in 
Fig. 6 is the average of the upper and lower envelopes of 
the oscillatory mode. Since the ratio of the heights of 
adjacent peaks is 0.85, the damping ratio in this case is 
0.05. 

The undamped natural frequency is calculated from 
the damping ratio and the measured period of oscilla- 
tion, T\osc.). A sinusoid repeats every 2 7 rad.; hence, 
from Eq. (29), 


= 24 /T 1 — (DR)? (31) 


In the example of Fig. 6, the period is § sec., so that the 
undamped natural frequency becomes 0.78 rad. per sec. 

The oscillatory mode is seldom zero initially (when 
t = 0), especially since the actuating pulse is of finite 
duration. For this reason, the phase angle of the oscilla- 
tion, (PA) osc., must be determined so that the time de- 
lay is properly taken into account. For example, in the 
case shown in Fig. 6, the phase angle is 44°, which is 
equivalent to a time delay of approximately one-eighth 
of the period of oscillation. Once the damping ratio, 
natural frequency, and phase angle have been found, the 
only unknown in Eq. (29) is the amplitude of the oscilla- 
tion, Qiout)(osea- The amplitude can be calculated so 
that Eq. (29) matches the response at any time ¢ and, 
hence, at every time /. 

It is convenient to express the Fourier Transform of 
the oscillatory mode as a function of the frequency ratio 
defined in Eq. (26). Then, substitution of Eq. (29) into 
Eq. (15) with s replaced by jw gives 


JS Tout) (ose Sin (PA 
|F7 | [@cout) (ose.) (t) = f x 


— (FR)? + j2(FR)(DR)§ 
where 
(CT) ose. = a characteristic time of the oscillatory mode 
h istic ti f ill 
Sin (PA ose. 
wnlV | — (DR)? (PA Jose. + (DR) sin (PA 
(33) 


The actual evaluation of the frequency function of Eq. 
(32) is facilitated by the use of nondimensional plots of 
the amplitude ratio and phase angle, as given in Figs. 14 
and 15. 

The performance function equals the sum of the 
Fourier Transforms of the oscillatory mode and the re- 
mainder function, divided by the Fourier Transform of 
the input quantity, in conformity with Eq. (28). The 
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Fic. 6. The triangular pulse response of a system with a highly oscillatory mode. 


remainder function, called the output remainder quan- 
tity, Jrout)(rem,) (¢), equals the difference between the to- 
tal output quantity and the oscillatory mode. Fig. 7 
shows the triangular pulse response of Fig. 6 separated 
into the oscillatory mode and the remainder. 


The Fourier Transform of the remainder is obtained 
by the approximate method of Section (III), excepting 
in one particular. The remainder is zero prior to time 
ro, but at time 7 it has a finite value equal in magnitude 
but opposite in sign to the oscillatory mode. The re- 
sulting instantaneous change in magnitude cannot be 


matched by a triangular pulse that is isosceles in shape; 
consequently, a right triangular pulse is required.* 


* Experience gained recently on the lateral response indicates 
that it may be simpler and more accurate to commence the oscil- 
latory mode at a time, ¢p, when its value is zero [gyoutyose.) (4) = 
QO]. Then the right triangle is unnecessary, the phase angle 
(PA)osc. is eliminated, and the remainder re- 
mains equal to the output quantity, gour(¢), minus the oscil- 
latory mode for ¢ > ¢,; but for ¢ < ¢,, the remainder and the 
output are equal. A study of techniques for testing lightly 
damped systems was conducted by Comdrs. Bain, Michaelis, 


and Wootton and is discussed in their Master’s Thesis entitled, 
“Investigation of Aircraft Lateral Motion Performance Func- 
tion by Pulse Technique,’’ M.I.T., September, 1949. 
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The right triangular pulse attains its maximum value 
instantaneously and then decreases linearly with time, 
reaching zero after time Ar has elapsed. Consequently, 
the right triangular pulse has a duration half that of the 
regular isosceles triangular pulse described in Section 
(III). The unit right triangular pulse is shown in Fig. 
17 of Appendix A, and the Fourier Transform of the unit 
right triangular pulse is presented in Fig. 18 of the same 
Appendix. Since the unit right triangular pulse, 
[((URT)(t)], is unsymmetrical, its transform includes 
both an amplitude and a phase angle. The phase angle 
approaches 90° as the frequency becomes infinitely 
great. 

The shaded part of Fig. 7 is nearly equal to the prod- 
uct of the ordinate, [Qiout)(rem.)(70)], and a unit right 
triangular pulse. As a result, the total remainder can 
be expressed 


Q(out)(rem.) (t) [9cout)(rem.) (70) ] [(URT) (t)] + 


[cout)(rem.)(7n) ] [((UT) Tn) (34) 


Then, with the input quantity represented as a series of 
isosceles triangular pulses, substitution of Eqs. (32) 
and (34) into Eq. (28) gives 


(PF); )(Gin%out) = 


[FT] [dout) (ose.) ) (t) [Yiout)(rem )(70) 


—jnwAt oo 

(tn) Je le=1,2,... 
Although the complete equation for the performance 
function appears formidable, the work is reduced to the 
determination of the constants for the oscillatory mode; 
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3. 8. Performance function obtained from the triangular 
pulse response of Fig. 6. 


the summation of two exponential series; the selection 
of proper values from Figs. 14, 15, 16, and 18; and the 
addition of the complex quantities of Eq. (35). This 
procedure for the example of Figs. 6 and 7 gives the per- 
formance function of Fig. 8, represented in terms of am- 
plitude and phase angle. The amplitude and phase 
angle are shown as functions of the frequency ratio, 
which is the ratio of the angular frequency to the un- 
damped natural frequency of the oscillatory mode. 
An amplitude peak occurs when these two frequencies 
are nearly equal. If a component with the char- 
acteristics of this example were included in a closed- 
loop system, this peak, unless compensated, would 
tend to make the system unstable. 

Some dynamic systems include a component that is 
inherently unstable in one of its modes. For example, 
the characteristic equation of certain aircraft contains a 
divergence-producing root, and the resulting motion of 
the aircraft is usually referred to as the “spiral diver- 
gence.’ It might seem that such aircraft could not 
be included in an automatically controlled system with- 
out producing instability. This is not the case; in fact, 
in certain systems, unstable components and subloops 
may be intentionally created to produce improved over- 
all system performance. 

Since slightly unstable components are not only to be 
expected, but on occasions to be preferred, adequate 
methods for determining their performance functions 
experimentally are required. The performance function 
is a multiplying factor that relates the output response 
to a complex exponential input after steady-state condt- 
tions are reached. However, an unstable operating com- 
ponent never attains steady-state unless the component 
is imbedded in a stabilizing closed loop. When this is 
done, the entire loop can be excited, and the resulting 
input and output of the unstable component can be 
measured. Flight tests have demonstrated the effec- 
tiveness of this technique for steady-state sinusoidal in- 
puts. It is suggested that the pulse inputs might be 
applied in place of sinusoidal inputs. The artificial sta- 
bilization would bring the aircraft back to equilibrium, 
and, hence, the Fourier Transform of the input and 
output would exist and the performance function could 
be evaluated either from Eq. (19) or from the ap- 
proximate relationship of Eq. (25). 


(V) Fricutr Test 


The most attractive features of the pulse method for 
obtaining airplane performance functions are its econ- 
omy in equipment, installation time, and required flight 
time. In all probability the greatest saving is real- 
ized by having the human pilot initiate the forcing 
function. This eliminates an automatic pilot installa- 
tion, which would be required if the frequency response 
characteristics were to be obtained directly. Moreover, 
experience has shown that the flying time can be cut to 
one-third of that required with direct oscillation of the 
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control surfaces. As a direct result of the economy in 
flying time, sufficient flight-test data can readily be ob- 
tained to cross-check and resolve any discrepancies that 
may appear in the test results. 

The pulse method is in many ways similar to both the 
short-period longitudinal test and the lateral rudder- 
kick test, which are standard flight tests to obtain quali- 
tative information as to the transient response of an air- 
craft. The principal difference, however, is that the 
latter tests are normally made for the control-free air- 
plane, whereas the pulse method is used to obtain the 
control-fixed oscillation about the original equilibrium 
flight configuration. 

In order to ensure oscillation about the original equi- 
librium condition, an equilibrium control position stop 
must be provided. Theoretically, this stop alone is suf- 
ficient ; however, unless a stop is,also provided, to.con- 
trol the magnitude of the pulse, it is extremely difficult 
to set the oscillograph sensitivities to give reasonable re- 
sponse traces. 
test dictate the form of the control surface jig. It may 
be a simple, double chain arrangement, or a more elabo- 
rate mechanism like that shown in Fig. 9. The stop 
mechanism shown in this figure was designed by C. E. 
Muzzey, of the Flight Test Section of the M.I.T. In- 
strumentation Laboratory. The mechanism consists of 
a removable rack that rides ona sliding rack. The slid- 
ing rack is pivoted at the control column and is free to 
slide (except when constrained as noted below) in a 
slide pivoted to the air frame. This slide contains a 
hinged gate When the gate is closed against the slide 
block, it provides a double stop that limits the control 
column motion. When the gate is open, the control 
column is completely free. The purpose of the hinged 
gate arrangement is to allow the pilot to completely 
free the airplane so that he can control the long-period 
oscillation that becomes predominant after the short- 
period oscillation has subsided. 

In making a pulse test, the pilot executes the follow- 
ing simple steps: 

(1) He trims the airplane in equilibrium flight at the 
test speed and pressure altitude. 

(2) When equilibrium flight has been obtained, he 
moves the trim stop into position, taking care not to dis- 
turb the equilibrium flight condition as the trim stop is 
moved into place. 

(3) He moves the control rapidly against the oppo- 
site stop and returns it to the trim stop in accordance 
with the pulse duration time desired. Experience has 
shown that a control pulse of 1° sec. will initiate a rea- 
sonable transient response in an airplane of medium size. 

(4) He holds the control securely against the trim 
stop until the oscillation has subsided. 

The following variables were recorded during each 
transient on an oscillograph, shown in Fig. 10: (1) 
normal acceleration, (2) rate of pitch angle, and (3) 
elevator position. The measuring system used in this 
flight test program was largely composed of components 
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BACK STOP. (ALSO SERVES AS 
DISENGAGING LEVER 
FRONT STOP is BEHIND 
INSTRUMENT PANEL) 


SLIDE PIVOTED To 
AIRFRAME 


VERNIER TRIM 


ADJUSTMENT REMOVABLE 
RACK 


RACK PIVOTED TO CONTROL 
COLUMN 


CONTROL WHEEL 


Fic. 9. Photograph of stop mechanism for control column, 


previously designed for other applications by the In- 
strumentation Laboratory of M.I.T. 

The normal accelerometer used for measuring the 
normal acceleration includes a temperature-controlled 
viscous damper. The damper temperature is regulated 
by a thermoswitch and is held to close temperature tol- 
erance. For these tests, the damping ratio was 1.0 and 
the undamped natural frequency 10 cycles per sec. The 
uncertainty of this accelerometer is + 0.005 gravity. 

The rate-of-pitch indicator is a single degree of free- 
dom gyro computer, which was modified for these 
tests to obtain a damping ratio of 0.7. The undamped 
natural frequency of this instrument is 5 cycles per sec.; 
its uncertainty is +1 millirad. per sec. 

The elevator displacement is indicated by the electro- 
magnetic type of signal generator called a microsyn. 
This instrument is.a variable-transformer type designed 
by Dr. R. K. Mueller, of the Instrumentation Labora- 
tory M.1.T., for fire control applications where the an- 
gular displacement of a shaft must be measured within 
'/, millirad. 


The normal accelerometer and the rate-of-pitch gyro — 


are also furnished with a microsyn signal generator. 
The output from each of the signal generators was am- 
plified before being recorded. Filters were used on the 
first flight tests but were found to be unnecessary. The 
oscillograph record shown in Fig. 11 is unfiltered data. 
Both static calibrations and dynamic corrections must 
be applied to all test results to correct for instrumenta- 
tion errors. The dynamic corrections can be obtained 
from direct sinusoidal oscillation of the instrument in 
the laboratory. The frequency response characteris- 
tics are more conveniently determined after installation 
by applying either a step function or pulse to the instru- 
ment, by recording the transient response, and by com- 
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ACCELEROMETER LOCATED FORWARD 
AT CENTER OF GRAVITY OF AIRPLANE 


RATE- OF-PITCH 
INDICATOR 


GALVANOMETER 
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ACCELERATION AND 
ELEVATOR DISPLACE 
MENT SIGNALS. 


ATTENUATORS FOR 
THREE RECORDING 
CHANNELS 
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OSCILLOGRAPH 


Fic. 10. Photograph of flight-test instrumentation located in 
aft compartment. 


~ RATE OF PITCH@, 


| 
NOSE UP 
044Qq 
UPWARD 
ACCELERATI 
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~—+ ELEVATOR DISPLACEMENT, 
TIME —= 
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FLIGHT 12 10/13/48 
175 mph LAS ALTITUDE 8900 FT 
C.G. AT 21% MAC APPROX WT 26,000 LBS 


MOMENT OF INERTIA 64,000 SLUG FT* 


B-25J NO 328 


Fic. 11. Oscillograph record showing rate of pitch and nor- 
mal acceleration response of an aircraft to an applied elevator 


pulse displacement. 


puting the frequency response. The means for provid- 
ing a sinusoidal input are not readily available, whereas 
a pulse or step function can be easily applied. 


Static calibrations of the measuring system caused 
more difficulty than the dynamic calibrations. Experi- 
ence gained in this flight-test program has shown that 
the only way to be certain of static calibrations is to use 
frequent in-flight calibrations. In this manner any drift 
or change in the sensitivity of a channel is determined, 
and the proper corrections can be applied to the meas- 
ured data. 
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airplane obtained from elevator pulse displacement with perform- 
ance function obtained from theoretical considerations. 
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Fic. 13. Comparison of normal acceleration performance func- 
tions of airplane obtained from elevation pulse displacement with 
performance function obtained from theoretical considerations. 


(VI) CORRELATION OF AIRCRAFT PERFORMANCE FUNC- 
TIONS OBTAINED EXPERIMENTALLY AND THEORETICALLY 


Correlation of the results obtained by pulse data 
with possible results obtained by direct sinusoidal oscil- 
lations of the elevator was not attempted because the 
automatic pilot installation and subsequent flying 
would have extended the program beyond reasonable 
limits. In place of direct sinusoidal oscillation data, it 
was decided to use the theoretical performance equation 
of the airplane and, as far as possible, use only static 
flight-test data in determining the stability derivatives. 
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The theoretical performance functions*® required 
for this discussion are derived in Appendix C. Eq. 
(C-9) for rate-of-pitch response and Eq. (C-19) for 
normal acceleration response are reproduced here: 


1 + j(FR)on(CT)a (sa) 


(PP )a = Sa omy (eet) 1— (FR)? + j2(DR)(FR) (C-9) 
(PF) a (on) = Sean) + Saceny(rem) X 

1 + j(FR)on(CT) acon) (C-19) 


1 — (FR)? + j2(DR)(FR) 


When written in this manner, the frequency-variant 
portions of the performance functions have a character- 
istic numerator and denominator. Consequently, sepa- 
rate nondimensional curves can be plotted of each 
function and the curves used to compute the aircraft 
performance functions. These nondimensional curves 
are presented in Figs. 14 and 15. To use the curves, 
the characteristic parameters w,, (CT), ), and (DR) 
must be calculated from the stability derivatives ac- 
cording to their equations given in Appendix C. The 
numerator is then evaluated as a function of the par- 
ameter w,(CT)4,); the denominator, as a function of 
the damping ratio (DR). The aircraft performance 
function is then found by a straight algebraic process 
involving complex numbers. These nondimensional 
curves are extremely useful to determine quickly the 
effects of changes in the center of gravity, air speed, and 
altitude upon the dynamic performance of the airplane. 
The smooth curves of Figs. 12 and 13 were calculated 
from these equations for the values of the stability 
derviatives given in Appendix C. 

The normal acceleration performance function [Eq. 
(C-19)] is especially interesting because it contains a 
separate sensitivity, Sz(5,). The inverse transform of 
this sensitivity gives an impulse in the time domain. 
However, the inverse transform of the product of this 
sensitivity and the direct transform of a finite pulse are 
equal to the pulse in the time domain scaled by this 
sensitivity. Therefore, the system response to a finite 
pulse will include a similar finite pulse occurring at the 
same time as the input pulse. Hence, one would ex- 
pect to find a “‘pip”’ in the normal acceleration response. 
This has been observed experimentally. Although Fig. 
11 shows only a small initial negative bulge in the nor- 
mal acceleration response, similar records for shorter 
pulse time durations have exhibited sharp ‘“‘pips.” 
Physically, this ‘‘pip’”’ has a simple explanation: In or- 
der to rotate the wing so as to achieve an increased an- 
gle of attack, it is necessary for the tail to be depressed 
by the elevator. The first effect of this action is to 
decrease the total lift of the airplane because of the loss 
of lift on the tail. From Eqs. (C-19) and (C-20) it can 
be seen that the area of this ‘‘pip’”’ is directly propor- 
tional to Z; and, hence, represents the loss of total lift 
on the airplane caused by the elevator deflection. 

The experimental data points plotted on Figs. 12 and 
13 were computed from Fig. 11 and corrected for instru- 
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Fic. 15. Relationship between the performance function 
[1 + j(FR)w,(CT)], and the frequency ratio. 
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Fic. 18. Fourier transform of unit right triangular pulse in 
: terms of amplitude and phase angle. 
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ment errors. The performance functions were reduced 
by the Fourier Analyzer described in Appendix B and 
were checked by the graphical method described in 
Section (IIT). 

The degree of correlation between the experimental 
points and the smooth curves should not be interpreted 
as a measure of the accuracy of predicting the airplane 
performance functions from theoretical equations. The 
major drawback to calculating the performance func- 
tion is that, in general, the aerodynamic parameters will 
not be known accurately enough to give good results 
and the equations themselves will have to be extended 
to account for aeroelastic and Mach Number effects. 
However, where the performance functions have already 
been determined experimentally for a number of flight 
configurations, fitted stability derivatives can be used 
to interpolate the flight-test configurations. 


CONCLUSIONS 


(1) The pulse technique of determining airplane per- 
formance functions is a feasible flight-test method. 

(2) The use of the pulse method for determining the 
airplane performance functions will save in equipment, 
installation time, and required flying time. 

(3) The instrumentation as used in the flight-test 
program has low enough uncertainty levels and suf- 
ficient sensitivity to determine the performance vari- 
ables to the degree of accuracy required for dynamic 
analysis. 

(4) The airplane theoretical performance equations 
are especially useful when interpolating dynamic 
flight-test data. 

(5) The Fourier Synthesizer and Analyzer has a defi- 
nite advantage over any of the conventional harmonic 
analyzers for two reasons: (a) Transforms can be 
made directly from oscillograph records, and, hence, 
data need not be plotted to any special scale; (b) the 
transforms are obtained as continuous functions of fre- 
quency instead of selected frequency points. 

(6) The approximation of a time function by a series 
of short chords will give accurate transforms at low fre- 
quencies and can be extended to any frequency desired 
by decreasing the time increment Ar. 


APPENDIX A——-FOURIER TRANSFORM OF UNIT 
TRIANGULAR PULSES 


In Section (IV), a method is described for approxi- 
mating the Fourier Transform of a response that has a 
lightly damped oscillatory mode. This method requires 
the Fourier Transform of the unit triangular pulse and 
the unit right triangular pulse. It is convenient to have 
the Fourier Transforms’ of the unit triangular pulse and 
the unit right triangular pulse available in curve form 
so that their values need not be computed for the fre- 
quencies of interest. These curves and a brief explana- 
tion of their derivation are presented in this appendix. 
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The unit triangular pulse is symmetrical about the 
time origin and has a time duration 2 Ar, so that, prior 
to time — Ar and after time Ar, its magnitude is zero. 
According to the definition of the unit triangular pulse 
as described in this appendix and as shown in Fig. 1, 


(UT)(t) = 0, fort < — Ar 

(UT)(t) = 1+ (t/Ar), for —Ar<t <0 
(UT)(t) = 1—(t/Ar), forO <t< Ar 
(UT)(t) = 0, for Ar<t 


The use of these conditions for the unit triangular pulse 
permits its Fourier Transform to be expressed as 


ao 


[FT][(UT)()] -f dt = 


sin w Ar/2 


wAr/2 


Since the unit triangular pulse is an even function, its 
transform has an amplitude but no phase angle. The 
amplitude is plotted in Fig. 16 as a function of the non- 
dimensional frequency wAr. The amplitude of the 
transform has a value of Ar at zero frequency and is re- 
duced to zero when the nondimensional frequency 
equals 27 rad. In order to extend the useful range of 
frequencies, the time increment Ar can be decreased. 

A quantity with an abrupt change in value at the 
origin cannot be matched by a symmetrical triangular 
pulse. Consequently, the unit right triangular pulse 
shown in Fig. 17 is defined. The unit right triangular 
pulse is zero for all negative times, has a height of 
unity at the origin, and decreases to zero at Ar. This 
definition may be stated mathematically : 


(URT)(t) = 0, fort <0 
(URT)(t) = 1 — (t/Ar), forO <t < Ar 
(URT)(t) = 0, for Ar <t 


With the unit right triangular pulse defined in this 
fashion, its Fourier Transform becomes 


] sin wAr/2] 
Jw wAr/2 


The Fourier Transform of Eq. (A-2) contains both a 
real and an imaginary part. Hence, the transform of the 
unit right triangular pulse has an amplitude and a phase 
angle, both of which vary with frequency. Values for 
the amplitude and phase angle can be obtained from 
Fig. 18. 


APPENDIX B—THE FOURIER SYNTHESIZER 


The resolution of a series of complex exponentials 
into the sum of their real and imaginary parts is an op- 
eration frequently encountered in the process of trans- 
forming either from the time domain to the frequency 
domain or from the frequency domain to the time do- 
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Fic. 19. Photograph of Fourier synthesizer. 

main. When large amounts of data are to be reduced, 
the need for a specially constructed synthesizer becomes 
apparent. The Instrumentation Laboratory at M.LT, 
has built the electromechanical synthesizer shown in 
Fig. 19. This unit consists of an assembly of twelve 
synchros. The rotors of these synchros are geared to- 
gether either in the ratio 1:2:3.... or in the ratio 
1:3:5.... The second of these ratios is used exclu 
sively for inverse transformations involving the response 
of a component to an applied step function, The stator 
of a synchro can be set in any angular position to ae- 
count for the phase angle of the complex exponential 
term it is representing. The output from each synchro 
is the voltage induced across a stator winding. This 
voltage is attenuated by means of a Helipot potenti- 
ometer to give the desired magnitude of the complex 
exponential term. The attenuated voltages from the 
synchros are added electrically. The resulting voltage 
is amplified, rectified, and impressed on the galvanome- 
ter of a Sanborn Recorder. The recorder produces a 
trace on 2-in. paper. 

The use of the synthesizer for evaluating the series oi 
complex exponentials of Eq. (24) may be described in 
the following manner: The series of Eq. (24) can be 
written 


q-series= 
= q(t) Cos nw Ar — 


Q(T) sin mw Ar 


n=1,2,3... 


R, |q-series] + [q-series] (B-1) 


The amplitude and phase angle of the series can be 
readily determined, since 

[q-series] = [g-series],e’ (PA) [q-series| (B-2) 
Then 


[q-series]q = V {R, [q-series] 2 In {g-series] (B-3) 


and 
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and w= (PF) 4 (C-6) 
(PA) {q-series) = tan { In|g-series]/ R, [g-series] } (B-4) (PF) 4 sin) (C-7) 


The stators of the synchros are oriented so that, with 
the phase angle setting at zero, the voltage induced 
across each synchro secondary is proportional to cos ®), 
where () 4s the angular position of the rotor. With a 
gear ratioof 1:2:3... m ..., the rotor position can be 
related to the rotor position of the fundamental rotor, 
@, simply as n@,. However, the real part of the nth 
term of Eq. (B-1) depends upon cos ” wAr; conse- 
quently, the nondimensional frequency, w Ar, is equiva- 
lent to the angular position of the rotor representing the 
fundamental complex exponential. Since both the syn- 
thesizer and the recorder are driven at constant speed, 
the timing line of the recording paper can be calibroted 
in terms of angular frequency, w. The galvanometer 
current is proportional to the sum of the cosine volt- 
ages, each properly attenuated, and, consequently, the 
ordinate can be calibrated in terms of R,[g-series]. 

Recalling that sin @ = —cos (90 + @), it is seen 
that, by simply rotating each stator by 90°, the J [q- 
series] is obtained. It is anticipated that the synchros 
will be replaced by magnetic resolvers havitig secondary 
windings in space quadrature, so that eventually both 
the real and imaginary parts can be obtained simulta- 
neously. 


APPENDIX C—AIRCRAFT LONGITUDINAL RESPONSE 


In this appendix the equations of motion are de- 
veloped that are used to calculate the theoretical re- 
sponse curves for comparison with the experimentally 
determined response curves. Since experimental pulse 
data have been obtained to date primarily for the longi- 
tudinal response, only the longitudinal equations are 
considered. The derivation of the classical linearized 
equations of motion are not pertinent to this paper and, 
consequently, are abstracted from Volume V of Durands’ 
Aerodynamic Theory. From page 129 of Volume V, the 
longitudinal equations may be written as: 


a+ g@ = uX, + wX, + @MX 
w— uM = uZ, + wZ, + 6Z; (C-2) 


@ = + + + wM, + 6M, (C-3) 


It has been shown by Rea" that, for frequencies greater 
than 0.5 rad. per se -., the variation in air speed is neg- 
ligible. Hence, the above set of differential equations 
can be reduced from three to two—1.e., foru = u = 0, 


® — @Ma —wM; —wM,=6M, (C-3) 
When the elevator input displacement is expressed in 


the complex exponential form, 6,e7“, the solution of 
Eqs. (C-4) and (C-5) gives 


where 


+ I(FR)on(CT) 4 (sw) 
PF w) = S, w)(ref.), 
(PP) = + j2(DR)(FR) 


Ai (C-8) 


1 + j(FR) on(CT) 
1 — (FR)? + j 2(DR)(FR) 


(AR) acs) ei (C-9) 


(PF) (sin) = S4 (ai) (ref) 


In Eqs. (C-8) and (C-9) it has been found convenient 
to define the following characteristic parameters: 


= VZ.M@ — (C-10) 
(DR) = —(Mj + Ze + moM;) + 
2VZ,M; — (C-11) 


(FR) = w/e, (C-12) 
Sacaayret) = (ZsMo — — (C-13) 
(CT) acn = (Mst+ Z,Mz)/(MeZ, — MiZe) (C-14) 
Sacawyeet) = (toMy — M@Zs)/(ZeM@ — (C-15) 
(CT) ata) = Zs/(toM, — M@Zs) (C-16) 


Eq. (C-8), however, involves linear velocity along the 
Z axis, which is an inconvenient output quantity for 
experimental measurements. A more convenient 
quantity to measure in flight is normal acceleration, m,. 
Eqs. (C-8) and (C-9) can be combined to give a per- 
formance function in terms of normal acceleration, since 


= (w— uD) (C-17) 


making 


(jw) 
g 


= 

(PF) sn) (PF) 4 (aw) g (PF) aw) (C-18) 
After carrying out the indicated operations and sim- 

plifying, the performance function can be written 


(PF) 
1 WFR) wn(CT) on) 


E(én + A (dn) ( ) 1— (FR)? + J 2(FR)(DR) 


where 


u Ms | (C-21) 
— wM, 


Sa (én)(rem.) = 


CT 4 (5n) = ¢ T (C-22) 
In addition, 
S4 (an) (ref) = SE (sn) SA (sn)(rem.) 
= £2) SA (ref) (C-23) 
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Eqs. (C-9) and (C-19) can be used to calculate the 
aircraft theoretical performance functions. The theo- 
retical response curves in Section (VII) are calculated 
from these equations using the following numeral data: 


Physical Data 


I, = 64,000 slug ft.? (measured) 
W = 26,000 Ibs. 

S(wing) = 610 ft.? 

Ciwing) = 9.67 ft. 


Aerodynamic Data (175 m.p.h., 21 per cent M.A.C.) 


S(tait) = 182.5 ft.* 
ital) = 25.9 ft. 


dC 
(3) = 5.25 per rad., from flight test 
A 


da 
SM = static margin = —0.15, from flight test 
dCn 
or = —0.79 per radian, from flight test 
a 
dC, 
(““ ‘) n, = 3.9 per rad., computed 
da t 
d dC, 
— - 0.4, computed "= —13 per rad., computed 
da dé, 


Computed Stability Derivatives 
175 m.p.h., 21 per cent M.A.C., 10,000 ft., Pressure Altitude 


Sele 


din 
Z, = — = — 30.4 ft. Ss 2 


= —9.8 (1/sec.?) 


dem\SwCw 
M, = (“=)5-6 q = —0.0192 (1/ft. sec.) 
da 


qSw 
(4+) = —1,065 (1/sec.) 
da J 
(1/ft.) 
dC Sil? 
Me = — (1.2) = =1.72 (1/sec.) 
da t 
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Characteristic Parameters for Use in Eqs. (C-9) and (C-19) 


wy, = 2.78 rad./sec. from Eq. (C-10) 
(DR) = 0.6 from Eq. (C-11) 


deg. /sec 


SA(8®)(ref.) = —1.2 —— from Eq. (C-13) 
deg. 
(CT)A(sm) = 0.995 sec. from Eq. (C-14) 
SE(sn) = —0.0165 g/deg. from Eq. (C-20) ‘ 


SA(én)(rem.) = 0.229 g/deg. from Eq. (C-21) 
(CT) A (sn) = 0.010 sec. from Eq. (C-22) 
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SUMMARY 

Results available on lateral-stability derivatives for wing 
plan forms at supersonic speeds are summarized. The effects of 
aspect ratio and Mach Number are studied primarily, and sub- 
sonic values for the derivatives are included in the investigation 
of the Mach Number effects. Nonlinearities in the sideslip 
yariation of the yawing and rolling moments are considered 
briefly, and the effects of combining a wing and body are indicated 
for two of the derivatives. 


INTRODUCTION 


_ PRIMARY FIELD of application for stability de- 
rivatives is the calculation of the motion of an air- 
plane. A stability derivative, expressed in simple 
terms, is the ratio of the force or moment produced by 
motion along or around one of the airplane axes to the 
magnitude of the motion. These derivatives appear 
as coefficients in the classical set of six linear differential 
equations derived for the motion of the airplane on the 
basis of small disturbances. The application of these 
equations, three of which apply to longitudinal motion 
and three to lateral motion, to the study of the stability 
of the airplane motions is generally referred to as dy- 
namic-stability analysis in contrast to the more widely 
applied static-stability analysis. 

In the past 10 or 15 years, dynamic-stability analysis 
has become more and more popular as a tool for the air- 
craft designer and the research analyst. The interest 
in this type of analysis seems to be reaching a new 
zenith in the design and study of guided missiles. The 
stepup in speed of aircraft which has occurred in this 
10- to 15-year period does not materially affect the 
classical form of the equations of motion. The increase 
in speed, however, does have a significant effect on the 
values for the stability derivatives that appear as the 
coefficients of these equations. In fact, in the transition 
from subsonic to supersonic speeds, changes in the sign, 
as well as the magnitude, of the derivatives occur in 
some cases. 

It is the purpose of this paper to summarize the 
lateral-stability derivatives that have been calculated 
for wing plan forms traveling at supersonic speeds. 
The wings considered have no dihedral or twist. Longi- 
tudinal-stability derivatives will not be included in the 
summary, although a considerable amount of investiga- 
tion to obtain these derivatives has been carried out. 


Presented at the Aerodynamics Session, Annual Summer Meet- 
ing, I.A.S., Los Angeles, July 21-22, 1949. 
* Aeronautical Research Scientist. 
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The Theoretical Lateral-Stability Derivatives 
for Wings at Supersonic Speeds 


ARTHUR L. JONES* 
Ames Aeronautical Laboratory, N.A.C.A. 


No experimental comparisons will be made, inasmuch 
as the experimental data that exist on supersonic sta- 
bility derivatives are classified at the present time. 
The main objective of the summary will be to study the 
effects of aspect ratio and of Mach Number on the 
derivatives. Other results of interest will also be in- 
cluded, and, in view of the fact that this paper is a sum- 
mary and the results are of paramount interest, only a 
brief discussion of the methods of calculation will be 
given. 

There are nine lateral-stability derivatives. Roll- 
ing, yawing, and sideslip each contribute a side force, 
a rolling moment, and a yawing moment. For the wing 
alone, the magnitudes of the three side-force derivatives 
Cy,, Cy,, and Cyg are extremely small and their con- 
tribution to the motion of the airplane is, in general, 
insignificant. Thus, the only stability derivatives with 
which we shall be concerned are, in coefficient form: 


C,, = the damping in roll 

C,, = the yawing moment due to rolling 

C,, = the rolling moment due to yawing 

C,, = the damping in yaw 

Ciz = the rolling moment due to sideslip, or ‘the 
effective dihedral angle’”’ 

Cig = the yawing moment due to sideslip or “the 
directional stability.” 


SYMBOLS AND COEFFICIENTS 


A = aspect ratio (b?/S) 

b = span of wing measured normal to plane of symmetry 
B 

¢- = chord of wing root 

C: = rolling-moment coefficient (rolling moment/gSb) 

Cyn = yawing-moment coefficient (yawing moment/g.Sb) 

Cy = side-force coefficient (side force/gS) 


= 0C,/0(pb/2V) 
Ci, = OC,/d(rb/2V) 
Cig = OC,/d8 

Cap = OC,/0(pb/2V) 
Cap = OCn/d(rb/2V) 
= OC,/o0B 

Cy, = OCy/d(pb/2V) 
Cy, = OCy/d(rb/2V) 


Cyg = OCy/os 

l = overall longitudinal length of sweptback wing 

m = slope of right wing tip relative to plane of symmetry 
(see Fig. 1) 

M = free-stream Mach Number 

p = rate of roll, rad. per sec. 

gq = free-stream dynamic pressure 


q 
| 
re 
7 
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r = rate of yaw, rad. per sec. 

S = area of wing 

V = free-stream velocity 

a = angle of attack, rad. 

8 = sideslip angle (positive when sideslipping to right), deg. 
Subscript 

w = wing alone 


METHOD OF ANALYSIS 


Of the theoretical approaches that are available for 
the analysis of problems in supersonic fluid dynamics, 
the one that has seen the greatest application in the 
calculation of supersonic stability derivatives for wing 
plan forms is the linearized theory of compressible flow 
as applied to thin airfoils. This is a potential theory, 
and its application reduces the study of supersonic lift- 
ing-surface problems to a process of fitting potential 
solutions of the well-known wave equation of mathe- 
matical physics to the somewhat simplified boundary 
conditions of thin-airfoil theory. Once the potential 
solutions corresponding to the boundary conditions are 
known for a given wing plan form, the pressures acting 
normal to the plate and the pressures on the subsonic 
leading edges,* usually referred to as the edge suction 
pressures, can be determined. From these pressures 
the lift and drag forces and the various moments can be 
obtained by straightforward integrations. 

With regard to the accuracy of the calculated results, 
it should be noted that, although the theory is based 
on the assumptions of small angles of attack and cor- 
respondingly small perturbation velocities, the similarly 
restricted theory in subsonic flow has yielded surpris- 
ingly accurate results. Furthermore, comparisons 
that have been made, such as in reference 1, between 
the available supersonic experimental data on lift- 
curve slope and the theoretical values based on the lin- 
earized compressible-flow theory are encouraging. 
The lift-curve slope is, of course, a longitudinal deriva- 
tive. Furthermore, it is dependent on the total or 
integrated pressure and not the distribution of the 
pressure. This favorable comparison, therefore, is 
somewhat indirect, but, when it is coupled with the fact 
that the lateral derivatives considered herein are not 
very sensitive to small shifts in the spanwise centers of 
pressure (being dependent on fairly large lateral mo- 
ment arms), the indication is that the linearized theory 
should provide fairly good estimates of the potential- 
flow contribution to the lateral derivative, at least in- 
sofar as general trends and orders of magnitude are con- 
cerned. Certain of the viscous effects, such as the con- 
tribution of skin friction to C,,, can be evaluated in- 
dependently and added to the potential flow contribu- 
tions. Thus, it is possible that the present results will 


* In the terminology that has become associated with super- 
sonic flow analysis, the term ‘‘subsonic edge’’ refers to an edge 
swept behind the Mach cone, from which it follows that the ve- 
locity normal to the edge is subsonic; the term ‘‘supersonic edge”’ 
refers to an edge swept ahead of the Mach cone. 
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be sufficiently accurate to preclude widespread appli- 
cation of a more elaborate theoretical method or com. 
plete recourse to experimental results. In an instance 
where a rather precise value for a derivative is desired, 
however, an experimental result or a more compre. 
hensive theoretical investigation may be required, 


As previously mentioned, one of the first steps in the 
solution of the general problem of lifting-surface analy. 
sis is the establishment of the boundary conditions 
Sideslip and roll provide no great difficulties in this 
respect. A sideslipped flat plate at an angle of attack 
can be represented by a uniform distribution of the 
downwash perturbation velocity over the plan form, 
The rolling wing is readily represented by a linear 
spanwise variation of the downwash perturbation ve. 
locity. The yawing wing, however, is not so easily 
prescribed in terms of the boundary conditions; in 
fact, the yawing wing is not truly adaptable to a 
steady-state solution through a substituted warping or 
some such artifice and should be considered as an 
unsteady-flow problem. At the present time, deriva- 
tives based on an unsteady flow analysis for the yawing 
wing are not available. Results for these derivatives 
calculated through the use of strip theory and other ap- 
proximations, however, are included in the results pre- 
sented. 


In view of the boundary-condition situation and 
taking into consideration the nonlinearities that are 
involved in the sideslip analysis and the complication 
of the yawing-moment derivatives due to the edge- 
suction force, it is not surprising that the preponderance 
of the reported calculations of lateral-stability deriva- 
tives are for the damping-in-roll derivative, C),.  Be- 
yond this one derivative, however, the bulk of the ma- 
terial generally available on lateral-stability derivatives 
has come from a group, headed by Dr. H. S. Ribner, 
in the Stability-Analysis Section at the Langley Aero- 
nautical Laboratory of the N.A.C.A. and a group in 
the Theoretical Aerodynamics Section of the Ames 
Aeronautical Laboratory. Most of the results on 
supersonic stability derivatives that will be presented 
and discussed in this paper have come from one of these 
two sources. 


The values given for the derivatives are referred to the 
stability axes system.” This is a right-handed system 
of three orthogonal axes—namely, the longitudinal 
axis, the vertical axis, and the lateral axis of the air- 
craft. The axes are fixed in the airplane and travel with 
it. They are oriented in a manner such that the lon- 
gitudinal axis is in the plane of symmetry, positive 
forward, and aligned with the free-stream vector at 
zero sideslip. The vertical axis is positive downward 
and the lateral axis is positive in the direction toward 
the right wing tip. The signs for the moments follow 
the right-hand screw rule. 

The angle of attack a is measured in radians, and the 
sideslip angle 8 in degrees, wherever they occur herein. 
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LATERAL-STABILITY DE 


RESULTS AND DISCUSSION 


The presentation and discussion of the results will be 
divided into four categories: (1) the variation of the 
derivatives with aspect ratio; (2) the variation of the 
derivatives with Mach Number for both subsonic and 
supersonic speeds; (3) the variation of C, and C, with 
sideslip; and (4) the effects of wing-and-body com- 


binations. 


Four types of plan forms will be considered in this 
study. These plan forms, shown in Fig. 1 for an aspect 
ratio of 4, are characteristic of the supersonic wing con- 
figurations that are frequently proposed. They fit the 
following descriptions: (1) rectangular; (2) trape- 
zidal with tips raked out at a slope of 0.5; (3) triangu- 
lar with apex forward; and (4) swept back and tapered 
to points at the wing tips, having subsonic leading edges, 
supersonic trailing edges, and a ratio of root chord to 
overall length of 0.5. 


The sweptback plan form just described was selected 
because more data were available for it than for any 
other of the sweptback type. It is unfortunate that 
this plan form has not been analyzed as yet for condi- 
tions that make the leading edge supersonic or the 
trailing edge subsonic. With regard to this situation 
it can be said, however, that the sweptback plan forms 
as a Class present a rather complex integration prob- 
lem and, consequently, require rather long and tedious 
calculation procedures. All this accounts for the 
rather limited extent of the curves presented for the 
sweptback plan form. 


In the following discussion, the stability derivative 
coefficients—that is, the nondimensional form of the 
derivatives that are plotted on the graphs—usually will 
be referred to as the ‘derivatives’ for simplicity and 
brevity. 


The Variation of the Derivatives with Aspect Ratio 


The variations with aspect ratio of the six deriva- 
tives under consideration are shown in Figs. 2 through 
7. These results are presented for B = 1 (B is the 


Mach Number parameter ~/./? — 1) and were ob- 


Fic. 1. Plan-form types investigated. 4A = 4. 
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tained from the data and expressions contained in refer- 
ences 3 through 9. 

The Damping in Roll, C,,.—The negative values for 
C,, evident in Fig. 2 indicate damping of the rolling 
motion which is to be expected for unstalled wings. 
As might also have been anticipated, the derivatives, in 
general, increase in magnitude with increasing aspect 
ratio. The one exception is the supersonic-edged tri- 
angular plan form for which C,, remains constant at 
—'/3. The trapezoidal and rectangular plan forms ap- 
proach asymptotically the two-dimensional value of 
2 


— /'3. 

For the ratio of c,/l of 0.5, the aspect ratio of the 
sweptback wings is given by 8m, where m is the slope 
of the leading edge, as compared with 4m for the tri- 
angular wing. Thus it can be determined from Fig. 2 
that, for equal leading-edge slopes, the sweptback wing, 
which has twice the aspect ratio of the triangular wing, 
provides the greater damping in roll. 

The Yawing Moment Due to Rolling, C,,/a.—The 
leading-edge suction forces, as well as the normal forces, 
contribute to the value of C;,/a. The normal-force 
contribution is, for practical purposes, —aC), and is 
independent of the longitudinal location of the moment 
center. The contribution of the edge-suction force, 
however, is quite dependent on the moment-center lo- 
cation. Thus, the magnitude and, in some cases, the 
sign of this derivative can be altered by a chordwise 
shift of the moment center. These statements apply 
also to the other derivatives that contribute yawing 
moments, C,, and C,,. For the particular plan forms 
under consideration, the values presented for the yaw- 
ing-moment derivatives correspond to the following 
moment-center locations: c,/2 for the rectangular and 
trapezoidal plan forms; 2c,/3 aft of the leading-edge 
apex for the triangular plan form; and 2//3 aft of the 
leading-edge apex for the sweptback plan forms. Sub- 
sequent calculations have shown, however, that, for 
any location of the moment center between those indi- 
cated and the leading edge, the results presented may 
be considered as fairly typical. 

It is evident from Fig. 3 that the variation of C,,/a 
with aspect ratio is similar to the C,, variation. The 
values of the derivative generally increase with increas- 
ing aspect ratio, only in this case it is in algebraic rather 
than absolute magnitude because of the shift in the 
horizontal axis and the inversion of the vertical axis. 

The Rolling Moment Due to Vawing, C,,/a.—The re- 
sults for this derivative are presented in Fig. 4. The 
analysis presented in reference 6 is based on a strip 
theory and, for B = 1, predicts zero as the value for 
C,,/a for the rectangular wing. This analysis extended 
to any closely related trapezoidal plan form yields the 
same result. 

The triangular and sweptback plan forms with sub- 
sonic leading edges show, in general, an increase in 
Ci,/a with an increase in aspect ratio. Furthermore 
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ASPECT RATIO, A 
Fic. 2. Variation of damping-in-roll derivative with aspect 
ratio. B = 1. 
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Variation of yawing-moment-due-to-roll derivative with 
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in this case, as in the C;, case, equal leading-edge slopes 
generally yield greater magnitudes of C,,/a for the 
sweptback plan forms than for the triangular plan 
forms. In references 4 and 7, the sources of the re- 
sults for triangular and sweptback plan forms, the 
belief is expressed that, if the analysis were to be carried 
out for supersonic leading edges, a reversal in the sign 
C, would occur. It is also stated, with regard to the 
approximations involved in the analysis, that, although 
the spanwise increase in speed is taken into account, the 
spanwise variation in Mach Number is not included. 

The Damping in Yaw, C,,/a?.—The values of C,,/a? 
shown in Fig. 5 are for a potential flow. No effects of 
profile drag are included. Thus, the major contribu- 
tions to this derivative are again the normal-force 
moment, represented by —a(C;,,, and the moment due 
to the edge-suction force. 

As might have been anticipated, the variations of 

,,/ fit the descriptions given for C,/a for each of 
the plan forms. 

In general, the value of C,, contributed by the wing 
will be dependent on the effects of viscosity. Further- 
more, for an aircraft having fuselage and tail compo- 
nents, the contribution of the wing alone to the value of 
C,, for the complete configuration is likely to be 
small. 

The Dihedral Effect, Cig/a.—There are two curves 
shown in Fig. 6 for the rectangular plan forms. The 
dotted curve from reference 6 represents a result based 
on the hypothesis that the Kutta condition does not 
hold at the trailing tip for small angles of sideslip. The 
solid curve from reference 5, on the other hand, repre- 
sents a value of the derivative based on the assumption 
that the Kutta condition holds throughout the side- 
slip range. Experimental results will be needed to 
determine eventually just what condition or changes 
in condition are prevalent at the trailing tip during side- 
slip. At the Mach Number considered, J = 1.41 or 
B = 1, the two curves differ in sign throughout most 
of the aspect-ratio range. Both curves, however, show 
a general decrease in magnitude as the aspect ratio in- 
creases. 

The trapezoidal plan form selected for this investiga- 
tion provides positive dihedral effect throughout the as- 
pect-ratio range at this Mach Number. As shown in 
reference 5, however, this result cannot be generalized 
to apply to all trapezoidal plan forms of this type. 
For a different tip slope, the sign of the derivative can 
change within the aspect-ratio range shown. 

The results for the triangular plan form indicate a 
reversal from positive to negative dihedral effect as the 
leading edge becomes supersonic. This is essentially a 
Mach Number, rather than an aspect ratio, effect, and 
the primary result of increasing the aspect ratio is a 
reduction in the magnitude of the derivative. 

Although only a limited aspect-ratio range is shown 
for the sweptback wing, the results indicate that this 
plan form is behaving similarly to the subsonic-edged 


triangular plan form and, thus, probably will continue 
to do so when the leading edges becomes supersonic. 
Curves for both subsonic-leading-edge and supersonic- 
leading-edge sweptback plan forms of different types, 
presented in reference 5, indicate the probability of this 
behavior. 

The Directional Stability, Cyg/a?.—The differences 
between the magnitudes of the plotted values of C,,/a? 
and Cyg/a are due to the yawing moments contributed 
by the leading-edge-suction forces. The suction forces 
generally tend to increase the directional stability over 
that contributed by the normal-force component 
—aCg alone. The variations with aspect ratio for 
Cyg/a® are much the same as those for Cis/a, with the 
exception that the suction effects cause the derivatives 
for the subsonic-edged triangular and sweptback plan 
forms to increase rather than to decrease in magnitude 
as the aspect ratio is increased. 

The effects of assuming that the Kutta condition does 
not hold on the trailing tip are shown again by the 
dotted curve for the rectangular plan form.’ Further 
details on the sideslip variations of C,,/a? and C,/a and 
the effect of the Kutta condition are given in a later 
section of this paper. 


Variation with Mach Number 


The four plan forms selected for the investigation of 
Mach Number effects are shown in Fig. 1. By choos- 
ing the aspect ratio of these plan forms to be 4, the lead- 
ing-edge slopes of the triangular and sweptback wings 
were fixed at 1.0 and 0.5, respectively. 

The variations of the stability derivatives with Mach 
Number have been extended to cover the subsonic, as 
well as the supersonic, velocity régime in order to pro- 
vide a more interesting and significant comparison 
The values of the derivative for J = 0 were obtained 
from references 9 through 14 for all but the trapezoidal 
plan form. The results for the rectangular plan form, 
however, will be approximately correct for the trape- 
zoidal plan form at subsonic speeds. The references 
just cited contain some experimental, as well as calcu- 
lated, results, and some of the experimental values are 
included in the subsonic portions of the curves of Figs. 
8 through 13. With the exception of C,,, the values for 
the derivatives at m = 0 were approximately corrected 
for the effects of compressibility, and thereby extended 
to M & 0.95, by the methods of reference 15. 

The Damping in Roll, C,,.—The curves for the C, 
variation with Mach Number, shown in Fig. 8, follow a 
similar pattern for each plan form. The general trend 
is a slight increase in the magnitude of C), as the Mach 
Number increases and approaches unity, followed by 
some transition through the sonic speed to a comparable 
value of C,, and a fairly rapid decrease in magnitude at 
the low supersonic Mach Numbers. The curves then 
tend to level off, however, and begin a rather asymptotic 
approach to the axis as the Mach Number is further 
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Variation of rolling-moment-due-to-sideslip~ deriva- 
tive with Mach Number. A = 4, 
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Variation of yawing-moment-due-to-sideslip deriva- 
tive with Mach Number. A = 4. 
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The Yawing Moment Due to Rolling, C,,/a.—The 
negative sign of this derivative, which is associated 
with these plan forms at subsonic speeds, changes to 
a positive one at various supersonic Mach Numbers 
jor all the plan forms. For the rectangular plan forms 
this reversal occurs at the sonic speed, JJ = 1. For 
the trapezoidal, triangular, and sweptback plan forms 
it occurs at supersonic speeds. 

With regard to the magnitudes of the derivatives, it 
is evident that, in general, a decrease in the magnitude 
of Cy,/a@ takes place as the supersonic Mach Number 
increases. It is likely that, once the leading edges of 
the sweptback plan form become supersonic, this plan 
form would follow the trend of the others. 

The Rolling Moment Due to Yawing, C:,/a.—The 
curves for the rectangular and trapezoidal plan forms, 
shown in Fig. 10, indicate that reversals in the sign of 
this derivative occur at the sonic speed and exist up to 
M= 141. There are no changes in the sign of C,,/a 
for the triangular and sweptback wings, however. As 
the Mach Number is increased, the magnitudes of the 
derivatives for all the plan forms increase as the sonic 
speed is approached and, in general, decrease at the 
supersonic speeds. 

The Damping in Yaw, C,,/a?.—The reversals of the 
sign of C,,/a* for the rectangular and trapezoidal plan 
forms which occur when the sonic speed is exceeded 
hold only until the Mach Number 1.41 is reached. 
The value for the triangular wing shows a sizable jump 
in magnitude at the sonic speed but no reversal in 
sign; whereas the value for the sweptback wing shows 
both. The lack of damping in yaw for the sweptback 
wing at subsonic speeds occurs only at large angles of 
sweep. This same condition in the cases of the rec- 
tangular and trapezoidal wings at the low supersonic 
speeds, however, is due to the Mach Number effects on 
the spanwise distribution of the loading.* 7 

As mentioned previously, the values of C, shown in 
Fig. 11 represent only the potential flow contribution 
of the wing alone. The other aircraft components, such 
as the fuselage and the vertical tail, and the viscous 
effects will generally contribute a large part of the mag- 
nitude of this derivative. 

The Dihedral Effect, Cig/a.--The curve for the tri- 
angular plan form, shown in Fig. 12, indicates that this 
plan form loses its positive dihedral effect when its 
leading edges become supersonic. The resulting nega- 
tive dihedral effect diminshes rapidly, however, as the 
Mach Number is further increased. The data avail- 
able in reference 5 on sweptback plan forms indicate 
that the sweptback wing considered would probably be- 
have in a similar fashion. 

For the rectangular plan form, the switch from posi- 
tive to negative dihedral effect occurs at J = 1. No 
curve has been shown for the case where the Kutta con- 
dition on the trailing tip is neglected, but it should suffice 
to say that this curve would indicate positive dihedral 
effect throughout most of the supersonic Mach Number 
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range. At the particular aspect ratio of 4, the trape- 
zoidal plan form with slightly raked-out tips maintains 
a positive dihedral effect throughout the Mach Num- 
ber range. As indicated previously in the aspect-ratio 
study, that result cannot be generally extended to all 
trapezoids of this type. The aspect ratio and the 
amount of tip rake are rather critical factors with re- 
gard to the value of this derivative at supersonic speeds. 

On the supersonic side of the Mach Number range, 
the magnitudes of Cig/a generally tend to decrease with 
increasing Mach Number for all plan forms. 

The Directional Stability, Cyg/a?.—The curves for 
Cyg/a® shown in Fig. 13 follow the same general trends 
as the curves for Cig/a. As explained previously in the 
discussion of the variation of these derivatives with 
aspect ratto, any differences in the curves for C,,/a?* 
and Cjz/a@ at supersonic speeds are due to the effect of 
leading-edge suction on C,,/ a’. 


Variation of C, and C, with Sideslip 


The nonlinearities in the sideslip variations, which 
are illustrated by the C,,/a? and C,/a curves for triangu- 
lar and rectangular plan forms in Figs. 14 and 15, are 
of special interest. The solid curves, taken from refer- 
ences 5 and §, illustrate the three types of nonlinearities 
commonly occurring: (1) the warped but still continu- 
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ous curve present in both the rectangular and triangular 
cases; (2) the break in the curve for the triangular 
wing; and (3) the sudden jump in the curve at zero 
sideslip for the rectangular wing. The first type is due 
mainly to the generally nonlinear edge-suction effects. 
The second type occurs usually as an edge passes 
through a Mach cone and changes from a subsonic 
edge with suction to a supersonic edge without suction, 
or vice versa. The third type is usually caused by an 
edge passing through the free-stream direction and 
changing from a leading edge with suction to a trailing 
edge without, or vice versa. This third type of non- 
linearity occurs if the Kutta condition is assumed to 
exist at the trailing tip. The dotted curve for the rec- 
tangular plan form in Fig. 14, obtained from reference 
6, represents the assumption that the Kutta condition 
does not hold. It can be seen that, in addition to 
eliminating the jump in the curve, the latter assumed 
condition changes the sign of the slope of the curve. 
This change in the sign of the derivative is also evident 
in the C, variations plotted in Fig. 15. 


Effects of Wing-and-Body Combinations 


It is certain that the values of the stability deriva- 
tives for the wing alone will be modified somewhat by 
the addition of a body such as a fuselage. Unfortu- 
nately, not much material is available as yet on the 
lateral-stability derivatives for wing-and-body combi- 
nations. Some idea, however, of the effects on Cig 
and C;, of adding a fuselage can be obtained for low- 
aspect-ratio plan forms from references 16 and 17, 
respectively. These results are based on a theory de- 
rived particularly for low-aspect-ratio wings and slender 
bodies. This theory is applicable throughout the 
whole Mach Number range, including the sonic speed. 
Fig. 16 contains the ratios of these derivatives for the 
wing and body combined to those for the wing alone, 
()/()w, as a function of the ratio of the diameter of the 
body to the wing span, D/b. A triangular wing and a 
circular cylindrical body are represented. The entire 
wing area, submerged as well as exposed, is the refer- 
ence area on which the coefficients are based. As might 
be expected, the ratio Cig/(Cig). diminishes steadily as 
the ratio D/b increases. The rate of decrease of C);/ 
(Cig)w also increases steadily. A rather remarkable 
result was obtained for the C,, ratio, however. For 
values of D/b less than 0.4, the ratio of Cip/(Cy)w is 
greater than 1. This indicates that the interference 
due to the body either increases or favorably redistrib- 
utes the wing loading for this range of D/b. At higher 
values of D/b, however, the C,, ratio tends toward 
zero as might be reasonably anticipated. 


CONCLUDING REMARKS 


From the foregoing collection of material on the 
lateral stability derivatives at supersonic speeds it is 
possible to pick out a few definite trends. The magni- 
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Fic. 16. Stability-derivative ratios for the damping in roll and 
the dihedral effect of a wing-and-body combination. 


tudes of the supersonic derivatives in coefficient form, 
for instance, are fairly comparable to those encountered 
at subsonic speeds. The effect of increasing the as- 
pect ratio is, in general, an increase in the magnitudes 
of the stability derivatives due to rolling and yawing and 
a decrease in the magnitudes of the derivatives due to 
sideslip. 

With regard to the effects of Mach Number, the 
magnitudes of the derivative coefficients generally tend 
to increase as the sonic speed is approached from the 
subsonic side and to decrease as the Mach Number 
is increased on the supersonic side. 

One of the most interesting aspects of the behavior of 
the results was the reversal in sign for some of the de- 
rivatives. These reversals usually occurred either at 
the sonic speed or at supersonic speeds as the subsonic 
leading edge changed to a supersonic leading edge. Of 
the derivatives contributed by the wing alone, those 
contributing rolling moments, C,,, C;,, and Cig are the 
most significant in their effect on the stability of the 
aircraft. Two of these, C, and Cig, show reversals in 
sign from their usual subsonic values. The change in 
having the right wing drop instead of pick up when the 
right wing speeds up in yaw or when sideslipping to the 
right will probably pose serious stability problems at 
supersonic flight speeds. Further analysis is needed, 
however, in order to bring out the complete physical 
significance of all the trends and magnitudes of the 
derivatives presented. 
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Influence of Acceleration on Aerodynamic 
Characteristics of Thin Airfoils in 


Supersonic and Transonic Flight’ 


integral expressions, 
Number and on an acceleration parameter. 

The evaluation of these solutions is carried out in a different 
manner, depending on whether the transonic or the supersonic 
range is considered. 
speeds is expanded in positive powers of the acceleration parameter 
3. The influence of acceleration turns out to be extremely small. 
In Section (B), the transonic range is treated, by expanding the 
expressions for lift and drag, at M = 
the negative one-fourth power of 8. This elucidates why, at M = 
|, the theoretical expressions for the aerodynamic parameters re- 
main finite in the case of accelerated motion but become infinite 
in the limit of stationary flight speed (8 = 0). 
a physical discussion of the results for wedges and diamond pro- 
The curves representing the aerodynamic param- 
eters as a function of the Mach Number have a more or less 
pronounced maximum at JJ = 1, depending on the value of the 
flight acceleration 8. The application of a validity criterion 
shows that for sufficiently 
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ABSTRACT 


An attempt has been made to derive the aerodynamic charac- 
teristics, in the supersonic and transonic range, of a two-dimen- 
sional airfoil moving at uniformly accelerated speed; it is assumed 
that the airfoil is thin and the shocks resulting from its motion are 
extremely weak. Under these conditions, the airfoil motion may 
be treated as a time-dependent perturbation problem, which can 
be solved by means of ‘‘Possio solutions.’’ The solutions for lift, 
drag, and pressure on the airfoil surface turn out to be elliptic 


whose arguments depend on the Mach 


In Section (A), the solution at supersonic 


1, into a series starting with 


In Section (C), 


“large”’ 6 the results can be taken quan- 


“small”? 8 the results have only qualitative 


NOTATION 


coordinates of point P in reference system 
fixed in still air 

coordinate in reference system, fixed in mov- 
ing airfoil 

strength of ‘‘temporary source’ in 
solution 

coordinates of 
solution 

sound speed in still air 

static pressure, density in still air 

instantaneous Mach Number 

chord length of airfoil 


Possio 


“temporary source’’ in Possio 


Ludloff at the Aero- 


Daniel Guggenheim School of Aero 


a = angle of attack of chord 

6 = vertex angle of wedge profile, semivertex angle 
of diamond profile 

b = acceleration 

8 = —2b//a? = dimensionless acceleration parameter 


INTRODUCTION 


5 pew PATH OF A uniformly accelerated, two-dimen- 
sional thin airfoil can be represented in the x-¢ 
plane as a strip between two parabolic arcs (see Fig. 1), 
the width of the strip being /, where / is the chord length. 
If the airfoil is assumed to produce small disturbances, 
any point inside the strip may be considered to be a 
source of disturbance waves propagating with sound 
speed. The velocity field can be derived from a per- 
turbation potential that has to fulfill the differential 
equation of ‘‘sound’’ waves, and the solution of this 
equation may be represented as an integral over ‘““Tem- 
porary Sources’ (Possio’s Solutions'), which are dis- 
tributed over the parabolic strip mentioned before: 


g(x, t) 


(é, 


At a given station P on the airfoil surface, only those 
disturbances will be felt which are within the shaded 
area of Fig. 1, that is limited by the “sound speed lines, 
passing through P. Thus the domain of integration has 
to be extended over this shaded area 

The characteristic difference between the supersonic 
and subsonic case is indicated in Fig. 1. While in super- 
sonic motion the pressure at P is influenced only by a 
small finite source area, whose vertex is at P and whose 
base is the parabolic path CD of the leading edge, in 
the subsonic case the sound speed-lines through P cut 
off the whole lower part of the parabolic, source-covered 
strip. This means that during the subsonic flight, 
the whole previous history of the airfoil motion will 
influence the pressure distribution, lift, and drag, while 
during the supersonic flight the airfoil flies through un- 
disturbed air. 

Thus it may be noted that the aerodynamic charac- 
teristics at JJ = 1 will be different, depending on 
whether the airfoil has been accelerating from subsonic 
speeds or decelerating from supersonic speeds. 
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The present investigation deals with the change of 
aerodynamic characteristics during a decelerating 
flight, starting at supersonic speeds (see author’s note,” 
December, 1948). It turns out that the influence of 
deceleration is hardly appreciable at supersonic speeds 
(slightly more appreciable, if the airfoil flies at a finite 
angle of attack). But in the transonic range (see Fig. 3), 
the shape of the pressure-, drag-, and lift-curves (i.e., 
Cp vs. M, etc.), as derived from this theory, depends 
strongly on the value of the acceleration parameter £. 
While the theory becomes invalid in the limit of station- 
ary flight speed (i.e., for 8 = 0), its degree of validity 
increases rapidly with increasing 8-value [see Section 
(C, 3)]. Already at 8 ~ 0.02, the theory is as good an 
approximation as the linearized theory in the supersonic 
range. For 8 ¥ 0 all curves have a marked maximum 
at M = 1, but all parameters remain finite and vary con- 
tinuously through the transonic range. 


The finiteness of cp, etc., at JJ = 1 for accelerated 
motion, as in contrast to the infinite cp-value attained 
for stationary speed, is made plausible in Fig. 2. It is 
seen that the parabolic strip of Fig. 1 becomes a straight 
strip for stationary speed; consequently the shaded 
integration area becomes of infinite extent for 8 = 0— 
i.e., the integral expressions for Cp, cz, and c, diverge 
for constant speed. 


In more precise terms: At supersonic speeds, the aero- 
dynamic quantities may be expanded in positive powers 
of B,solongas —1> V/B. However, in the vicinity 
of M = 1 (i.e., for WV — 1 < V8) the elliptic integral, 
representing the pressure, drag, or lift value, has to be 
expanded in a different manner—i.e., by a series start- 
ing with a negative one-fourth power of 8. This eluci- 
dates the fact that in the limiting case of constant 
flight speed (8 = 0), the well-known divergent values 
for lift and drag are obtained, while for uniformly 
accelerated speed, all aerodynamic parameters, as lift, 
drag, and pressure on the airfoil surface, remain finite 
at = 1. 
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Path of motion of decelerating airfoil in x-t plane. 


A singular behavior is evidenced by the theoretical 
pressure value right at the leading edge, which can be 
shown to be finite at supersonic speeds, but becomes 
infinite as x~' at sound speed, and logarithmically 
infinite at subsonic speeds. The latter can be antici- 
pated, in analogy to the stationary thin airfoil theory, 
Of course, this singularity of the pressure at the leading 
edge leaves the lift and drag values finite. 


At the conclusion of these considerations, it has been 
checked, to what extent the results of the linearized 
theory around J = 1 may be expected to be in fair 
agreement with actual conditions. Therefore a validity 
criterion is discussed in Section (C, 3). Its application 
shows that under certain conditions—i.e., for small 
thickness ratio [and small A.R.; see Section (C, 3)]— 
the results derived here should be valid approximations, 
while in other cases the conclusions drawn for the tran- 
sonic range are merely of qualitative nature. In this 
connection, one may also consult a recent paper by Lin, 
Reissner, and Tsien,* in which it is explained why line- 
arized treatment of slender bodies in the transonic range 
becomes legitimate for sufficiently large accelerations. 
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(A) COMPUTATION OF DRAG AND LIFT OF A 
DECELERATING THIN AIRFOIL IN SUPERSONIC FLIGHT 


(1) Formulation of the Analytic Problem 


Suppose a thin airfoil moves with uniformly accel- 
erated, supersonic speed U + Ot into still air, in the 
direction of the negative x-axis. In coordinates, fixed 
in the airfoil, the upper profile surface may be given by: 
y = f(@); here f(0) = f(/) = 0,/ being the chord length 
isee Fig. 4). In coordinates fixed in still air, y = f[x + 
Ut + (1/2)bt?] where x, y, and ¢ are so chosen that, at 
t = 0, the leading edge is at the origin x = y = 0. 
As it is assumed that the airfoil is thin and flat—i.e., f 
<< land f’ << 1—the disturbance, caused by the 
motion, is a small perturbation of the state of rest, and 
can be described by a perturbation potential g which 
obeys the differential equation of disturbance waves: 


Ay — (1/a*)gu = 0 (1) 
a being the sound velocity in still air. 
The pressure increment due to the disturbance is 
given by 
Ap = —p(¢,) (2) 
where p is the density of still air. 
To determine the solution of Eq. (1), initial and 


boundary conditions must be specified. Assuming that 
the motion started at some past time ¢ = —7, the ini- 


tial conditions are 
¢g¢=0, =0 att = (3) 


The boundary condition on the upper surface of the 
airfoil is in linearized approximation : 


gy = (U + + Ut + (1/2)bt?] 


t+ Utt (1/2)b2 <1 


forall? > —T. 

Because of the symmetry of the flow pattern about 
y = 0, it follows that for y = 0 at a point of the x-/ 
plane exterior to the airfoil, g, must be zero. Thus if 
the definition of f is extended so that f(#) = 0 for # < 
> /, it may be assumed that 


gy = (U + + Ut + (1/2)bt*] (4) 


holds at y = O for all values of x and ail values of t > 
-T. This extension of the definition of f(x) will be 
useful in the course of the derivation [see Section 


(A, 2)]. 


Co vs M- WEDGE PROFILE 
Cu vs M— FLAT PLATE PROFILE, DIAMOND PROFILE 


Fic. 3. cp vs. M and cy vs. M for wedge profile, at several ,- 
values. 


U + bt 


AIRFOIL VELOCITY 


K= x + Ut + got? 


Fic. 4. Cross section of airfoil in x-y plane. 


| 


Fic. 5. Composition of diamond by wedges. 


The solution of the initial-boundary problem defined by Eqs. (1), (3), and (4) can be given as the potential of a 
layer of temporary “‘sources”’ distributed over the x-t plane, the potential of a “source’’ of strength o, located at 


(&, 0, 7) being 


7) Vit — 7)? — (i/a)[(x — + y?] 
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where the domain of integration is the region where —7’ < + < ¢ and the denominator is real. [Any point insid 
the so-defined -7 domain may contribute to a disturbance felt at the point (x, y, f).] | It has been shown by many 
authors'® that, if yg is represented in this way, then 


Thus 


g(x, y, t) 


lim ¢,(x, y, = —7o(x, t) (6 
yoot+ 


Thus substituting from Eqs. (4) and (6) into Eq. (5) and, since the airfoil surface is of primary interest, setting 


y = 0: 
fa fa (U + br)f'lé + Ur + 
x —a(t—r) V(t — 7)? — (1/a?) — 


It is convenient to regard ¢ as a third coordinate, so that the projection of the airfoil motion on to y = 0 traces out 
a parabolic strip in the x-¢ plane (see Fig. 1). If P is the point (x, 0, f), it can be seen that the potential at P js 
given by an integral over the domain PAB cut out by the line ¢ = —7 and the lines PA, PB along which é = » 
* a(t — r)—.e., the lines representing the motion of a disturbance starting at P(x, 0) and propagating with 
sound speed (“‘sound speed lines’). It may be seen that after a sufficient time has elapsed, the lines PA, PB 
both intersect the path of the leading edge before they intersect / = —7. It will be supposed that this is the case 
for all points P which are being considered. (Otherwise—i.e., if the line ¢ = —T is part of the boundary of the 
integration area—a transient effect connected with the start of the motion enters the problem.) The integrand 
then vanishes except over the shaded region PCD, since f’ vanishes ahead of the leading edge. In particular, the 


integrand vanishes, as can be seen from Fig. 1, in a certain r-interval say between 7 = —7 and +r = —T +1. 
Thus, upon introducing a new integration variable 7 = ¢ — 7, Eq. (7) can be written as 
b¢ — + UE — 7) + 2)b(t — 7)*] 
0, t) = = (8) 
— GF = t), 


According to Eq. (2), the pressure at the surface of the airfoil is proportional to ¢,. To find g, differentiate Eq. 
(8) and set afterward ¢ = 0, which is a convenient transformation of the time coordinate, but no restriction of the 
generality, as the time zero point may be chosen at convenience: 


star (U — b7)2f"[— — Ut + (1/2)b — UF + (1/2)b7" 


(2) Expansion in Powers of a Deceleration Parameter 8 


Since, in general, the acceleration of an airfoil is small 7+ 
compared to a standard acceleration a?/l, a series ex- is 
pansion in terms of an acceleration parameter should be 6+ a . ar 
attempted. To obtain an expansion of ¢,(x, 0, 0) in atl cncnail 
powers of b, the integrand in Eq. (9) may be expanded 5+ ie 
in a Taylor series and the series integrated term by &e ae 
term. Of course, this is possible only if f is an analytic 4t 
function of its argument. However, this is in reality no 
restriction, as any f can be approximated uniformly 3+ 
throughout the domain of integration as exactly as £ 
desired, by an analytic function. The extended defini- 2 
tion of f(x) as introduced in Eq. (4) makes it possible 


to keep the integration limits in Eq. (9) free of b, so 00s 6~«CkS:*é<CSS:S:C SC 
that only the integrand has to be expanded ;. this simpli- B= 


fies the procedure greatly. Fic. 6. cp vs. 8, for M = 1 and M = 1.3, of wedge profile. 
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Thus, expanding the numerator in Eq. (9): 
— — Ut + (1/2)b7?] + bf’ [— — UF + (1/2)b7?] 


n = ofl. n = 0 (10) 
| 


2n(2n — 1)?" ~ — UF)] 


where | as denotes the kth derivative of f. This expansion (10) must now be considered substituted into the 


numerator of Eq. (9). 
In order to carry out the termwise integration, required by Eq. (9), the following procedure is used. If g(x) 


isan arbitrary function, let L(g) be defined by 


Ug'(t — UF 
s-@ — [(x — é) a]? 


Itcan be seen that by setting b and ¢ equal to zero in Eq. (8), L(g) is the potential on the surface of an airfoil, whose 


profile is given by y = g(x) and which moves at constant speed U in the direction of the negative x-axis. Thus, 


by the Ackert theory, 
Lig) = —Ug(x)/V M2 — 1, where M = U/a 


Now, it is to be noted that the kth integrated term of series (10) can be obtained by differentiating Eq. (11) & times 
with respect to U. Thus 


_ o* 
k k — (12) 


I | k! |p ( M ) (x) (12a) 
—| ——— —)o(x 2: 
y a*(M? — 92)" * VM? — 1 . 


where P,(x) is the Legendre Polynomial of degree k. In particular, if k = 2” and g(x) = f‘~ "+" (x), mean- 
ing the n-fold repeated indefinite integral of f’(x), Eq. (12) yields: 


x + a7 + 2) (é UF) (2n)! 
Ton fas at [(x £) al? a?"(M? 1)° 2n + 1)/2 * 2 f (x) 


or, since * = [(x — '/(n — 1)! (Ode, 


Ton F | Me 1)" + 2(n (x t) f' (dt (13) 


As shown in Appendix (1) 


Similarly, 


(2n — 1)! P (- (x t)" f'(t)dt (14) 
9 = = 9 ) 2n — 


and 


(2n — 2)! ( M 
In - = x — t)” "(t)dt 15 
< a” 2(M? D/2(4 1)! 2 2 VANE 0 (x ) f ( ( >) 


In the termwise integration of the series expansion to be substituted in Eq. (9), Eqs. (13), (14), and (15) may be 


used and in addition the relation 


(n + 1)P,, 4 — (2n + 1)sP,(z) + nP, = 
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may be taken into consideration. Thus 


M? — 4n)M 
= = + 2n) Poy - (16 


From Eq. (16), the pressure coefficient of a wedge profile of vertex angle 6 (i.e., f(x) = 6x) whose lower side 
parallel to the flight direction, can be obtained as: 


where the acceleration parameter 8 = —2bx/a?. Then the drag coefficient 
1 2 M?+2 1 M* + (27/4) M2 + 1 


For symmetrical airfoils of chord length /, the drag coefficient is given by means of Eq. (2) as: 


Cp = 2Ap(x)f'(x)dx = — af vif’ (x)dx (19 


M? — 2n)Po, t) 


where now 8 = —26//a*. Here it is assumed that the profile function f(«) is subtended above the flight direction; 


therefore f(/) % 0, depending on whether the angle of attack is zero or not. 


Indicating the first terms more explicitly, Eq. (20) yields for zero angle of attack: 


(24.7% + 348M4 + 303M? + 18) 
0 


Here, according to definition, f(/) = 0, so that the term proportional to the first power of 8 vanishes in this case 


In case of an angle of attack, f(x) is resolved geometrically into a profile function of f(x), subtended above the 
chord, and the angle of attack a of the chord; thus for symmetric airfoils: 


2 


1 M* + (27/4) M? ] 
4 M?(M? — 1 at - + .. (2 


The lift coefficient of a symmetrical airfoil can be obtained by use of the relation 


here = 


2/1) So cpdx = — (4/1U) So! dx (23 
Thus at an angle of attack a: 
4 8 Mt+2 M4 + (27/4 
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M 

Fic. 7. ¢p vs. WV for diamond profile, at 8 = 0.02 and 8B = 0.05, 
Eqs. (18), (21), and (24) have been used to plot Figs. 


3and 7, in which cp and c, vs. J are represented at 
given values of 8, for wedge ant diamond _ profiles. 
The form of the above series indicates that the influence 


of the second and higher terms becomes appreciable 


only if 8 ~ OC. — 1). 
0.05 [see end of Section (C, 2)]. 


Under actual conditions 8 < 


(B) COMPUTATION OF AERODYNAMIC CHARACTERISTICS 
IN THE TRANSONIC RANGE 


(1) Formulation of the Analytic Problem 


From Eqs. (18) to (24) it can be seen that, near JJ = 
1, the series for cp, Cp, and c, diverge, as each term be- 
comes infinite. It is shown in Appendix (2) that the 
series will converge, so long as 8 < (JJ — 1)”. Ifa 
representative value [see end of Section (C, 2)] of B— 
eg. 8 ~ 0.02—is taken, the series expansions fail for 
M < 1.14, which means that the whole transonic 


regime must be represented in another manner. To 
investigate the transonic range, it is necessary to evalu- 
ate integrals (S) and (9), without making use of the 
benefits which an expansion in 8 may involve. In 
order to simplify the evaluation, the investigation will 
be restricted to diamond-shaped airfoil profiles, which 
can be considered composed of two wedges (see Fig. 5). 

The present objective is, therefore, to compute the 
pressure distribution and the drag on a wedge profile 


[U + b(t — 


02 04 06 08 02 04 06 08 
B B 
Co vs B-WEDGE PROFILE Co vs 8B ~DIAMOND PROFILE 


Fic. 8. cp vs. 8, at M = 1, for wedge and diamond. 

of vertex angle 6 and length x. Under these conditions, 
it is no longer profitable to make use of the extended 
definition of f(x); but the actual integration limits 
have to be substituted in Eq. (8), explicitly. Thus, 
the integration in and + must be extended over the 
shaded area limited by the two ‘‘sound speed lines”’ 
going through P and the parabolic segment CD lying 
between them (see Fig. 1). The parabolic path of the 
leading edge is given by & = —Ur — (b/2)r’, and its 
intersections with the two sound speed lines & = 
x + a(t — r) may be denoted by &, 71, and &, 72, respec- 
tively. It is appropriate to divide the area of integra- 
tion into an upper triangular area and the rest of the 
domain, as indicated in Fig. 1. Thus, setting « — & = 
~ and f’ = 6, Eq. (8) can now be written: 


x +U(t — ('/2) b(t — 
dé X 


U+b(t—#) “ar U + W(t 3) 
at 

a? — £*/a* 

(25) 


In order to get the pressure distribution, relation (2) 
may again be used; ¢; is evaluated by differentiating 
Eq. (25) with regard to the upper and lower limits as 
well as under the integral sign; as a result of this, four 
line integrals and two surface integrals are obtained. 
Two of the line integrals cancel and a third one vanishes, 
so that finally: 


~- — dz X 
Vv — (1/a%) [x + Ut — + (/ — 


+ U(t — 7) + ('/2) b(t — 7)? dé — at dé 
di 


(26) 


While the third integral in Eq. (26) yields an elementary function, the first two can be reduced to combinations of 
complete and incomplete elliptic integrals of the first and second kinds. 
therefrom will be discussed later. 
with regard to x, as indicated in Eq. (18). 


The pressure distribution following 


Proceeding now to obtain the drag of the wedge profile, ¢, has to be integrated 
It can be shown that, putting ¢ = 0, as before: 


* 
| 17 
7 a 
6 = 
r | 
5 Mei Mei 
" 
14 
13 
3 
| j 
(18 
| 
(19 — 
} ; 
| 
| 
q 
is 
| 
(21) 
Case 
the 
| 
| 


54 JOURNAL OF THE AERONAUTICAL SCIENCES—JANUARY, 1950 


af (U — 67)? + dix — — Ur + (0/2) 


Differentiating Eq. (27) with regard to x again yields four line integrals and two surface integrals that after son 
calculation can be reduced to Eq. (26). This verifies Eq. (27), except for an integration constant which is dete. 
mined by the condition that the drag has to vanish for x = 0. 


(2) Evaluation of Solution (27) 


To evaluate Eq. (27), the numerator of the integrand is split into two parts, one part depending on 7, the othe 
on £ only. In this way f;’¢, dx is obtained as the sum of four double integrals, each of them containing £ as wel 
as 7 as integration variables. The integrations after — can be immediately evaluated and yield elementary fun. 
tions. Upon substitution, the sum of the four outer integrals (over 7) yields, after introduction of the dimension. 
less quantities: 


z= (b/a)r; 2, = (b/a)r, 4 = 1,2,3,4; B = —2bx/a?; M=U/a (28 
the following result: 
— dx = (J — L)(2a*x/p) (29 
here: 
[M26 — (1/2)8?] + [(3/2)MB]z + + (8/2) Mz* + (1/2)z! 
ac fru [— (1/2)67] + + 2[1 + (8, 2) — — (2M)z* — (1/2)z* 
VV (5 — 2)(z — — 23)(2 — 24) 
while = —(M — 1) + V(M— 1)? +8, —(M+1)+ V(M + 1)? + 8, = —(M —1)- 


V(M — 1)?+ Banda = 1) —- V(M41)?4+ 6; 21 >2>23>%. 

Both J and L are elliptic integrals of the same type, the numerator, say R(z), and the denominator, say P(:, 
being polynomials of fourth degree in z. The procedure by which either of the integrals is reduced consists of tw 
steps. In the first step, the degree of the numerators of either J or L is reduced from four to two, by observing 


that 
* R(z) d T(z) 
-dg = S(z) V P(z)\d === de 
V P(z)‘ ( (z) + V P(e) 


so that: 
R = (1/2)SP’ + S’P + T (31 


where 7 is an unknown polynomial of degree two and S is an unknown polynomial of degree one. Since R and? 
are given polynomials of fourth degree, Eq. (31) yields five equations for the three unknown coefficients of T and 
the unknown two coefficients of S. 

Thus, in Eq. (29), J — L can be replaced by It, — Me, where 


A,+ Hz+ C2? 
V (21 — 2)(z — 2) (z — 23)(2 — 2%) 
In the second step, IN; can be reduced by standard methods* to an expression of the form: 
Om, = RK(k) + S,E(k) + T,{(r/2) + [K(k) — E(k)] F(R’, 0) — K(R)E(R’, 0)} (32b) 


Here the K or F, respectively, and the £ are elliptic integrals of the first and second kinds. The R,, S,, and 7, 
may be expressed in terms of A,, H,, and C,, and the latter ones again in terms of the coefficients in Eq. (30). 
With these substitutions IN, and I, may be computed by means of Eq. (32b). Then the drag becomes, by use of 
Eqs. (29) and (32): 
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D' = Ap(x) f’(x)dx = 6 
0 0 


and 


Cp = 


V/opU2x 


For arbitrary / and 8, the quantities appearing in Eq. (32b) are given in Appendix (C). 
As for all practical purposes 8 < 1, the evaluation of Eq. (33) is considerably expedited, if R,, S;, Ti, and the 
dliptic integrals, which according to Eq. (32b) compose 9M, are expanded in powers of 6'*. 


At M = 1, the expansion yields: 
y, | 
2 12 
Si’ = B + 


1 l ’ 
| AL | 


19 
= 8 Att et 


Expanding the elliptic integrals at k = 1 ‘4/2 and @ = cot~! (1/+/2), one gets: 


K=K+ i(k — +...0(8) 


,1 


+ (K — E)F(k’, 0) — KE(k’, 0)] 


here 


(K 


...0(8) 


= — RK) (1/6)Kp'* + ...] 


K = K(1/V2) and E = E(1/V2) 


Inserting Eq. (34) into Eqs. (32) and (33), one obtains at J = 1: 


and 


since = —2bdx/a*. Since 


substitution of Eq. (35b) yields: 


and 


K 


Pung Po = a?2""*KB gil? 4. 


At any , the pressure coefficient can be obtained by evaluating Eq. (26). 


and introducing dimensionless variables yields: 


Cc. = 


(33) 


(34a) 


(84h) 


(35a) 


(35b) 


(35c) 


(36a) 


(36b) 


Carrying out the inner integrations 


(36c) 


Since Eq. (36c) has the same form as Eq. (32a), it can also be reduced to standard form [Eq. (32b)] and be eval- 


uated in the same manner. 
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(C) PuysicaL Discussion oF RESULTS 


(1) The Theoretical Results at and near M = 1 


From Eq. (36a) it can be seen that at JJ = 1, the 
pressure becomes infinite, at the leading edge of the 
wedge profile as x~'" for a given deceleration (—d). 
Also, at any given station on the wedge surface, x + 0, 
p increases as (—b)~'“*; for nearly constant speed (very 
small b), the representation becomes invalid. 


It might be guessed from analogy to the stationary 
case that for subsonic velocities the pressure at the lead- 
ing edge of a thin airfoil has a logarithmic infinity with 
x, and it can be seen from Eq. (16) that at supersonic 
velocities remains finite. These features of 


for MI = 1, can be derived from Eq. (26) after the inner 


integrals in Eq. (26) have been evaluated, and the 
outer integrals reduced to standard form. 


From Eq. (35a), it can be seen that at 1/ = 1, the 
drag of the wedge decreases with the length x of the 
wedge surface as x". From Eq. (18) it follows that, 
for supersonic velocities, the drag decreases, propor- 
tional to x. It can also be shown that for subsonic 
velocities, D ~ x (log x). 


From Fig. 1 it can be seen by inspection that for any 
point P between L.E. and T.E., 0 < x < J, the magni- 
tude of the shaded area, changes continuously, as P 
moves through the x-axis, where J = 1 (i.e., where the 
sound speed lines are parallel to the tangent of the 


vs M— FLAT PLATE 


cp vs. M for finite wings of different A.R.’s. 


parabola). This is confirmed by investigating the 
expression (26) for g;—i.e., p, in which integrand and 
integration limits change continuously, by passing 
through the sound speed. 


The singular point x = 0 is seen to be the only ex- 
ception, as in the denominator of the line integral in 
Eq. (26), the term proportional 7* vanishes for U =a 
andx = 0. Therefore, in the immediate vicinity of the 
leading edge of the wedge (or airfoil), the pressure dis 
tribution gradually changes its character, by passing 
through JJ = 1. This does not affect, however, the 
continuity of drag and lift, that as integrals over the 
p-distribution remain everywhere finite and change 
continuously by passing through the sound speed. But 
the slope of the cp or ¢, vs. \/-curve experiences a dis- 
continuity as shown in Figs. 3 and 7. 


Eqs. (35b) or (33) and (18) have been used to plot in 
Fig. 3, Cp vs. M for a wedge profile of vertex angle 4, 
at several @-values. Since, in supersonic flow, the 
pressure distribution on the upper side of a wedge is the 
same as on the upper (or, except for the sign, on the 
lower) side of a flat plate of the same inclination, Fig. 3 
represents at the same time the c, vs. /-curve of a flat 
plate profile, at angle of attack a = 4, if the ordinate 
scale is changed. The curves in Fig. 3 expose the above- 
described features for cp and c,; for comparison the 
or ¢, vs. \/-curve at constant speed (8 = 0) is added. 


Eqs. (36c) and (19) have been used to plot the pres 
sure distribution along the surface of the wedge profile 
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(a) for a supersonic speed: MW = 1.3, (b) for the sonic 
speed : M = 1, with 8 = —2bx/a* ranging from 0... . 
0,03 (see Fig. 6). It is interesting to note how the pres- 
sure increases, near 7 ~ 1, at all stations of the surface, 
but particularly in the vicinity of the leading edge. 
These results should be regarded merely as qualitative 
indications. 

Thus, these considerations have shown that the aero- 
dynamic characteristics at supersonic speeds (J — 1 > 
V2) are hardly affected by 6. But in the transonic 
range (MM — 1 < 1/8), the theoretical values of c¢,, 
(p, and ¢,, as computed by the linearized theory, de- 
pend strongly on the magnitude of the acceleration 
parameter. Since this theory is seen to become invalid 
in the limit of stationary motion at constant speed (6 = 
0), the results at small values of acceleration cannot be 
taken at face value. However, at high accelerations 
(8 ~ 0.02... 0.1) these approximations become in- 
creasingly more valid, as shown in Section (C, 3). 
Thus, the spread of the cp-curves, at various @, in the 
transonic range, may indicate that there the influence 
of acceleration is more essential than in the supersonic 


range. 
(2) The Drag of a Diamond Profile 


It remains to compute the drag of a diamond-shaped 
profile by adding the drag contributions originating 
from the two wedges of which the diamond is composed 
(see Fig. 5). This drag is given by: 


D' Ap f'(x)dx = —26p dx + 2p dx 
Dividing by '/2 pa*l: 
Cp = —(46/la®) dx + (46/la®) (37) 


Now ¢,—1.e. c,-—-can be obtained as the sum of the con- 
tributions to ¢; of two infinite wedges, wedge (1) start- 
ingat x = 0 and having vertex angle 6, and wedge (I1) 
starting at x = //2 and having vertex angle —26 (see 
Fig. 5). For 0 < x < 1/2, only wedge (I) contributes. 

Thus for 0 < vx < //2: = ¢, due to wedge (I), so 
that: 


4"? 
dx = — 
la? So 3 
2090/2’ K 
here, as well as in the following 8 = —2bl/a*, where / 


is the chord length of the diamond. 
For 1/2 < x < 1, = due to wedge (I) + ¢,, due 
to wedge (II), so that 


Site dx for wedge (I) + dx for 
wedge (IT) 


Substituting for (¢,); and (¢;);7: 
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‘ 2E — RK 
— (.214) ( + (39) 


Adding Eqs. (38) and (39): 


Cp = 1 + 2.76 
RK 


(40) 


Eqs. (40) and (21) have been used to plot cp vs. ./, 
at 8 = 0.02 and 0.1, for diamond profiles. Assuming 
that the greatest accelerations, which are nowadays 
available with rockets or otherwise, are of the order of 
5.10%. ..104 ft. per sec.,? the wings (e.g., / ~ 3 ft.) of 
such rockets will experience 6-values of the indicated 
order. It is seen from Fig. 7 that also for diamond- 
shaped profiles, (cp/6?) has a marked maximum at = 
1, but passes continuously through the transonic range, 
according to these considerations. It will be shown in 
another paper that in the subsonic range the drag falls 
off with decreasing .1/, until at very low speeds it reaches 
a value which is equal to the product of deceleration 
times “‘apparent mass.” 

Naturally the (c,/a) vs. \/-curve, plotted in Fig, 3, 
holds also for diamond profiles, as the lift coefficient of 
symmetrical airfoils does not depend on the profile 
form. 

Fig. 8 shows how cp varies with the acceleration 
parameter 6, at J = 1. 


(3) Concluding Remarks 


Extension to Finite Aspect Ratios. A Validity Cri- 
terion.—At the end of these considerations, an attempt 
shall be made to find out to what extent these results 
may be expected to give a quantitative account of 
actual conditions. Before this is done, it may be 
appropriate to mention an extension of the present re- 
sults to wings of finite aspect ratio and to discuss the 
results together. 

The aerodynamic characteristics for rectangular 
wings of finite aspect ratio have been derived by the 
authors* by using the method of “retarded potentials” 
in an x-y-t-space. While the derivations of this investi- 
gation must be left for a separate paper, some analytical 
results have been plotted in Fig. 9, representing Cp 
vs. \/ for different aspect ratios, at an acceleration 
parameter 8 = 0.02. The maximum of the curves at 
M = 1 is seen to become less and less pronounced, as 
the aspect ratio decreases. This is further elucidated 
by Fig. 10, representing cp vs. aspect ratio, at J = I, 
for different B-values; it is seen that cp — 0 for A.R. 
— 0, so long as 8 # 0. For 6 = 0, cp becomes loga- 
rithmically infinite for all finite A.R. # 0. 


* See a forthcoming report of the authors to O.N.R. 
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Fic. 10. cp vs. A.R., at M = 1, for different B-values. 


Since the results of the present investigation are 


obtained by a linearized approximation, they may be 
considered valid so long as the (velocity or) pressure- 
disturbance produced by the airfoil is small compared 
to the undisturbed (velocity or) pressure; in the pres- 
ent case the undisturbed quantities are (the sonic flight 
speed or) the dynamic pressure at the leading edge. 
In other words, it should be (Ap/'/spa?) = cp < 1. 

This criterion cannot be expected to be fulfilled at 
any station on the airfoil surface, separately. (At the 
leading edge, the pressure has been shown to grow be- 
yond bound at J/ = 1, similarly as it does in the thin 
airfoil theory of the subsonic regime.) However, the 
criterion should be indicative, if an average of the 
pressure distribution over the chord length is taken or, 
in the three-dimensional case, an average over the plan 
form. Therefore, a more a form of the cri- 
terion is 

op <1 
where ¢, = (1/1) %(x) dx or & = (1/A) Ji dx 
dy, respectively. But according to definition (35c), 
€p = Cp/5, which is shown in Fig. 9 for A.R. = 2, 7, and 
. 

For instance, if a very thin supersonic airfoil of 5 
per cent thickness ratio is taken, the curves in Fig. 9 
indicate at M = 1 for 8 = 0.02: ¢c, ~0.25 for A.R. ~; 
~0.2 for A.R. 7; and ~0.14 for A.R. 2. 

Thus, the above results might be considered as rea- 
sonably valid, even at 1J = 1. The smaller the aspect 
ratio is, the better is the approximation. Or, to put it 
more carefully: if quantitative results are desired, it 
might be worth while to develop a second order theory— 
e.g., along lines suggested by the paper of Lin, Reissner, 
and Tsien.* 


APPENDIX (1)—PROOF OF RELATION (12a) 
If «is a sufficiently small quantity, then: 
1 ak of (= 
— —— — 1 
V(M + a)?—1 = — 1) 


On the other hand: 


JOURNAL OF THE AERONAUTICAL SCIENCES—JANUARY, 1950 


1 
V(M+a)?—1 


id 


By means of the formula 


(l— 22h DY 


where Pk is the Legendre polynomial of degree } 
Eq. (2) can be expanded as follows: 


- 


l 


at ) 
V(M + a)? 1 


Comparing coeflicients in Eqs. (1) and (3): 


of ( ) 


P 


APPENDIX (2)-——RADIUS OF CONVERGENCE OF THE 
SERIES IN POWERS OF 8, Eas. (17) To (24) 


If in Eq. (30) the new variable of integration y is 
introduced, defined by y = ('/28 — Mz — !/227)/s,a 


straightforward calculation shows that 


J== 


bo 


ay?) VM + y)* +8 
vi — 


(My + 
dy 


and 


+1 
L= rM—2 fay + 


V1 
V(M +3 y)? + 


may be expanded eed a power series in 8 v alid pis: 
the range of integration. The power series for J and L 
and therefore the power series for the drag may then be 
obtained by means of term-by-term integration. Thus 
it is seen that the power series, as given in Section (A, 2), 
must converge for |a| < (M — 1)?. 
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APPENDIX (3) 
Computation of the quantities appearing in Eq. (32) 
yields: 
R, Rp = — 1)? + BIC + 1)? + x 
{Ar-2 + (WM + 8 (M+ 1)» + 
(M+ 1)? — 2MV(M + + B- 
+ 1)? + — 1)? + BIC-2} 


5, — Se = — 1)? + BIC +1)? 


where: 
— —2M* + 101/68 + 2M 
6 6 
—M?+ 28+4 
Ci-2 = 
6 
Furthermore: 


(M?-—1)+8 
2V — 1)? + + 1)? 4+ 8] 
2+ V(M+1)?+8 
2V(M+1)2+8 


=~ 
to 
| 


n = cotg@ = 
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Supersonic Airfoils Simplified 


James E. Broadwell 

Department of Aeronautical Engineering, University of Michigan, 
Ann Arbor, Mich. 

October 25, 1949 


EFERENCE IS MADE to the simplified ‘‘Approach to Supersonic 

Airfoil Theory” of W. E. Strohmeyer and D. R. Gero in the 
August issue, as corrected in the October issue (JOURNAL OF THE 
AERONAUTICAL SCIENCES). The following simpler and still more 
direct approach has been found particularly well suited to under- 
graduate instruction because of the easy visualization of the 
geometrical concepts involved. The customary symbols are 
employed, as noted in the accompanying sketch. 

It is assumed as usual that the general flow conditions at the 
given airfoil section are known and two-dimensional in character. 
The figure shows a simple example at M = V/2; a thin, flat 
plate moving horizontally through the air at a small angle of 
attack a. Consider what has happened since the LE was one 
chord length back of where the 7E is now. At that moment, 
pressure impulses traveled out from the LE’ in all directions, 
relative to the air at a rate of V/.1/ ft. per sec.; and, in the time 
required for the wing to go a distance 2c, will have reached the 
circumference of a circle of radius ¢ 4/2. Similar circles could be 
drawn from all other points through which the LE has passed, 
defining a zone limited by a common tangent line (the so-called 
Mach line) in this case at 45° to the direction of motion, as shown. 

Consider the air now as in motion relative to the airfoil’s lower 
surface. By construction, the pressure and deflection of air (be- 
cause of its reaction on the airfoil surface) are confined to the re- 
gion above and to the right of the line LE-A extended, because 
only within that zone is it possible for a disturbance to be trans- 
mitted. A particle of air striking directly under the LE will be 
deflected out of its original path through the angle a, but will 
have no immediate effect on particles still farther below. Parti- 
cles of air which have gone before, however, have made it possible 
for all the air in the triangle LE-TE-A to be similarly deflected 
and, for negligible changes in density, it is uniformly so deflected. 
Furthermore, due to the weak oblique shock wave and small de- 
flection, the longitudinal component of speed in this triangle is 
still approximately V. 

The total section of air (normal to the direction of motion) 
acted upon over a unit of span is then TE-A or approximately c; 
the mass of air acted upon is p Vc slugs per sec., and the downward 
velocity imparted to it is Va. The rate of change of momentum 
per sec., constituting the normal force for the lower surface (Ibs. 
per ft. of span), is 


F, = pV*ca = 2qca (1) 


The surface pressure per unit area is then Fj/c, and the gage 
Pressure coefficient 


b/q = Fi/cq = 2a (2) 
Similar analysis holds for the upper surface, the pressure there 

being numerically equal but negative. Thus the section lift 

coefficient (substantially equal to the normal force coefficient) : 


Ci = (Fi — Fu/cq) = 4a (3) 


It must be noted that Eqs. (2) and (3) are valid for the assumed 
M = ~/2 only. However, it is clear from the same principles 
that in any case the length of the line 7E-A, and hence the react- 
ing force, will be proportional to tan yw, which in this case happened 
to be unity. As tan» = 1/ V M? — 1, the final equations for 
pressure and lift (for a< yu) become: 


p/q = (2a/ VM? — 1) (4) 
CG = (4a/ VM? — 1) (5) 


In Eq. (4), @ is considered negative when in a direction to pro- 
duce negative p. The corresponding drag (neglecting surface 
friction) is then simply the rearward component of the normal 
force, or approximately in this case: 


Cig = Cra = (4a*/ — 1) (6) 


Here the subscript w indicates that it is wave drag, because the 
actual energy loss making up this element of the drag occurs in 
two oblique shock waves, one below the leading edge (practically 
coincident with the line LEZ-A) and a similar one above the trail- 


ing edge. 


The same principles and the same result, expressed in Eq. (4) 
hold for any small increment of surface in two-dimensional flow. 
Hence, aerodynamic forces on finitely thick or curved surfaces 
complying with these general conditions can be found by simple 
summation or integration of the pressure times each increment of 
area involved. For a cambered airfoil of finite thickness with 
sharp LE and TE, C; integrates to the same form as Eq. (5) if @ 
is taken as the angle of attack of a straight line joining the LE 
and TE. 

Corrections for sweep or tip loss are negligible for small LE 
sweep angle (up to about 10°) and for tips cut back at approxi- 
mately the Mach angle. The theory becomes increasingly in- 
accurate at / — 1, and in any case where the distance from the 


given airfoil section to the tip is less than c/ V M? — 1. Subject 


to these and other minor qualifications, there is no induced angle 
in the sense employed in low-speed airfoil theory, and no direct 
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effect of aspect ratio. Thus wave drag pertains to thickness as _ bers (predicting a large reduction in damping) cited by us jy rel. ‘Mil 
well as to lift (angle-of-attack), and can be superimposed from _ erence 3 are far from realistic and that Statler’s use of strip theory | sonic S 


results of separate computations. From test data thus far avail- cannot be a valid basis for any general conclusion. Nevertheles | 24.4 

able, frictional drag appears to be of the same order of magnitude an investigation of a conventional, subsonic configuration Oat, “7 
as under subsonic conditions. low aspect ratio tail aft of a straight wing based on the results 9 6 Jon 


references 3 and 7, reveals that the only important unsteady floy | 1.8. N 


4 effect appears in the “complex downwash lag” at the tail, jp 

agreement with Statler’s conclusion. (This investigation wa ae 

carried out as a part of a more general study being conducted g | yo, 18 
Unsteady Flow Theory in Dynamic Stability the North American Aerophysics Laboratory.) However, thi 
conclusion cannot be extended to such configurations as th 
John W. Miles canard or the flying wing or to any configuration in supersonic 


Department of Engineering, University of California at Los Angeles flow without further analysis, and, in particular, it appears thy 
August 22, 1949 unsteady flow considerations will exercise a profound effect 


the pitch damping of a flying wing in either subsonic or supe. 


N A PAPER delivered at a recent I.A.S. meeting,! I. C. Statler th 
I concluded that ‘‘unsteady flow effects are not important in complex sidewash correction, may be expected to exercise a pro 
the calculation of the dynamic stability and response of high-speed vertical p. P. 
With this we agree, and we fin, as implied by the tendency of many high-speed, jet aircraft 


(with fins designed on the basis of steady flow theory) to ‘snake’ es 


wish to point out (what we regard as) certain fallacies in both the > : : i 
Thus, while we should agree with Statler’s results for the par. 


interpretation of the experimental data and the theoretical : 
arguments presented in his paper. At the same time, we wish ticular configuration which he considered, we cannot agree with 
a ‘ mie his conclusion that this happy state of affairs can b icipated 

to modify some earlier predictions of our own.? ebeie 8 happy state of affairs can be anticipate ’ 
In calculating the dynamic stability derivatives, we may gener- 8 i [ \ | 


Septer 


As previously pointed out,” * one of the errors in applying 


ally assume that the reduced frequency k( = wc/2U) is less than ; s THE J 
0.1. It follows that we may neglect terms of order & or higher steady flow theory to unsteady flows bea the failure of the former subjec 
in the seredynamic forces in phase with displacements or ac- to include the term p(0%/d1) in the Euler equation for the pres ion 
celerations, since they are to be compared with aerodynamic S“e- This same objection must be made to Statler’s analysis of tel 
terms of order, zero, and the very large mechanical inertia terms, the oscillating fuselage. The potential — by the Munk the e1 
respectively. In the case of the aerodynamic forces in phase with theory® for a slender body of revolution is indeed valid for non high 

the velocities, however, we must consider all terms of order k steady flow, but the lift per unit length on a body of revolutionis ol 
(including terms of order k log k), and it follows that the first- ' OL (x, t) U2 ra) 1 0 \} S(x)w(x, 2) which 
order effects of unsteady flow considerations will manifest them- + U ot U subst: 
selves only in the damping terms. Moreover, in a system which d Int 
is sufficiently underdamped to exhibit marked resonances, damp- w(x, ind ¢ 
ing exerts only a second-order effect on the response of the system — "° (x) 1s the cross-sections area atx. (We hope toc iscuss t i Vain 
away from resonance, which will be governed by displacement or situation — extensively in a future paper.) The pene srove 
“spring” forces at lower frequencies and inertia forces at higher time differentiation on that part of w(x, ¢) in phase with the dis rea 
frequencies, whereas the magnitude of the resonant response is placement evidently to force. Th: 
critically dependent on damping. Hence, the effects of discrep- angular velocity, / the length, and So the vase area, the damping ilies 
ancies between the predictions of steady and unsteady flow oe due to an oscillation about a point a aft of the nose, hold | 
theory as to the dynamic response of a system must be inferred ound to be given by comp 
from its behavior in the vicinity of resonance. Examination shear 


2 
dM/[d(ql/U)] = (1 *) 


of Statler’s results reveals that this is indeed the case, albeit he 
draws his conclusions from the behavior at the higher frequencies. 

In two recent discussions in tnese columns, we have presented 
the first-order (in frequency) results for two-dimensional, thin 
airfoils in subsonic flow.” * (The first-order results for super- 
sonic flow had been given earlier, for both the two- and three- 
dimensional cases.” ®) These results indicate that the damping 
moment due to rotary motion (i.e., the total component of mo- 
ment in phase with the angular velocity) is of order k log k for 
sufficiently small k, so that the pitch-damping derivative be- 
comes infinite as log k when k approaches zero, indicating two- 
dimensional, unsteady flow considerations to be increasingly 
important (in determining the damping per cycle) as k is decreased, 
at least if strip theory is used, as in Statler’s calculations. More- 
over, the results of reference 3 indicate that this logarithmic term 


The result obtained by Statler (which, incidentally, retains some — 
terms not justified by slender body theory) is appropriate to flight wich 
along a curved flight path at constant angle of attack, but not te Spit 
the case of oscillation about a straight flight path. Indeed, itis 2). 
the confusion of these two situations which has been primarily % on 
responsible for the failure to account correctly for the effects 5; fol 
of unsteady flow in the calculation of dynamic stability deriva: 
tives 

We do not wish to appear unduly critical of Statler’s paper ané 
remark that the errors we have discussed are extant in much of With 
the antecedent work on the same subject, There has been 4 
good deal of specious reasoning relative to nonstationary process J 3. . 
in moving reference frames; in many cases this may not have 
led to serious error, but it may be recalled that Zeno’s sophistry 


requires a compressibility correction factor of (1 — 
rather than the usual (1 — \/?)~'/?, Unfortunately, as shown with respect to time-dependent phenomena did not, in fact, pre- 
by Jones® and Reissner,’ strip theory is completely inadequate om Ge the while 
for the treatment of unsteady flow effects at extremely low fre- catching Achilles. 
quencies. Indeed, the reduction of Reissner’s results for small R 
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must be subtracted from the term k& log k in the two-dimensional ' Statler, I. C., Dynamic Stability at High Speeds from Unsteady Flot | the 
result, whence the term log & in the pitch-damping derivative is | /heory, New York, N.Y. (January 24-27, 1949); L.A.S. Preprint 187. 
replaced by the considerably smaller (in magnitude) term log 2 Miles, J. W., Quasi-Stationary Thin Airfoil Theory, Journal of the Acro 

re ‘a ° : nautical Sciences, Vol. 16, No. 7, p. 440, July, 1949. 4 
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Stability of Sandwich Plates in Combined 
Shear and Compression 


P. Bijlaard 
Professor, Institute ot Technology, Delft, Holland, on leave; 
Associate Professor, Cornell University, Ithaca, N.Y. 


September 23, 1949 


— HE HAD written his Reader’s Forum item published 
on pages 573-574 of the September issue of the JOURNAL OF 
THE AERONAUTICAL SCIENCES, two other papers on the same 
subject which substantiate his conclusions came to the writer’s 
knowledge. In references 1 and 2 the shear buckling of infi- 
nitely long, simply supported sandwich plates is investigated by 
the energy method. Though this method gives somewhat too 
high critical stresses, these authors find a close agreement be 
tween their values and those of the present writer’s method, 
which is either exact or conservative, so that both solutions are 
substantially accurate. 

In reference 1 interaction curves for combined loading in shear 
ud compression were derived from an earlier paper by van 
Vijngaarden. For shear and longitudinal compression they 
yrove to be approximated by a parabola, which was known 
iuready to obtain for homogeneous plates. 

That the parabola also holds true for sandwich plates is demon- 
strated directly by the writer’s method, and it can be shown to 
hold exactly for case (2). With a ratio P;/S = y between the 
compression ?, in longitudinal (X-) direction of the plate and the 
shear S, both per unit width, the critical shear is 

in which .Sp is the critical shear of the faces; S; obtains for a sand- 
wich plate with infinite transverse shearing rigidity (case 1) and 
5: if only the shear modulus G, of the core has a finite value (case 
2). All these values refer to the same ratio y = P,/S. Values 
Sjand S; follow from the theory of homogeneous plates. Value 
5; follows from the differential equation for case (2), being here* 


[(he + ty)?Ge/he]V — Pre + 25S, 


» - = () (2) 
ox? Oxoy 


With w = Ye'** and Y = e'®¥ Eq. (2) yields 


| h, All 
(he +t)?Ge +4)? G 


he Py 1/3 
| (3) 


(he + Ge 
while 
Y = (4) 
The boundary conditions along the long sides y = +)/2 lead to 
the condition 
(B2b/2) = + nw (5) 


With « = x/L, in which L is the half-wave length in X-direction, 
and with » = 1, Eq. (5) yields 
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(he + ty)? Ge (he + 


which, according to Eqs. (33). p. 86, and (15), p 189 of reference 
5, may be written as 


(S2/Seer,)? = 1 — (Ps2/P reer.) (7) 


in which So,,, and Pz, obtain for shear alone or compression 
alone. Hence, for case (2) the parabolic interaction curve holds 
exactly, if buckling stresses for equal wave lengths are considered. 
The writer found that in the latter case this relation holds also 
sufficiently accurately for homogeneous plates (case 1). Hence, 
with Eq. (1) it may be easily shown that it holds also for the 
complete sandwich plate. 

From Eq. (6) it follows that the critical shear S.,, may also 
be calculated directly for combined shear and longitudinal com- 
pression with Eqs. (5), (6), and (7) of the writer’s Reader’s Forum 
item in the September issue, if in these equations value 


a = —(y/2) + [(y?/4) + 1 + (8) 


In those Eqs. (5), (6), and (7) value 7 is the ratio between plastic 
and elastic buckling stress for the assumed wave length. 

In the Reader’s Forum item in the September issue on page 
574, a typographical error occurs in the last formula of Eq. (7), 
where the middle division bar should vanish, so as to read: 


r = (w2D,/b°G.) (le + ty)®] 


Furthermore, in the last sentence of that item, a line was omitted 
in typing. The writer wanted to say there that the expression 
(D;/D, + ¢) in Eq. (6) is involved in the minimization of P.,. or 
S.,. by variation of the wave length, so that zero wave lengths, as 
found with vanishing D;, do not occur here. 
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Errata—Note on the Velocity of Sound 


Morris Morduchow 
Polytechnic Institute of Brooklyn, Brooklyn, N.Y. 
June 25, 1949 


AM OF ERRORS in the article ““Note on the Velocity of 
Sound,” by Morris Morduchow (JouRNAL OF THE AERO- 
NAUTICAL SCIENCES, Vol. 16, No. 10, p. 635, October, 1949) fol- 
lows: 

Page 635, line 21, the sentence ‘The derivative dp/dp is then 
zero everywhere except at the discontinuity where its value is 
indeterminate,” should be replaced by the sentence “The deriva- 
tive dp/dp is then indeterminate everywhere including the dis- 
continuity.” 

The following misprints may also be noted. On line 28, the 
word “‘of” should be omitted between the words ‘‘limit’’ and ‘‘the 
equations.” On the top of p. 635, right column, the equation 
should read: 


udu = (1/2)yd(p/p) 
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On the Reduction of Unsteady Supersonic 
Flow Problems to Steady Flow Problems 


John W. Miles 
Department of Engineering, University of California at Los Angeles 


September 26, 1949 


W' CONSIDER a flow of supersonic velocity LU, directed along 
the positive x-axis, over an airfoil whose projection on the 
plane z = 0 occupies the region S, and which has the prescribed 
normal velocity w(x, y) exp (iwf). The linearized, boundary value 
problem is then represented by the equations 


1V2 — M?[(0/dx) + ik]*} o(x, y, 2) = 0 (1) 


g(x, y,0 +) = w(x, y), (x, y) in Sy (2) 


o(x, y,0 +) = 0 +) exp[—tk(x — x)], 
(x, y) in Sp (3) 
(x, y, O) = 0, (x, y) in S3 (4) 
k = whb/U (5) 
where M is the Mach Number, & the reduced frequency, ) the 
unit of length, x = x,(y) the trailing edge, S. the wake aft of the 
airfoil, including the trailing edge, and S; the remainder of the 
plane z = 0 (i.e., exclusive of S;, S:). The potential implicitly 

possesses the harmonic time dependence exp(/wf). 

We now restrict our problem by assuming the trailing edge to be 
everywhere supersonic, so that the zone of action of the wake has 
no region in common with the zone of influence of the wing. 
The boundary condition (2) may then be suppressed.  Intro- 
ducing the transformation! 

x’ = x,y’ = By, 2 = Bz, b’ = (6) 
(x’, y’, 2’) = exp Mx) $(x, y, 2) (7) 
w’(x’, y’) = B- exp (ix Mx)w(x, y) (8) 
= B-2Mk (9) 

the boundary value problem of Eqs. (1) to (4) is reduced to 
(10) 
(11) 
(12) 


(x', y’, O +) = w(x’, y’), (x’y’) in Si’ 


'(x’, 0) = 0, y’) in S,’ 


where S,’ and S;’ differ from S; and S; only by the scale transfor- 
mation of Eq. (6). 

We now consider the simplified boundary value problem of 
Eqs. (10) to (12) and drop the primes on the coordinates. Mag- 
naradze and Galin® have shown that, if ¥(x, y, s) is a solution to 
the steady flow equation 

= 0 (13) 


a solution to Eq. (10) is given by 


x 
2) = ¥(x, 2) — 2) (0/0E) X 


— (14) 


where x = 0 is taken at the most forward point on the leading 
edge, and Jo is Bessel’s function of order zero. The derivation 
given by Galin is ingenious but does not clearly exhibit the re- 
strictions on the result. However, having the result, we may 
verify it by direct substitution and state the sufficient (and, for 
all practical cases, necessary) conditions as: 


1950 


(a) ¢ and its first derivatives with respect to x, y, z must be 
continuous in the open domain z > 0. 

y, 2) = be(%, 2) = 3) = ¥, 2) = 0, whem 
x = X(y, 2) is the upstream boundary of the zone of action of thy 
airfoil (e.g., = + (y? + 3*)'/? for a delta wing with subsonie 
leading edges). 

In view of the condition (b), we may rewrite Eq. (14) in the 
form 


x / 
o(x, y, 3) = Jo «x(x? — €2)' y, 2) (15) 


In general, the conditions (a) and (b) will not be satisfied by 
the velocity potential, but this difficulty may be circumvented 
by recasting the boundary value problem in terms of the modi 
fied potential 


x 
P(x, y, 2) = (16) 


as done by Galin. In the case where the leading edge is every. 
where subsonic, it does not appear that this would be necessary, 
In order to apply the transformation (15), we require its in. 


version. Assuming the conditions (a) and (6), we obtain 


¥(x, y, 5) = — 
2) 

where Jo is Bessel’s function of imaginary argument. This re 

sult may be established by direct substitution in Eqs. (10) and 

(13). 

We may use Eq. (17) to transform the boundary condition (1), 
or, if is introduced, the integral of Eq. (11), in accordance with 
Eq. (16). We then have the problem of finding a solution to 
Eq. (13) subject to a prescribed, normal (z) derivative over §, 
and vanishing over S;. Since both the boundary conditions 
and the potential equation are those for steady flow, albeit the 
prescribed, normal velocity may exhibit a frequency dependenee, 
our problem is reduced to one in steady flow. The restriction 
to a supersonic trailing edge is necessary for this reduction to am 
equivalent, steady flow problem, since the boundary condition 
(3) would not transform to the analogous steady flow boundary 
condition, 


y,0 +) = +y(x, y, +) in Sy 


Galin? has applied the foregoing transformation to the oseik 
lating, rectangular wing, but the results appear to be consider 
ably more complex than those obtainable by other methods? 
The principal difficulty in its application to the usual problems 
is the evaluation of the required integrals, although the results 
may be expanded in terms of Sonine’s integrals.‘ Following this 
pattern we have checked the method by deriving the two-dimet 
sional result. We note that a related method has been proposed 
by Germain and Bader in a paper not yet available to us.® 
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(7) above. 
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